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Abstract 

Fundamental constants are a cornerstone of our physical laws. Any constant varying in 
space and/or time would reflect the existence of an almost massless field that couples to mat- 
ter. This will induce a violation of the universality of free fall. It is thus of utmost importance 
for our understanding of gravity and of the domain of validity of general relativity to test for 
their constancy. We thus detail the relations between the constants, the tests of the local posi- 
tion invariance and of the universality of free fall. We then review the main experimental and 
observational constraints that have been obtained from atomic clocks, the Oklo phenomenon, 
Solar system observations, meteorites dating, quasar absorption spectra, stellar physics, pul- 
sar timing, the cosmic microwave background and big bang nucleosynthesis. At each step we 
describe the basics of each system, its dependence with respect to the constants, the known 
systematic effects and the most recent constraints that have been obtained. We then describe 
the main theoretical frameworks in which the low-energy constants may actually be varying 
and we focus on the unification mechanisms and the relations between the variation of differ- 
ent constants. To finish, we discuss the more speculative possibility of understanding their 
numerical values and the apparent fine-tuning that they confront us with. 
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1 Introduction 



Fundamental constants appear everywhere in the mathematical laws we use to describe the phe- 
nomena of Nature. They seem to contain some truth about the properties of the physical world 
while their real nature seem to evade us. 

The question of the constancy of the constants of physics was probably first addressed by 
Dirac |154l I155j who expressed, in his "Large Numbers hypothesis" , the opinion that very large 
(or small) dimensionless universal constants cannot be pure mathematical numbers and must not 
occur in the basic laws of physics. He suggested, on the basis of this numerological principle, that 
these large numbers should rather be considered as variable parameters characterizing the state of 
the universe. Dirac formed five dimensionless ratios among which S = H^h/m-pC^ ~ 2/i x 10^*^ 
and e = Gpq/Hq ~ 5/i~^ x 10"'* and asked the question of which of these ratio is constant as the 
universe evolves. Usually, S and e vary as the inverse of the cosmic time. Dirac then noticed that 

/ /ictj.^ , representing the relative magnitude of electrostatic and gravitational forces between a 
proton and an electron, was of the same order as H^e^ / rricC^ = 5a^^^^ representing the age of the 
universe in atomic units so that his five numbers can be "harmonized" if one assumes that and 
5 vary with time and scale as the inverse of the cosmic time. 

This argument by Dirac is indeed not a physical theory but it opened many doors in the 
investigation on physical constants, both on questioning whether they are actually constant and 
on trying to understand the numerical values we measure. 

First, the implementation of Dirac's phenomenological idea into a field-theory framework was 
proposed by Jordan [267 who realized that the constants have to become dynamical fields and 
proposed a theory where both the gravitational and fine-structure constants can vary (Ref. |492j 
provides some summary of some earlier attempts to quantify the cosmological implications of 
Dirac argument). Fierz |196j then realized that in such a case, atomic spectra will be spacetime- 
dependent so that these theories can be observationally tested. Restricting to the sub-case in which 
only G can vary led to definition of the class of scalar-tensor theories which were further explored 
by Brans and Dicke |66j . This kind of theory was further generalized to obtain various functional 
dependencies for G in the formalization of scalar-tensor theories of gravitation (see e.g. Ref. [123) 1. 

Second, Dicke [149] pointed out that in fact the density of the universe is determined by its age, 
this age being related to the time needed to form galaxies, stars, heavy nuclei... This led him to 
formulate that the presence of an observer in the universe places constraints on the physical laws 
that can be observed. In fact, what is meant by observer is the existence of (highly?) organized 
systems and this principle can be seen as a rephrasing of the question "why is the universe the way 
it is?" (see Ref. [251] ). Carter [8T1|82], who actually coined the term "anthropic principle" for it, 
showed that the numerological coincidences found by Dirac can be derived from physical models 
of stars and the competition between the weakness of gravity with respect to nuclear fusion. Carr 
and Rees [79] then showed how one can scale up from atomic to cosmological scales only by using 
combinations of a-^„, and rric/mp. 

To summarize, Dirac's insight was to question whether some numerical coincidences between 
very large numbers, that cannot be themselves explained by the theory in which they appear, was 
a mere coincidence or whether it can reveal the existence of some new physical laws. This gives 
three main roads of investigations 

• how do we construct theories in which what were thought to be constants are in fact dynamical 
fields, 

• how can we constrain, experimentally or observationally, the spacetime dependencies of the 
constants that appear in our physical laws 

• how can we explain the values of the fundamental constants and the fine-tuning that seems 
to exist between their numerical values. 
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While "Varying constants" may seem, at first glance, to be an oxymoron, it has to be considered 
merely as jargon to be understood as "revealing new degrees of freedom, and their coupling to 
the known fields of our theory" . The tests on the constancy of the fundamental constants are 
indeed very important tests of fundamental physics and of the laws of Nature we are currently 
using. Detecting any such variation will indicate the need for new physical degrees of freedom in 
our theories, that is new physics. 

The necessity of theoretical physics on deriving bounds on their variation is, at least, threefold: 

1. it is necessary to understand and to model the physical systems used to set the constraints. 
In particular one needs to determine the effective parameters that can be observationally 
constrained to a set of fundamental constants; 

2. it is necessary to relate and compare different constraints that are obtained at different 
spacetime positions. This often requires a spacetime dynamics and thus to specify a model 
as well as a cosmology; 

3. it is necessary to relate the variations of different fundamental constants. 

We shall thus start in § [2] by recalling the link between the constants of physics and the theories 
in which they appear, as well as with metrology. From a theoretical point of view, the constancy 
of the fundamental constants is deeply linked with the equivalence principle and general relativity. 
In § [5] we will recall this relation and in particular the link with the universality of free fall. 
We will then summarize the various constraints that exist on such variations, mainly for the fine 
structure constant and for the gravitational constant in § |3] and |4] respectively. We will then turn 
to the theoretical implications in §[5] in describing some of the arguments backing up the fact that 
constants are expected to vary, the main frameworks used in the literature and the various ways 
proposed to explain why they have the values we observe today. We shall finish by a discussion 
on their spatial variations in § [6] and by discussing the possibility to understand their numerical 
values in § [71 

Various reviews have been written on this topic. We will refer to the review |495) as FVC and 
we mention the following later reviews [27l SHI [HI [ml |228l |27Ql [2731 EMI SMI [500] and we refer 
to Ref. [354] for the numerical values of the constants adopted in this review. 

2 Constants and fundamental physics 
2.1 About constants 

Our physical theories introduce various structures to describe the phenomena of Nature. They 
involve various fields, symmetries and constants. These structures are postulated in order to 
construct a mathematicaly consistent description of the known physical phenomena in the most 
unified and simple way. 

We define the fundamental constants of a physical theory as any parameter that cannot be 
explained by this theory. Indeed, we are often dealing with other constants that in principle can 
be expressed in terms of these fundamental constants. The existence of these two sets of constants 
is important and arises from two different considerations. From a theoretical point of view we 
would like to extract the minimal set of fundamental constants, but often these constants are not 
measurable. From a more practical point of view, we need to measure constants, or combinations 
of constants which allow to reach the highest accuracy. 

These fundamental constants are thus contingent quantities that can only be measured. Such 
parameters have to be assumed constant in this theoretical framework for two reasons: 
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• from a theoretical point of view: the considered framework does not provide any way to 
compute these parameters, i.e. it does not have any equation of evolution for them since 
otherwise it would be considered as a dynamical field, 

• from an experimental point of view: these parameters can only be measured. If the theories 
in which they appear have been validated experimentally, it means that, at the precisions of 
these experiments, these parameters have indeed been checked to be constant, as required 
by the necessity of the reproductibility of experimental results. 

This means that testing for the constancy of these parameters is a test of the theories in which 
they appear and allow to extend our knowledge of their domain of validity. This also explains the 
definition chosen by Weinberg |521| who stated that they cannot be calculated in terms of other 
constants "... not just because the calculation is too complicated (as for the viscosity of water) 
but because we do not know of anything more fundamental" . 

This has a series of implications. First, the list of fundamental constants to consider depends 
on our theories of physics and thus on time. Indeed, when introducing new, more unified or 
more fundamental, theories the number of constants may change so that this list reflects both our 
knowledge of physics and, more important, our ignorance. Second, it also implies that some of these 
fundamental constants can become dynamical quantities in a more general theoretical framework 
so that the tests of the constancy of the fundamental constants are tests of fundamental physics 
which can reveal that what was thought to be a fundamental constant is actually a field whose 
dynamics cannot be neglected. If such fondamental constants are actually dynamical fields it also 
means that the equations we are using are only approximations of other and more fundamental 
equations, in an adiabatic limit, and that an equation for the evolution of this new field has to be 
obtained. 

The reflections on the nature of the constants and their role in physics are numerous. We 
refer to the books [Ml EIH [Ml [502] as well as Refs. [1641 [383 MJ\ MH [534] for various 
discussions on this issue that we cannot develop at length here. This paragraph summarizes some 
of the properties of the fundamental constants that have attracted some attention. 

2.1.1 Characterizing the fundamental constants 

Physical constants seem to play a central role in our physical theories since, in particular, they 
determined the magnitudes of the physical processes. Let us sketch briefly some of their properties. 

How many fundamental constants should be considered? The set of constants which are conven- 
tionally considered as fundamental [214 consists of the electron charge e, the electron mass TOc, the 
proton mass , the reduced Planck constant h, the velocity of light in vacuum c, the Avogadro 
constant iV^^, the Boltzmann constant , the Newton constant G, the permeability and permittiv- 
ity of space, eo and hq. The latter has a fixed value in the SI system of unit (/io = Airx 10"'' H m"-'^) 
which is implicit in the definition of the Ampere; eo is then fixed by the relation eoMo — c^^. 

It is however clear that this cannot corresponds to the list of the fundamental constants, as 
defined earlier as the free parameters of the theoretical framework at hand. To define such a list 
we must specify this framework. Today, gravitation is described by general relativity, and the 
three other interactions and the matter fields are described by the standard model of particle 
physics. It follows that one has to consider 22 unknown constants (i.e. 19 unknown dimensionless 
parameters): the Newton constant G, 6 Yukawa couplings for the quarks (/iu, h^, hs, ht, h^) and 
3 for the leptons (/ic, h^, hr), 2 parameters of the Higgs field potential (/i. A), 4 parameters for the 
Cabibbo-Kobayashi-Maskawa matrix (3 angles 9ij and a phase (^ckm); 3 coupling constants for the 
gauge groups SU{3)c x SU{2)l x U{1)y {91,92,93 or equivalently 52,33 and the Weinberg angle 
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Constant 


Symbol 


Value 


Speed of light 


c 


299792458 m-s"' 


Planck constant (reduced) 


h 


1.05457f 628(53) x 10"^* J -s 


Newton constant 


G 


6.67428(67) x 10"" m^ • kg"^ • s^^ 


Weak coupling constant (at mz) 




0.6520 ± O.OOOf 


Strong coupling constant (at mz) 


gsi^z) 


1.221 ±0.022 


Weinberg angle 


sin^6'„(91.2GEV)„-, 


0.23120 ± 0.00015 


Electron Yukawa coupling 


he 


2.94 X 10"" 


Muon Yukawa coupling 


ha 


0.000607 


Tauon Yukawa coupling 


hr 


0.0102156 


Up Yukawa coupling 


hn 


0.000016 ± 0.000007 


Down Yukawa coupling 


hd 


0.00003 ± 0.00002 


Charm Yukawa coupling 


he 


0.0072 ± 0.0006 


Strange Yukawa coupling 


hs 


0.0006 ± 0.0002 


Top Yukawa coupling 


ht 


1.002 ± 0.029 


Bottom Yukawa coupling 


hh 


0.026 ± 0.003 


Quark CKM matrix angle 


sin 612 


0.2243 ±0.0016 




sin 623 


0.0413 ±0.0015 




sin 6*13 


0.0037 ± 0.0005 


Quark CKM matrix phase 


r 

^CKM 


1.05 ±0.24 


Higss potential quadratic coefficient 
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Higss potential quartic coefficient 


A 


? 


QCD vacuum phase 


^'qcd 


< 10-8 



Table 1: List of the fundamental constants of our standard model. The numerical values are given 
in the Planck system of units (see below) defined by the requirement that the numerical value ofG, 
c and h is 1 in this system of units. 



9w)j and a phase for the QCD vacuum (Oqcb), to which one must add the speed of light c and 
the Planck constant h. See Table [T] for a summary and their numerical values. 

Again, this list of fundamental constants relies on what we accept as a fundamental theory. 
Today we have many hints that the standard model of particle physics has to be extended, in 
particular to include the existence of massive neutrinos. Such an extension will introduce at 
least seven new constants (3 Yukawa couplings and 4 Maki-Nakagawa-Sakata (MNS) parameters, 
similar to the CKM parameters). On the other hand, the number of constants can decrease if some 
unifications between various interaction exist (see ii l5.3.1l for more details) since the various coupling 
constants may be related to a unique coupling constant au and an energy scale of unification Mu 
through 

where the hi are numbers which depend on the explicit model of unification. Note that this would 
also imply that the variations, if any, of various constants shall be correlated. 

Relation to other usual constants. These parameters of the standard model are related to 
various constants that will appear in this review (see Tabled]). First, the quartic and quadratic 
coefficients of the Higgs field potential are related to the Higgs mass and vev, tuh — \J —[J? and 

V = yJ—fi^/X. The latter is related to the Fermi constant — (w^-\/2)-^ which imposes that 

V = (246.7 ± 0.2) GeV while the Higgs mass is badly constrained. The masses of the quarks and 
leptons are related to their Yukawa coupling and the Higgs vev hy m — hv/ The values of the 
gauge couplings depend on energy via the renormalisation group so that they are given at a chosen 
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energy scale, here the mass of the Z-boson, mz- gi and 32 are related by the Weinberg angle as 
gi = 52 tan 0^. The electromagnetic coupling constant is not gi since SU{2)l x U{1)y is broken 
to C/(l)cicc so that it is given by 

5'em("^z) = e = 52(TOz)sin6l„. (1) 

Defining the fine-structure constant as a^.^, — g'^„/f>'C, the (usual) zero energy electromagnetic 
fine structure constant is = 1/137.03599911(46) is related to a^^^['mz) = 1/(127.918 ± 

0.018) by the renormalisation group equations. In particular, it implies that ^ a{mz) + 

;r- In ( —4 — 4 — 3 — 3 — 3- 1 . We define the QCD energy scale, A„„„, as the energy at which the 

strong coupling constant diverges. Note that it implies that A^^^^ also depends on the Higgs and 
fermion masses through threshold effects. 

More familiar constants, such as the masses of the proton and the neutron are, as we shall 
discuss in more details below (see ^ I5.3.2p . more diflBcult to relate to the fundamental parameters 
because they depend not only on the masses of the quarks but also on the electromagnetic and 
strong binding energies. 

Are some constants more fundamental? As pointed-out by Levy-Leblond [326) . all constants 
of physics do not play the same role, and some have a much deeper role than others. Following 
Ref. [326], we can define three classes of fundamental constants, class A being the class of the 
constants characteristic of a particular system, class B being the class of constants characteris- 
tic of a class of physical phenomena, and class C being the class of universal constants. Indeed, 
the status of a constant can change with time. For instance, the velocity of light was initially 
a class A constant (describing a property of light) which then became a class B constant when 
it was realized that it was related to electromagnetic phenomena and, to finish, it ended as a 
type C constant (it enters special relativity and is related to the notion of causality, whatever 
the physical phenomena). It has even become a much more fundamental constant since it now 
enters in the definition of the metre |408j (see Ref. |503j for a more detailed discussion). This 
has to be constrasted with the proposition of Ref. |534| to distinguish the standard model free 
parameters as the gauge and gravitational couplings (which are associated to internal and space- 
time curvatures) and the other parameters entering the accomodation of inertia in the Higgs sector. 

Relation with physical laws. Levy-Leblond |326' thus proposed to rank the constants in terms 
of their universality and he proposed that only three constants be considered to be of class C, 
namely G, h and c. He pointed out two important roles of these constants in the laws of physics. 
First, they act as concept synthetizer duving the process of our understanding of the laws of nature: 
contradictions between existing theories have often been resolved by introducing new concepts that 
are more general or more synthetic than older ones. Constants build bridges between quantities 
that were thought to be incommensurable and thus allow new concepts to emerge. For example c 
underpins the synthesis of space and time while the Planck constant allowed to related the concept 
of energy and frequency and the gravitational constant creates a link between matter and space- 
time. Second, it follows that this constants are related the domains of validity of these theories. 
For instance, as soon as velocity approaches c, relativistic effects become important, relativistic 
effects cannot be negligible. On the other hand, for speed much below c, Galilean kinematics is 
sufficient. Planck constant also acts as a referent, since if the action of a system greatly exceeds 
the value of that constant, classical mechanics will be appropriate to describe this system. While 
the place of c (related to the notion of causality) and h (related to the quantum) in this list are 
well argumented, the place of G remains debated since it is thought that it will have to be replaced 
by some mass scale. 
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Constant 

Electromagnetic coupling constant 

Higss mass 

Higss vev 

Fermi constant 

Mass of the 

Mass of the Z 

Fine structure constant 

Fine structure constant at mz 

Weak structure constant at mz 

Strong structure constant at mz 

Gravitational structure constant 

Electron mass 

Mu mass 

Tau mass 

Up quark mass 

Down quark mass 

Strange quark mass 

Charm quark mass 

Bottom quark mass 

Top quark mass 

QCD energy scale 

Mass of the proton 

Mass of the neutron 

proton-neutron mass difference 

proton-to-electron mass ratio 

electron-to-proton mass ratio 

d — u quark mean mass 

d — u quark mass difference 

proton gyromagnctic factor 

neutron gyromagnetic factor 

Rydberg constant 



Symbol 



: e = 52 sin 



mu 

V 

mw 
mz 

O'lLmimz) 
(mz) 

= Gmp/hc 



me 



m^ = 
mr = 

mu = 
md = 
ms = 
mc = 
mb = 
m.t = 
Aqcd 
m„ 



mn 

Qnp 

H = mp/me 
p = me/mp 
mq = (mu -I- md)/2 



5mq 
9p 

9n 
Roo 



■ ma 



Value 

0.313429 ± 0.000022 
> 100 GeV 
(246.7 ± 0.2) GeV 
1.16637(1) X lO^-^GeV-^ 
80.398 ± 0.025 GeV 
91.1876 ± 0.0021 GeV 
1/137.035999679(94) 
1/(127.918 ±0.018) 
0.03383 ± 0.00001 
0.1184 ± 0.0007 
~ 5.905 X 10"^^ 
510.998910 ± 0.000013 kcV 
105.658367 ± 0.000004 MeV 
1776.84 ± 0.17 MeV 
(1.5-3.3) MeV 
(3.5 - 6.0) MeV 



105: 



MeV 



i.27+;;;';i cev 

4:.20tl\ll GeV 

171.3 ±2.3 Gev 

(190 - 240) MeV 

938.272013 ± 0.000023 MeV 

939.565346 ± 0.000023 MeV 

1.2933321 ± 0.0000004 MeV 

1836.15 

1/1836.15 

(2.5 - 5.0) MeV 

(0.2 ~ 4.5) MeV 

5.586 

-3.826 

10973731.568527(73) m'^ 



Table 2: List of some related constants that appear in our discussions. 



Evolution. There are many ways the list of constants can change with our understanding of 
physics. First, new constants may appear when new systems or new physical laws arc discovered; 
this is for instance the case of the charge of the electron or more recently the gauge couplings of 
the nuclear interactions. A constant can also move from one class to a more universal class. An 
example is that of the electric charge, initially of class A (characteristic of the electron) which 
then became class B when it was understood that it characterizes the strength of the electromag- 
netic interaction. A constant can also disappear from the list, because it is either replaced by 
more hmdamental constants (e.g. the Earth acceleration due to gravity and the proportionality 
constant entering Kepler law both disappeared because they were "explained" in terms of the 
Newton constant and the mass of the Earth or the Sun) or because it can happen that a better 
understanding of physics teaches us that two hitherto distinct quantities have to be considered 
as a single phenomenon (e.g. the understanding by Joule that heat and work were two forms of 
energy led to the fact that the Joule constant, expressing the proportionality between work and 
heat, lost any physical meaning and became a simple conversion factor between units used in the 
measurement of heat (calories) and work (Joule)). Nowadays the calorie has fallen in disuse. In- 
deed demonstrating that a constant is varying will have direct implications on our list of constants. 
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In conclusion, the evolution of the number, status of the constants can teach us a lot about 
the evolution of the ideas and theories in physics since it reflects the birth of new concepts, their 
evolution and unification with other ones. 

2.1.2 Constants and metrology 

Since we cannot compute them in the theoretical framework in which they appear, it is a crucial 
property of the fundamental constants (but in fact of all the constants) that their value can be 
measured. The relation between constants and metrology is a huge subject to which we just draw 
the attention on some selected aspects. For more discussions, see Refs. [2721 1273] . 

The introduction of constants in physical laws is also closely related to the existence of systems 
of units. For instance, Newton's law states that the gravitational force between two masses is pro- 
portional to each mass and inversely proportional to the square of their separation. To transform 
the proportionality to an equality one requires the use of a quantity with dimension of m^-kg~^ -s"^ 
independent of the separation between the two bodies, of their mass, of their composition (equiva- 
lence principle) and on the position (local position invariance). With an other system of units the 
numerical value of this constant could have simply been anything. Indeed, the numerical value of 
any constant crucially depends on the definition of the system of units. 

Measuring constants. The determination of the laboratory value of constants relies mainly on 
the measurements of lengths, frequencies, times,... (see Ref. |409j for a treatise on the measurement 
of constants and Ref. |214| for a recent review). Hence, any question on the variation of constants 
is linked to the definition of the system of units and to the theory of measurement. The behavior 
of atomic matter is determined by the value of many constants. As a consequence, if e.g. the 
fine-structure constant is spacetime dependent, the comparison between several devices such as 
clocks and rulers will also be spacetime dependent. This dependence will also differ from one clock 
to another so that metrology becomes both device and spacetime dependent, a property that will 
actually be used to construct tests of the constancy of the constants. 

Indeed a measurement is always a comparison between two physical systems of the same di- 
mensions. This is thus a relative measurement which will give as result a pure number. This trivial 
statement is oversimplifying since in order to compare two similar quantities measured separately, 
one needs to perform a number of comparisons. In order to reduce this number of comparisons 
(and in particular to avoid creating every time a chain of comparisons), a certain set of them has 
been included in the definitions of units. Each units can then be seen as an abstract physical sys- 
tem, which has to be realised effectively in the laboratory, and to which another physical system 
is compared. A measurement in terms of these units is usually called an absolute measurement. 
Most fundamental constants are related to microscopic physics and their numerical values can be 
obtained either from a pure microscopic comparison (as is e.g. the case for me/mp) or from a 
comparison between microscopic and macroscopic values (for instance to deduce the value of the 
mass of the electron in kilogram). This shows that the choice of the units has an impact on the 
accuracy of the measurement since the pure microscopic comparisons are in general more accurate 
than those involving macroscopic physics. 

It is also important to stress that in order to deduce the value of constants from an experiment, 
one usually needs to use theories and models. An example |273j is provided by the Rydberg 
constant. It can easily be expressed in terms of some fundamental constants as Roo — a^^mcc/2h. 
It can be measured from e.g. the triplet Is — 2s transition in hydrogen, the frequency of which 
is related to the Rydberg constant and other constants by assuming QED so that the accuracy 
of i?oo is much lower than that of the measurement of the transition. This could be solved by 
defining Roo as 4i^H(ls — 2s) /3c but then the relation with more fundamental constants would be 



10 



more complicated and actually not exactly known. This illustrates the relation between a practical 
and a fundamental approach and the limitation arising from the fact that we often cannot both 
exactly calculate and directly measure some quantity. Note also that some theoretical properties 
are plugged in the determination of the constants. 

As a conclusion, let us recall that (i) in general, the values of the constants are not determined 
by a direct measurement but by a chain involving both theoretical and experimental steps, (ii) 
they depend on our theoretical understanding, (iii) the determination of a self-consistent set of 
values of the fundamental constants results from an adjustment to achieve the best match between 
theory and a defined set of experiments (which is important because we actually know that the 
theories are only good approximation and have a domain of validity) (iv) that the system of units 
plays a crucial role in the measurement chain, since for instance in atomic units, the mass of the 
electron could have been obtained directly from a mass ratio measurement (even more precise!) 
and (v) fortunately the test of the variability of the constants does not require a priori to have a 
high-precision value of the considered constants. 

System of units. One thus need to define a coherent system of units. This has a long, complex 
and interesting history that was driven by simplicity and universality but also by increasing stability 
and accuracy [251 1502J . 

Originally, the sizes of the human body were mostly used to measure the length of objects (e.g. 
the foot and the thumb gave feet and inches) and some of these units can seem surprising to us 
nowaday (e.g. the span was the measure of hand with fingers fully splayed, from the tip of the thumb 
to the tip of the little finger!). Similarly weights were related to what could be carried in the hand: 
the pound, the ounce, the dram. . . Needless to say that this system had a few disadvantages since 
each country, region has its own system (for instance in France there was more than 800 different 
units in use in 1789). The need to define a system of units based on natural standard led to several 
propositions to define a standard of length (e.g. the mille by Gabriel Mouton in 1670 defined as 
the length of one angular minute of a great circle on the Earth or the length of the pendulum 
that oscillates once a second by Jean Picard and Christiaan Huygens) . The real change happened 
during the French Revolution during which the idea of a universal and non anthropocentric system 
of units arose. In particular, the Assemblee adopted the principle of a uniform system of weights 
and measures on the 8th of May 1790 and, on March 1791 a decree (these texts are reprinted in 
Ref. |503] ) was voted, stating that a quarter of the terrestrial meridian would be the basis of the 
definition of the metre (from the Greek metron, as proposed by Borda) : a metre would henceforth 
be one ten millionth part of a quarter of the terrestrial meridian. Similarly the gram was defined 
as the mass of one cubic centimetre of distilled water (at a precise temperature and pressure) and 
the second was defined from the property that a mean Solar day must last 24 hours. 

To make a long story short, this led to the creation of the metric system and then of the 
signature of La convention du metre in 1875. Since then, the definition of the units have evolved 
significantly. First, the definition of the metre was related to more immutable systems than our 
planet which, as pointed out by Maxwell in 1870, was an arbitrary and inconstant reference. He 
then suggested that atoms may be such a universal reference. In 1960, the BIPM established a new 
definition of the metre as the length equal to 1650763 wavelengths, in a vacuum, of the transition 
line between the levels 2pio and 5c?5 of krypton-86. Similarly the rotation of the Earth was not 
so stable and it was proposed in 1927 by Andre Danjon to use the tropical year as a reference, 
as adopted in 1952. In 1967, the second was also related to an atomic transition, defined as the 
duration of 9162631770 periods of the transition between the two hyperfine levels of the ground 
state of caesium-133. To finish, it was decided in 1983, that the metre shall be defined by fixing 
the value of the speed of light to c = 299792458m- s^^ and we refer to Ref. [478] for an up to date 
description of the SI system. Today, the possibility to redefine the kilogram in terms of a fixed 
value of the Planck constant is under investigation |271j . 
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This summary illustrates that the system of units is a human product and all SI definitions are 
historically based on non-relativistic classical physics. The changes in the definition were driven 
by the will to use more stable and more fundamental quantities so that they closely follow the 
progress of physics. This system has been created for legal use and indeed the choice of units is 
not restricted to SI. 

SI systems and the number of basic units. The International System of Units defines seven 
basic units: the metre (m), second (s) and kilogram (kg), the Ampere (A), Kelvin (k), mole (mol) 
and candela (cd), from which one defines secondary units. While needed for pragmatic reasons, 
this system of units is unnecessarily complicated from the point of view of theoretical physics. In 
particular, the Kelvin, mole and candela are derived from the four other units since temperature 
is actually a measure of energy, the candela is expressed in terms of energy flux so that both can 
be expressed in mechanical units of length [L] , mass [M] and time [T] . The mole is merely a unit 
denoting numbers of particule and has no dimension. 

The status of the Ampere is interesting in itself. The discovery of the electric charge [Q] led 
to the introduction of a new units, the Coulomb. The Coulomb law describes the force between 
two charges as being proportional to the product of the two charges and to the inverse of the dis- 
tance squared. The dimension of the force being known as [MLT~^], this requires the introduction 
of a new constant Eq (which is only a conversion factor), with dimensions [Q^M~^L~'^T^] in the 
Coulomb law, and that needs to be measured. Another route could have been followed since the 
Coulomb law tells us that no new constant is actually needed if one uses [M-'^/^L'^/^T"-'^] as the 
dimension of the charge. In this system of units, known as Gaussian units, the numerical value of 
eo is 1. Hence the Coulomb can be expressed in terms of the mechanical units [L], [M] and [T], and 
so will the Ampere. This reduces the number of conversion factors, that need to be experimentally 
determined, but this choice of units assumes the validity of the Coulomb law so that keeping a 
separate unit for the charge may be a more robust attitude. 

Natural units. The previous discussion tends to show that all units can be expressed in terms 
of the three mechanical units. It follows, as realized by Johnstone-Stoney in 1874, that these three 
basic units can be defined in terms of 3 independent constants. He proposed [25l 1266] to use three 
constants: the Newton constant, the velocity of light and the basic units of electricity, i.e. the 
electron charge, in order to define, from dimensional analysis a "natural series of physical units" 
defined as 



where the Eq factor has been included because we are using the SI definition of the electric charge. 
In such a system of units, by construction, the numerical value of G, e and c is 1, i.e. c = Ix -tg^ 
etc. 

A similar approach to the definition of the units was independently proposed by Planck [413] 
on the basis of the two constants a and b entering the Wien law and G, which he reformulated 





1.85 X 10"" ki 
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later |414j in terms of c, G and h as 




5.39056 X 10 



1.61605 X 10 



-35 



44 



m, 



mp 




- 2.17671 X 10"^ kg. 



The two systems are clearly related by the fine-structure constant since e'^/Aireo = a^j, /ic. 

Indeed, we can construct many such systems since the choice of the 3 constants is arbitrary. For 
instance, we can construct a system based on (e, mo, h), that we can call the Bohr units which will 
be suited to the study of the atom. The choice may be dictated by the system which is studied (it 
is indeed far fetched to introduce G in the construction of the units when studying atomic physics) 
so that the system is well adjusted in the sense that the numerical values of the computations are 
expected to be of order unity in these units. 

Such constructions are very useful for theoretical computations but not adapted to measure- 
ment so that one needs to switch back to SI miits. More important, this shows that, from a 
theoretical point of view, one can define the system of units from the laws of nature, which are 
supposed to be universal and immutable. 

Do we actually need 3 natural units? is an issue debated at length. For instance. Duff, Okun 
and Veneziano |164j respectively argue for none, three and two (see also Ref. |531j ). Arguing for no 
fundamental constant leads to consider them simply as conversion parameters. Some of them are, 
like the Boltzmann constant, but some others play a deeper role in the sense that when a physi- 
cal quantity becomes of the same order of this constant new phenomena appear, this is the case 
e.g. of h and c which are associated respectively to quantum and relativistic effects. Okun [386] 
considered that only three fundamental constants are necessary, as indicated by the International 
System of Units. In the framework of quantum field theory -I- general relativity, it seems that this 
set of three constants has to be considered and it allows to classify the physical theories (with the 
famous cube of physical theories). However, Veneziano '510' argued that in the framework of string 
theory one requires only two dimensionful fundamental constants, c and the string length A^. The 
use of h seems unnecessary since it combines with the string tension to give . In the case of the 
Goto-Nambu action S/h= (T/h) J d{Area) = Xj^ J d{Area) and the Planck constant is just given 
by Aj^. In this view, h has not disappeared but has been promoted to the role of a UV cut-off that 
removes both the infinities of quantum field theory and singularities of general relativity. This sit- 
uation is analogous to pure quantum gravity [3841 where H and G never appear separately but only 
in the combination i^^ = ^ Gh/c^ so that only c and £p, are needed. Volovik |516j made an analogy 
with quantum liquids to clarify this. There, an observer knows both the effective and microscopic 
physics so that he can judge whether the fundamental constants of the effective theory remain fun- 
damental constants of the microscopic theory. The status of a constant depends on the considered 
theory (effective or microscopic) and, more interestingly, on the observer measuring them, i.e. on 
whether this observer belongs to the world of low-energy quasi-particles or to the microscopic world. 

Fundamenal parameters. Once a set of three independent constants has been chosen as natural 
units, then all other constants are dimensionless quantities. The values of these combinations of 
constants does not depend on the way they are measured, |109[ I163| 1431) . on the definition of 
the units etc... It follows that any variation of constants that will let these numbers unaffected is 
actually just a redefinition of units. 

These dimensionless numbers represent e.g. the mass ratio, relative magnitude of strength etc... 
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Changing their values wih indeed have an impact on the intensity of various physical phenomena, 
so that they encode some properties of our world. They have specific values (e.g. a^^^ ~ 1/137, 
mp/nie ~ 1836, etc.) that we may hope to understand. Are all these numbers completely contin- 
gent, or are some (why not all?) of them related by relations arising from some yet unknown and 
more fundamental theories. In such theories, some of these parameters may actually be dynamical 
quantities and thus vary in space and time. These are our potential varying constants. 

2.2 The constancy of constants as a test of general relativity 

The previous paragaphs have yet emphasize why testing for the consistency of the constants is a 
test of fundamental physics since it can reveal the need for new physical degrees of freedom in our 
theory. We now want to stress the relation of this test with other tests of general relativity and 
with cosmology. 

2.2.1 General relativity 

The tests of the constancy of fundamental constants take all their importance in the realm of 
the tests of the equivalence principle |536) . Einstein general relativity is based on two indepen- 
dent hypotheses, which can conveniently be described by decomposing the action of the theory as 

*S* = »S*grav 5'mattcr. 

The equivalence principle has strong implication for the functional form of S'grav. This principles 
include three hypothesis: 

• the universality of free fall, 

• the local position invariance, 

• the local Lorentz invariance. 

In its weak form (that is for all interactions but gravity), it is satisfied by any metric theory of 
gravity and general relativity is conjectured to satisfy it in its strong form (that is for all interactions 
including gravity). We refer to Ref. |536) for a detailed description of these principles. The weak 
equivalence principle can be mathematically implemented by assuming that all matter fields are 
minimally coupled to a single metric tensor g^^. This metric defines the length and times measured 
by laboratory clocks and rods so that it can be called the physical metric. This implies that the 
action for any matter field, tp say, can be written as 

'S'mattcr(V^, Sm'')- 

This so-called metric coupling ensures in particular the validity of the universality of free-fall. 
Since locally, in the neighborhood of the worldlinc, there always exists a change of coordinates so 
that the metric takes a Minkowskian form at lowest order, the gravitational field can be locally 
"effaced" (up to tidal effects). If we identify this neighborhood to a small lab, this means that 
any physical properties that can be measured in this lab must be independent of where and when 
the experiments are carried out. This is indeed the assumption of local position invariance which 
implies that the constants must take the same value independent of the spacetime point where they 
are measured. Testing the constancy of fundamental constants is thus a direct test of this principle 
and thus of the metric coupling. Interestingly, the tests we are discussing in this review allow 
to extend them much further than the Solar scales and even in the early universe, an important 
information to check the validity of relativity in cosmology. 
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As an example, the action of a point-particle reads 



mc^-g^^{-K)vi'v''At, (2) 



with v^^ = dx^/di. The equation of motion that derives from this action is the usual geodesic 
equation 

a'' = u''V^u^' = 0, (3) 

where — dx^'/cdr, r being the proper time; is the covariant derivative associated with the 
metric 17^1/ and a'^ is the 4-acceleration. Any metric theory of gravity will enjoy such a matter 
Lagrangian and the worldline of any test particle shall be a geodesic of the spacetime with metric 
g^jy, as long as there is no other long range force acting on it (see Ref. [191] for a detailed review 
of motion in alternative theories of gravity) . 

Note that in the Newtonian limit 500 = — 1 + 2<i>Ar/c^ where $Ar is the Newtonian potential. It 
follows that, in the slow velocity limit, the geodesic equation reduces to 

V = a= -V$Ar = gTV, (4) 

hence defining the Newtonian acceleration g^. Remind that the proper time of a clock is related 
to the coordinate time by dr = ^/—goa<^t. Thus, if one exchanges electromagnetic signals between 
two identical clock in a stationary situation, the apparent difference between the two clocks rates 
will be 

_ $jv(2)-$jv(l) 

at lowest order. This is the so called universality of gravitational redshift. 

The assumption of a metric coupling is actually well tested in the Solar system: 

• First, it implies that all non-gravitational constants are spacetime independent, which have 
been tested to a very high accuracy in many physical systems and for various fundamental 
constants; this the subject of this review. 

• Second, the isotropy has been tested from the constraint on the possible quadrupolar shift 
of nuclear energy levels [98l I303| 1418) proving that different matter fields couple to a unique 
metric tensor at the 10~^^ level. 

• Third, the universality of free fall can be tested by comparing the accelerations of two test 
bodies in an external gravitational field. The parameter 7712 defined as 

^712 = 2 , (5) 

|ai + a2| 

can be constrained experimentally, e.g. in the laboratory by comparing the acceleration of a 
Beryllium and a Copper mass in the Earth gravitational field 4 to get 

77Bc,Cu = (-1.9 ±2.5) X 10-12. (6) 
Similarly the comparison of Earth-corc-like and Moon-mantle-like bodies gave [34] 

r/Earth,Moon = (0.1 ± 2.7 ± 1.7) X 10'^^ (7) 

The Lunar Laser ranging experiment [540] , which compares the relative acceleration of the 
Earth and Moon in the gravitational field of the Sun, also set the constraints 

77Earth,Moon = ("l-O ± 1.4) X 10'^^ (8) 
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Constraint 


Body 1 


Body 2 


Ref. 


(-1.9 ±2.5) X 10"'^ 


Be 


Cu 


m 


(0.1 ±2.7 ± 1.7) X 10"" 


Earth-like rock 


Moon-hke rock 


m 


(-1.0 ± 1.4) X 10"" 


Earth 


Moon 


[540j 


(0.3 ± 1.8) X 10~" 


Te 


Hi 




(-0.2 ±2.8) X 10"^^ 


Be 


Al 


p9] 


(-1.9 ± 2.5) X 10"^^ 


Be 


Cu 


pg] 


(5.1 ±6.7) X 10"^^ 


Si/Al 


Cu 


149] 



Table 3: Summary of the constraints on the violation of the universality of free fall 



Note that since the core represents only 1/3 of the mass of the Earth, and since the Earth 
mantle has the same composition of the Moon (and thus shall fall in the same way), one 
looses a factor 3 so that this constraint is actually similar as the one obtained in the lab. 
Further constraints are summarized in Table [H The latter constraint also contains some 
contribution from the gravitational binding energy and thus includes the strong equivalence 
principle. When the laboratory result of Ref. [33] is combined with the LLR results of 
Refs. |539) and [362] . one gets a constraints on the strong equivalence principle parameter, 
respectively 

773,, = (3 ±6) X 10-13 and ry^,, = (-4 ± 5) x IQ-i^. 

Large improvements are expected thanks to existence of two dedicated space mission projects: 
Microscope [4881 and STEP [353 . 

Fourth, the Einstein effect (or gravitational redshift) has been measured at the 2 x 10"'' 
level [HSj. 



We can conclude that the hypothesis of metric coupling is extremely well-tested in the Solar system. 

The second building block of general relativity is related to the dynamics of the gravitational 
sector, assumed to be dictated by the Einstein-Hilbert action 

c3 



■^gi-av = t^zf: I V-9*R*<^ X- (9) 



leTrG 

This defines the dynamics of a massless spin-2 field g*^, called the Einstein metric. General 
relativity then assumes that both metrics coincide, g^^i, = g*^, (which is related to the strong 
equivalence principle), but it is possible to design theories in which this indeed not the case (see 
the example of scalar-tensor theories below; § 15.1.11) so that general relativity is one out of a large 
family of metric theories. 

The variation of the total action with respect to the metric yields the Einstein equations 

1 SttG 

Rp,u - -^Rgtiv = —^Tf,^, (10) 

where T^" = {2 / y^~g) 6 S matter /^gfii^ is the matter stress-energy tensor. The coefficient SttG/c^ is 
determined by the weak-field limit of the theory that should reproduce the Newtonian predictions. 

The dynamics of general relativity can be tested in the Solar system by using the parameterized 
post-Newtonian formalism (PPN). Its is a general formalism that introduces 10 phenomenological 
parameters to describe any possible deviation from general relativity at the first post-Newtonian 
order |536l 1537] (see also Ref. [SSI for a review on higher orders) . The formalism assumes that 
gravity is described by a metric and that it does not involve any characteristic scale. In its simplest 
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form, it reduces to the two Eddington parameters entering the metric of the Schwartzschild metric 
in isotropic coordinates 



Indeed, general relativity predicts p^^" = = i. 

These two phenomenological parameters are constrained (1) by the shift of the Mercury per- 
ihelion |452| which implies that 127^^^^ — (3^^^ — 1| < 3 x 10^^, (2) the Lunar laser ranging 
experiments jMD] which implies that {A/S^^^ - -f^™ - 3| = (4.4 ± 4.5) x IQ-* and (3) by the 
deflection of electromagnetic signals which are all controlled by 7^^^^. For instance the very long 
baseline interferometry [455] implies that 17^^^ — 1| = 4 x 10"** while the measurement of the time 
delay variation to the Cassini spacecraft (51] sets 7^^^ — 1 = (2.1 ± 2.3) x 10~^. 

The PPN formalism does not allow to test finite range effects that could be caused e.g. by a 
massive degree of freedom. In that case one expects a Yukawa-type deviation from the Newton 
potential. 



that can be probed by "fifth force" experimental searches. A characterizes the range of the Yukawa 
deviation of strength a. The constraints on (A, a) are summarized in Ref. [257] which typically 
shows that a < 10^^ on scales ranging from the millimetre to the Solar system size. 

General relativity is also tested with pulsars |1241 1190J and in the strong field regime |419| . For 
more details we refer to Refs. [1261 14911 15361 1537] . Needless to say that any extension of general 
relativity has to pass these constraints. However, deviations from general relativity can be larger 
in the past, as we shall see, which makes cosmology an interesting physical system to extend these 
constraints. 

2.2.2 Varying constants and the universality of free fall 

As the previous description shows, the constancy of the fundamental constants and the univer- 
sality are two pillars of the equivalence principle. Dicke |151] realized that they are actually not 
independent and that if the coupling constants are spatially dependent then this will induce a 
violation of the universality of free fall. 

The connection lies in the fact that the mass of any composite body, starting e.g. from nuclei, 
includes the mass of the elementary particles that constitute it (this means that it will depend on the 
Yukawa couplings and on the Higgs sector parameters) but also a contribution, -Bbinding/c^, arising 
from the binding energies of the different interactions (i.e. strong, weak and electromagnetic) but 
also gravitational for massive bodies. Thus the mass of any body is a complicated function of all 
the constants, rn[ai]. 

It follows that the action for a point particle is no more given by Eq. ^ but by 



where aj is a list of constant including a^.^ but also many others and where the index A in itia 
recalls that the dependency in these constant is a priori different for body of different chemical 
composition. The variation of this action gives the equation of motion 






(11) 




(12) 
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It follows that a test body will not enjoy a geodesic motion. In the Newtonian limit 1700 = 
— 1 + 2$jv/c^, and at first order in w/c, the equation of motion of a test particle reduces to 



a = gN + Sa.A, Ja^ = -c^^/a,.* (Vtti + Q;i^^ (13) 

i 

SO that in the slow velocity (and slow variation) limit it reduces to 

SaA -c^ ^ ^A.i^OLi- 

i 

where we have introduce the sensitivity of the mass A with respect to the variation of the constant 

oii by 

91nrn^ 

This simple argument shows that if the constants depend on time then there must exist an anoma- 
lous acceleration that will depend on the chemical composition of the body A. 

This anomalous acceleration is generated by the change in the (electromagnetic, gravita- 
tional,...) binding energies [1511 12451 I382j but also in the Yukawa couplings and in the Higgs 
sector parameters so that the ofi-dependencies are a priori composition-dependent. As a conse- 
quence, any variation of the fundamental constants will entail a violation of the universality of free 
fall: the total mass of the body being space dependent, an anomalous force appears if energy is 
to be conserved. The variation of the constants, deviation from general relativity and violation of 
the weak equivalence principle are in general expected together. 

On the other hand, the composition dependence of 5a.A and thus of iiab can be used to opti- 
mize the choice of materials for the experiments testing the equivalence principle |118J but also to 
distinguish between several models if data from the universality of free fall and atomic clocks are 
combined [143] . 

From a theoretical point of view, the computation of rjAB will requires the determination of the 
coefficients jAi- This can be achieved in two steps by first relating the new degrees of freedom of 
the theory to the variation of the fundamental constants and then relating them to the variation 
of the masses. As we shall see in § [5l the first issue is very model dependent while the second is 
especially difficult, particularly when one wants to understand the effect of the quark mass, since 
it is related to the intricate structure of QCD and its role in low energy nuclear reactions. 

As an example, the mass of a nuclei of charge Z and atomic number A can be expressed as 

m(yl, Z) = Znip + {A - Z)m^ + ZtUc + + i^^Mi 

where E^ and are respectively the strong and electromagnetic contributions to the binding 
energy. The Bethe-Weizacker formula allows to estimate the latter as 

E^,, = 98.25^^^^a,„ MeV. (15) 

If we decompose the proton and neutron masses as [231] TO(p,n) = ""3 + ^(u,d)"^u + ^(d.u)^'T'd + 
5(p.n)Q!EM where U3 is the pure QCD approximation of the nucleon mass (6u, &d and i3(n_p)/u3 
being pure numbers), it reduces to 

m{A, Z) = {Auz + Es) + {Zbu + Nbd)m^ + {Zbd + Nbu)md 

ZBp + NBn + 98.25^^^^ Mev\a^^, (16) 
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with N — A — Z , the neutron number. For an atom, one would have to add the contribution 
of the electrons, Zm^. This form depends on strong, weak and electromagnetic quantities. The 
numerical coefficients S(n,p) are given explicitly by [23 IJ 



= 0.63MeV B„a, 



■EM 



= -0.13 MeV. 



(17) 



Such estimations were used in the first analysis of the relation between variation of the constant 
and the universality of free fall |130l 1165) but the dependency on the quark mass is still not well 
understood and we refer to Refs. |119l I156| I158| 1207] for some attempts to refine this description. 
For macroscopic bodies, the mass has also a negative contribution 



from the gravitational binding energy. As a conclusion, from (jl6p and (|18p. we expect the mass to 
depend on all the coupling constant, m(aj,Mi Q^w Q^s ^ '^g j •■•)■ 

We shall discuss this issue in more details in § 15.3.21 
2.2.3 Relations with cosmology 

Most constraints on the time variation of the fundamental constants will not be local and related to 
physical systems at various epochs of the evolution of the universe. It follows that the comparison 
of different constraints requires a full cosmological model. 

Our current cosmological model is known as the ACDM (see Ref. ^404j for a detailed de- 
scription). It is important to recall that its construction relies on 4 main hypotheses: (HI) a 
theory of gravity; (H2) a description of the matter components contained in the Universe and their 
non-gravitational interactions; (H3) symmetry hypothesis; and (H4) a hypothesis on the global 
structure, i.e. the topology, of the Universe. These hypotheses are not on the same footing since 
HI and H2 refer to the physical theories. These hypotheses are however not sufficient to solve 
the field equations and we must make an assumption on the symmetries (H3) of the solutions de- 
scribing our Universe on large scales while H4 is an assumption on some global properties of these 
cosmological solutions, with same local geometry. But the last two hypothesis are unavoidable 
because the knowledge of the fundamental theories is not sufficient to construct a cosmological 



The ACDM model assumes that gravity is described by general relativity (HI), that the Uni- 
verse contains the fields of the standard model of particle physics plus some dark matter and a 
cosmological constant, the latter two having no physical explanation at the moment. It also deeply 
involves the Copernican principle as a symmetry hypothesis (H3), without which the Einstein 
equations usually cannot been solved, and assumes most often that the spatial sections are simply 
connected (H4) . H2 and H3 imply that the description of the standard matter reduces to a mixture 
of a pressureless and a radiation perfect fluids. This model is compatible with all astronomical data 
which roughly indicates that fiAo — 0.73, rimato — 0.27, and $7^0 — 0. Cosmology thus roughly 
imposes that |Ao| < Bl, that is < Bq^ - lO^^ m - lO^i GeV"^ 

Hence, the analysis of the cosmological dynamics of the universe and of its large scale structures 
requires the introduction of a new constant, the cosmological constant, associated with a recent 
acceleration of the cosmic expansion, that can be introduced by modifying the Einstein-Hilbert 
action to 




(18) 



model gnnj- 
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Parametre 


Symbol 


Value 


Reduced Hubble constant 


h 


0.73(3) 


baryon-to-photon ratio 


Tj = nh/rij 


6.12(19) X 10'^° 


Photon density 




2.471 X 10"^ 


Dark matter density 




0.105(8) 


Cosmological constant 




0.73(3) 


Spatial curvature 




0.011(12) 


Scalar modes amplitude 


Q 


(2.0 ± 0.2) X 10~® 


Scalar spectral index 


ns 


0.958(16) 


Neutrino density 


flji^ 


(0.0005 - 0.023) 


Dark energy equation of state 


w 


-0.97(7) 


Scalar running spectral index 


as 


-0.05 ±0.03 


Tensor-to-scalar ratio 


T/S 


< 0.36 


Tensor spectral index 


riT 


< 0.001 


Tensor running spectral index 


aT 


? 


Baryon density 




0.0223(7) 



Table 4: Main cosmological parameters in the standard K-CDM model. There are 7 main parame- 
ters (because ^ fJi = OJ <o which one can add 6 more to include dark energy, neutrinos and gravity 
waves. Note that often the spatial curvature is set to Q.k = 0. 



Note however that it is disproportionately large compared to the natural scale fixed by the Planck 
length 

Classically, this value is no problem but it was pointed out that at the quantum level, the vacuum 
energy should scale as Af^, where M is some energy scale of high-energy physics. In such a case, 
there is a discrepancy of 60-120 order of magnitude between the cosmological conclusions and the 
theoretical expectation. This is the cosmological constant problem [524] . 

Two solutions are considered. Either one accepts such a constant and such a fine-tuning and 
tries to explain it on anthropic ground. Or, in the same spirit as Dirac, one interprets it as an 
indication that our set of cosmological hypotheses have to be extended, by either abandoning the 
Copernican principle [505] or by modifying the local physical laws (either gravity or the matter 
sector). The way to introduce such new physical degrees of freedom were classified in Ref. [497] . 
In that latter approach, the tests of the constancy of the fundamental constants are central since 
they can reveal the coupling of this new degree of freedom to the standard matter fields. Note 
however that the cosmological data still favor a pure cosmological constant. 

Among all the proposals quintessence involves a scalar field rolling down a runaway potential 
hence acting as a fluid with an effective equation of state in the range —l<w< 1 if the field 
is minimally coupled. It was proposed that the quintessence field is also the dilaton [230] 1428] 
1494] ■ The same scalar field then drives the time variation of the cosmological constant and of the 
gravitational constant and it has the property to also have tracking solutions [494] . One of the 
underlying motivation to replace the cosmological constant by a scalar field comes from superstring 
models in which any dimensionful parameter is expressed in terms of the string mass scale and the 
vacuum expectation value of a scalar field. However, the requirement of slow roll (mandatory to 
have a negative pressure) and the fact that the quintessence field dominates today imply, if the 
minimum of the potential is zero, that it is very light, roughly of order m ^ 10^'^'^ eV [50] , 

Such a light field can lead to observable violations of the universality of free fall if it is non- 
universally coupled to the matter fields. Carroll [8^ considered the effect of the coupling of this 
very light quintessence field to ordinary matter via a coupling to the electromagnetic field as 
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(pF'^" Ff^^. Chiba and Kohri [95] also argued that an ultra-light quintessence field induces a time 
variation of the coupling constant if it is coupled to ordinary matter and studied a coupling of the 
form (pF^'^Ff^^, as e.g. expected from Kaluza-Klein theories (see below). This was generalized to 
quintessence models with a couplings of the form ZicfyFt'" F^^ [TOl [TTT] fWI 1^2] 1^31 IM9l l399l [528] 
and then to models of runaway dilaton [1351 1136J inspired by string theory (see § I5.4.ip . The 
evolution of the scalar field drives both the acceleration of the universe at late time and the variation 
of the constants. As pointed in Refs. [951 1165| 15271 such models are extremely constrained from 
the bound on the universality of free-fall (see ? 16. 3p . 
We thus have two ways of investigation 

• The field driving the time variation of the fundamental constants does not explain the accel- 
eration of the universe (either it does not dominate the matter content today or its equation 
of state is not negative enough). In such a case, the variation of the constants is disconnected 
from the dark energy problem. Cosmology allows to determine the dynamics of this field dur- 
ing the whole history of the universe and thus to compare local constraints and cosmological 
constraints. An example is given by scalar-tensor theories (see § I5.1.ip for which one can 
compare e.g. primordial nucleosynthesis to local constraints [129] . In such a situation, one 
should however take into account the effect of the variation of the constants on the astro- 
physical observations since it can affect local physical processes and bias e.g. the luminosity 
of supernovae and indirectly modify the distance luminosity-redshift relation derived from 
these observations [ST] 1429] . 

• The field driving the time variation of the fundamental constants is also responsible for the 
acceleration of the universe. It follows that the dynamics of the universe, the level of variation 
of the constants and the other deviations from general relativity are connected [345] so that 
the study of the variation of the constants can improve the reconstruction of the equation 
state of the dark energy [20l (HH 1385 ] l399 ] . 

In conclusion, cosmology seems to require a new constant. It also provides a link between 
the microphysics and cosmology, as forseen by Dirac. The tests of fundamental constants can 
discriminate between various explanations of the acceleration of the universe. When a model is 
specified, cosmology also allows to set stringer constraints since it relates observables that cannot 
be compared otherwise. 

3 Experimental and observational constraints 

This section focuses on the experimental and observational constraints on the non-gravitational 
constants, that is assuming remains constant. 

The various physical systems that have been considered can be classified in many ways. We 
can classify them according to their look-back time and more precisely their space-time position 
relative to our actual position. This is a summarized on Fig. [1] Indeed higher redshift systems 
offer the possibility to set constraints on an larger time scale, but this is at the expense of usually 
involving other parameters such as the cosmological parameters. This is in particular the case 
of the cosmic microwave background or of primordial nucleosynthesis. The systems can also be 
classified in terms of the physics they involve. For instance, atomics clocks, quasar absorption 
spectra and the cosmic microwave background require only to use quantum electrodynamics to 
draw the primary constraints while the Oklo phenomenon, meteorites dating and nucleosynthesis 
require nuclear physics. 

For any system, setting constraints goes through several steps. First we have some observable 
quantities from which we can draw constraints on primary constants, which may not be fundamen- 
tal constants (e.g. the BBN parameters, the lifetime of /3-decayers,...). This primary parameters 
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Time-redshift relation 




Q i I I I \ I I I \ I I I \ I I I \ I I I \ I I I \ I I I L 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Log [1+z] 



Figure 1: (Top): Summary of the systems that have been used to probe the constancy of the 
fundamental constants and their position in a space-time diagram in which the cone represents our 
past light cone. The shaded areas represents the comoving space probed by different tests. (Bottom): 
The look-back time-redshift relation for the standard ACDM model. 
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System 



Observable 



Primary constraints 



Other hypothesis 



Atomic clock 
Oklo phenomenon 
Meteorite dating 
Quasar spectra 
Stellar physics 
21 cm 
CMB 
BBN 



isotopic ratio 
isotopic ratio 
atomic spectra 
element abundances 

Ti/TcMB 

AT/T 



9i, '^EM,^^ 

A 

5p.m. "em 
Bd 



geophysical model 

cloud physical properties 
stellar model 
cosmological model 
cosmological model 



light element abundances Qnp. Tn, me, mN, c^em . cosmological model 



Table 5: Summary of the systems considered to set constraints on the variation of the fundamental 
constants. We summarize the observable quantities, the primary constants used to interpret the 
data and the other hypothesis required for this interpretation. 



must then be related to some fundamental constants such as masses and couplings. In a last step, 
the number of constants can be reduced by relating them in some unification schemes. Indeed 
each step requires a specific modelisation and hypothesis and has its own limitations. This is 
summarized on Table [5l 



3.1 Atomic clocks 

3.1.1 Atomic spectra and constants 

The laboratory constraints on the time variation of fundamental constants are obtained by com- 
paring the long-term behavior of several oscillators and rely on frequency measurements. The 
atomic transitions have various dependencies in the fundamental constants. For instance, for the 
hydrogen atom, the gross, fine and hyperfinc-structures are roughly given by 



2p — Is : v on cRo 



Is : cx cR 



'oottEMfp^. 



respectively, where the Rydberg constant set the dimension, gp is the proton gyromagnetic factor 
and /2 = mc/mp. In the non-relativistic approximation, the transitions of all atoms have similar 
dependencies but two effects have to be taken into account. First, the hyperfine-structures involve 
a gyromagnetic factor gi (related to the nuclear magnetic moment by /ii = QiHn, with = 
eh/2mpc) which are different for each nuclei. Second, relativistic corrections (including the Casimir 
contribution) which also depend on each atom (but also on the type of the transition) can be 
included through a multiplicative function F^ci{o:em)- '^^^ ^ strong dependence on the atomic 
number Z , which can be illustrated on the case of alkali atoms, for which 



21-1/2 



1 - 3(^"em)^ 



l + y(^«EM)'- 



The developments of highly accurate atomic clocks using different transitions in different atoms 
offer the possibility to test a variation of various combinations of the fundamental constants. 



It follows that at the lowest level of description, we can interpret all atomic clocks results in 
terms of the g- factors of each atoms, gi, the electron to proton mass ration fi and the fine-structure 
constant a^,^. We shall thus parameterize the hyperfine and fine-structures frequencies as follows. 

The hyperfine frequency in a given electronic state of an alkali-like atom, such as ^'^^Cs, *^Rb, 

199Hg+, is 

J^hfs ^ RocC X Ams X ffj X al^ X fix ^hfs(a) (20) 
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Atom Transition sensitivity Kc 

"Ti Is - 2s OOO 

*^Rb hf 0.34 

"^Cs ^S'i/2(F = 2)-(F = 3) 0.83 

"^Yb+ ^Sy2~^D,/2 0.9 

"^Hg+ ^Si/2-'A5/2 -3.2 

*^Sr - ^Po 0.06 

^^Al+ - ^Po 0.008 



Table 6: Sensitivity of various transitions on a variation of the fine- structure constant. 



where gi = fJ-i/fi^} is the nuclear g factor. A^fs is a numerical factor depending on each particular 
atom and we have set F^^i = Fms{(^)- Similarly, the frequency of an electronic transition is well- 
approximated by 

t'clcc ^ RooC X Aelec X i^olcc(-Z', a), (21) 

where, as above, Adcc is a numerical factor depending on each particular atom and i^cioc is the func- 
tion accounting for relativistic effects, spin-orbit couplings and many-body effects. Even though 
an electronic transition should also include a contribution from the hyperfine interaction, it is 
generally only a small fraction of the transition energy and thus should not carry any significant 
sensitivity to a variation of the fundamental constants. 

The importance of the relativistic corrections was probably first emphasized in Ref. [417' and 
their computation through relativistic A^-body calculations was carried out for many transitions 
in Refs. |169l 11731 11741 1198] . They can be characterized by introducing the sensitivity of the 
relativistic factors to a variation of a^.^ , 

'^""^ (22) 



Table m summarizes the values of some of them, as computed in Refs. |1741 12091 . Indeed a reliable 
knowledge of these coefficients at the 1% to 10% level is required to deduce limits to a possible 
variation of the constants. The interpretation of the spectra in this context relies, from a the- 
oretical point of view, only on quantum electrodynamics (QED), a theory which is well tested 
experimentally |272J so that we can safely obtain constraints on (a^.^^ , /i, g^), still keeping in mind 
that the computation of the sensitivity factors required numerical A^-body simulations. 



From an experimental point of view, various combinations of clocks have been performed. It is 
important to analyze as much species as possible in order to rule-out species-dependent systematic 
effects. Most experiments are based on a frequency comparison to caesium clocks. The hyperfine 
splitting frequency between the F = 3 and F = A levels of its ground state at 9.192 GHz has 

been used for the definition of the second since 1967. One limiting effect, that contributes mostly 
to the systematic uncertainty, is the frequency shift due to cold collisions between the atoms. On 
this particular point, clocks based on the hyperfine frequency of the ground state of the rubidium 
at 6.835 GHz, are more favorable. 



3.1.2 Experimental constraints 

We present the latest results that have been obtained and refer to § III.B.2 of FCV [495, for earlier 
studies. They all rely on the developments of new atomic clocks, with the primarily goal to define 
better frequency standards. 
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Clock 1 


Clock 2 

d In / ■'clock 1 \ 
dt V "'clock, / 


Constraint (yr ^) 


Constants dependence 


Reference 




^^^Cs 


(0.2 ± 7.0) X lO"'" 


9Rb EM 


[344| 




^^^Cs 


(-0.5 ±5.3) X 10"^^ 


m\ 


iH 


i^^Cs 


(-32 ± 63) X 10"^^ 




[197] 


199jjg+ 


i^^Cs 


(0.2 ±7) X lO"^'^ 




[Ml 


199Hg+ 


i^^Cs 


(3.7 ±3.9) X 10-^"^ 




Em 


171 Yb+ 


i^^Cs 


(-1.2 ±4.4) X lO"^'^ 


- 1 93 

SCsMOem 


403] 


171 Yb+ 


"3CS 


(-0.78 ± 1.40) X lO"^'^ 




.402] 




133CS 


(-1.0 ± 1.8) X lO"^'^ 


- 2.77 
SCsMOem 


m 






(-2.7 ± 2.6) X 10"^^ 




[99] 


27^1+ 


199Hg+ 


(-5.3 ± 7.9) X lO"^'' 


-3.208 

EM 


[434] 



Table 7: Summary of the constraints obtained from the comparisons of atomic clocks. For each 
constraint on the relative drift of the frequency of the two clocks, we provide the dependence in the 
various constants, using the numbers of Table\^ 



• Rubidium: The comparison of the hyperfine frequencies of the rubidium and caesium in their 
electronic ground state between 1998 and 2003, with an accuracy of order 10~^^, leads to the 
constraint j344] 

^ Inf^) =(0.2±7.0)xlO-iVr-^ (23) 



dt \ vc 

With one more year of experiment, the constraint dropped to [57] 



^. In f ^ ) = (-0.5 ± 5.3) x 10-16 yr"!. (24) 



dt \ vc, 

From Eq. (I^U)) . and using the values of the sensitivities Kq, we deduce that this comparison 
constrains 

• Atomic hydrogen: The Is — 2s transition in atomic hydrogen was compared tp the ground 
state hyperfine splitting of caesium |197j in 1999 and 2003, setting an upper limit on the 
variation of z^h of (— 29± 57) Hz within 44 months. This can be translated in a relative drift 

^In f— ^ = (-32±63) X 10-i6yr-i. (25) 
dt V '^Cs J 

Since the relativistic correction for the atomic hydrogen transition nearly vanishes, we have 
i^H ^ Roo so that 

'^Cs - 2 83 

— cx gcs^^a^Z . 

• Mercury: The i^^Hg+ — ^£'5/2 optical transition has a high sensitivity to a^-^ (see 
Table [5]) so that it is well suited to test its variation. The frequency of the i^^Hg+ electric 
quadrupole transition at 282 nm was thus compared to the ground state hyperfine transition 
of caesium during a two year period, which lead to [56] 

Ain = (0.2±7) X 10-i5yr-i. (26) 

d^ V '^Cs J 
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2002 2003 2004 2005 2006 2007 2008 2009 

Year 

Figure 2: Evolution of the comparison of different atomic clocks summarized in Table^ 

This was improved by a comparison over a 6 year period [215] to get 

A In f = (3.7 ± 3.9) X IQ-^^ yr-^. (27) 

at \ vcs J 

While i^cs is stih given by Eq. (pOj) . i/Hg is given by Eq. (PT|) . Using the sensitivities of Tabled 
we conclude that this comparison test the stability of 

^^Cs _ 6.05 

— oc gcsMa™ • 

• Ytterbium: The ^6*1/2 — electric quadrupole transition at 688 THz of -^^^Yb+ was 

compared to the ground state hyperfine transition of caesium. The constraint of Ref. |403j 
was updated, after comparison over a six year period, which lead to j402J 

^ In ^ (-0.78 ± 1.40) X IQ-^^ yr^i. (28) 



At \ VC; 

Proceeding as previously, this tests the stability of 

^^Cs _ 1.93 

— ocffc.M«,„ ■ 

Strontium: The comparison of the ^Sq — "^Pq transition in neutral ^''Sr with a caesium clock 
was performed in three independent laboratories. The combination of these three experi- 
ments [59^ leads to the constraint 

A In f = (-1.0 ± 1.8) X 10-15 yr-i. (29) 
at \ vcs J 

Proceeding as previously, this tests the stability of 

^^Cs - 2 77 
oc qn^ua 
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• Atomic dyprosium: It was suggested in Refs. |1741ll73j (see also Ref. }172j for a computation 
of the transition amplitudes of the low states of dyprosium) that the electric dipole (El) 
transition between two nearly degenerate opposite-parity states in atomic dyprosium should 
be highly sensitive to the variation of aj^^ . It was then demonstrated |380j that a constraint of 
the order of 10^^*/yr can be reached. The frequencies of nearly of two isotopes of dyprosium 
were monitored over a 8 months period [99^ showing that the frequency variation of the 
3.1-MHz transition in ^'^^Dy and the 235-MIIz transition in ^''^^y ^^.g 9.o±6.7 Hz/yr and 
-0.6±6.5 Hz/yr, respectively. These provides the constraint 

^ = (-2.7 ±2.6) X 10-15 yr-\ (30) 

at la level, without any assumptions on the constancy of other fundamental constants. 

• Aluminium and mercury single-ion optical clocks: The comparison of the ^Sq — ^Pq transition 
in ^^Al+ and ~ ^^^12 in i^^Hg+ over a year allowed to set the constraint |434] 

^ In f = (^5.3 ± 7.9) X 10^" yr-\ (31) 

At \iAigJ 

Proceeding as previously, this tests the stability of 

^ "em ' 

which directly set the constraint 

^ = (-1.6±2.3) X 10-"yr-\ (32) 
since it depends only on a^-^ . 



While the constraint (|32p was obtained directly from the clock comparison, the other studies need 
to be combined to disentangle the contributions of the various constants. As an example, we first 
use the bound (j32p on a^.^ , we can then extract the two following bounds 



iL In ( ^ ) = (0.48 ± 6.68) x 10"^^ yr-\ In (gcsA) = (4-67 ± 5.29) x 10"^^ yr-\ (33) 

di V5Rb/ dt 

on a time scale of a year. Wc cannot lift the degeneracies further with this clock comparison since 
that would require a constraint on the time variation of All these constraints are summarized 
in Table El and Fig.[21 



A solution is to consider diatomic molecules since, as first pointed out by Thomson [482] . 
molecular lines can provide a test of the variation of fi. The energy difference between two adja- 
cent rotational levels in a diatomic molecule is inversely proportional to A/r^^, r being the bond 
length and M the reduced mass, and the vibrational transition of the same molecule has, in first 
approximation, a -y/M dependence. For molecular hydrogen M = mp/2 so that the compari- 
son of an observed vibro-rotational spectrum with a laboratory spectrum gives an information on 
the variation of and m-a_. Comparing pure rotational transitions with electronic transitions 
gives a measurement of pL. It follows that the frequency of vibro-rotation transitions is, in the 
Born-Oppenheimer approximation, of the form 

v^Ei (c^j^^ + c^^^ + c„, y) (34) 
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where 0^,,^,^ , c^;^ and c^^^ are some numerical coefficients. 



The comparison of the vibro-rotational transition in the molecule SF6 was compared to a 
caesium clock over a two-year period, leading to the constraint |460j 

'i. In ^ (1.9 ± 0.12 ± 2.7) x IQ-^^ yr^^, (35) 



where the second error takes into account uncontrolled systematics. Now, using again Table HI we 
deduce that 

i^SF6 -1/2 -2.83/ 

It can be combined with the constraint psp which enjoys the same dependence to caesium to infer 
that 

1^ = (-3.8±5.6) X lO^^Vr"^- (36) 

Combined with Eq. (jSSp . we can thus obtain independent constrains on the time variation of (;cs, 
ffRb and 11. 

3.1.3 Physical interpretation 

The theoretical description must be pushed further if ones wants to extract constraints on constant 
more fundamental than the nuclear magnetic moments. This requires to use quantum chromody- 
namics. In particular, it was argued than within this theoretical framework, one can relate the 
nucleon factors in terms of the quark mass and the QCD scale [198] . Under the assumption of a 
unification of the three non- gravitational interaction (see ? 16. Sp . the dependence of the magnetic 
moments on the quark masses was investigated in Ref. |209| . The magnetic moments, or equiva- 
lently the ^-factors, are first related to the ones of the proton and a neutron to derive a relation of 
the form 

Refs. |198II209) argued that these g-factors mainly depend on the light quark mass mq = ^ {rrin+m^) 
and ms, respectively for the up, down and strange quarks, that is in terms of Xq = mq/Aqco and 
Xs = mg/AqcD- Using a chiral perturbation theory, it was deduced that 

n ii--0.087x--0.013 „ _ y-0.118 1^-0.0013 

SO that for a hyperfine transition 

Both coefficients can be computed, leading to the possibility to draw constraints on the independent 
time variation of Xq, X^ and X^- 

To simplify, we may assume that Xq cx Xs, which is motivated by the Higgs mechanism of mass 
generation, so that the dependence in the quark masses reduces to k = ^(^q + Ks)- For instance, 
we have 

KCs = 0.009, KRb = -0.016, KH = -0.10. 

For hyperfine transition, one further needs to take into account the dependence in ji that can be 
described [213] by 

mp ^ 3AQCDXq" °3^Xf 



q 

so that the hyperfine frequencies behaves as 



VMS CX Xe 
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in the approximation Xq oc Xg and where Xc = mc/AqcD- This allows to get independent 
constraints on the independent time variation of X^, Xq and Q!j,m- Indeed, these constraints are 
model-dependent and, as an example. Table III of Ref. [209 compares the values of the sensitivity 
K when different nuclear effects are considered. For instance, it can vary from 0.127, 0.044 to 0.009 
for the caesium according to whether one includes only valence nucleon, non-valence non-nucleon 
or effect of the quark mass on the spin-spin interaction. It is thus a very promising framework 
which still needs to be developed and the accuracy of which must be quantified in details. 

3.1.4 Future evolutions 

Further progresses in a near future are expected mainly through three types of developments: 

• New systems: Many new systems with enhanced sensitivity |170l 12011 12031 12041 1416) to some 
fundamental constants have recently been proposed. Other atomic systems are considered, 
such as e.g. the hyperfine transitions in the electronic ground state of cold, trapped, hydrogen- 
like highly charged ions [IB] 12001 1444) . or ultra-cold atom and molecule systems near the 
Feshbach resonances [97], where the scattering length is extremely sensitive to 

Concerning diatomic molecules, it was shown that this sensitivity can be enhanced in tran- 
sitions between narrow close levels of different nature [13]. In such transitions, the fine 
structure mainly depends on the fine-structure constant, i^is (^Q^em)^-^ooC, while the vi- 
brational levels depend mainly on the electron-to-proton mass ratio and the reduced mass 
of the molecule, ^ Mr ^^'^ ft^^^RooC. There could be a cancellation between the two fre- 
quencies when V = vm — ^ with n a positive integer. It follows that Sv/v will be 
proportional to K = vm/i^ so that the sensitivity to a-^-^ and /i can be enhanced for these 
particular transitions. A similar effect between with hyperfine-structures, for which the sen- 
sitivity to Q!jj,^j can reach 600 for instance for ^■^^La'^^S or silicon monobrid [H] that allows 
to constrain a^.^/!"^/^. 

Nuclear transitions, such as an optical clock based on a very narrow ultraviolet nuclear tran- 
sition between the ground and first excited states in the ^^^Th, are also under consideration. 
Using a Walecka model for the nuclear potential, it was concluded [200' that the sensitivity 
of the transition to the fine-structure constant and quark mass was typically 



which roughly provides a 5 order of magnitude amplification, which can lead to a constraint at 
the level of 10^^^ yr^^ on the time variation of Xq. Such a method is promising and would 
offer different sensitivities to systematic effects compared to atomic clocks. However, this 
sensibility is not clearly established since different nuclear calculations do not agree [43j 1246) . 

• Atomic clocks in space (ACES): An improvement of at least an order of magnitude on cur- 
rent constraints can be achieved in space with the PHARAO/ACES project |427[ 1438) of 
the European Spatial Agency. PHARAO (Projet d'Horloge Atomique par Refroidissement 
d'Atomes en Orbite) combines laser cooling techniques and a microgravity environment in 
a satellite orbit. It aims at achieving time and frequency transfer with stability better than 



The SAGAS (Search for anomalous gravitation using atomic sensor) project aims at flying 
highly sensitive optical atomic clocks and cold atom accelerometres on a Solar system tra- 
jectory on a time scale of 10 years. It could test the constancy of the fine-structure constant 
along the sattelite wordline which, in particular, can set a constraint on its spatial variation 




10-16. 



of the order of lO'^ [42711543] . 
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• Theoretical developments: We remind one more time that the interpretation of the experi- 
ments requires a good theoretical understanding of the systems but also that the constraints 
we draw on the fundamental constants such as the quark masses are conditional to our the- 
oretical modelling, hence on hypothesis on a unification scheme as well as nuclear physics. 
The accuracy and the robustness of these steps need to be determined, e.g. by taking the 
dependence in the nuclear radius [153J . 

3.2 The Oklo phenomenom 
3.2.1 A natural nuclear reactor 

Oklo is the name of a town in the Gabon republic (West Africa) where an open-pit uranium mine 
is situated. About 1.8 x 10^ yr ago (corresponding to a redshift of ~ 0.14 with the cosmological 
concordance model), in one of the rich vein of uranium ore, a natural nuclear reactor went critical, 
consumed a portion of its fuel and then shut a few million years later (see e.g. Ref. j50?) for 
more details). This phenomenon was discovered by the French Commissariat a I'Energie Atomique 
in 1972 while monitoring for uranium ores |378] . Sixteen natural uranium reactors have been 
identified. Well studied reactors include the zone RZ2 (about 60 bore-holes, 1800 kg of ^35^ 
fissioned during 8.5 x 10^ yr) and zone RZIO (about 13 bore-holes, 650 kg of ^35^ fissioned during 
1.6 X 10^ yr). 

The existence of such a natural reactor was predicted by P. Kuroda [301] who showed that 
under favorable conditions, a spontaneous chain reaction could take place in rich uranium deposits. 
Indeed, two billion years ago, uranium was naturally enriched (due to the difference of decay rate 
between ^35^ g^^j ^■^^U) and 235u represented about 3.68% of the total uranium (compared with 
0.72% today and to the 3-5% enrichment used in most commercial reactors). Besides, in Oklo the 
conditions were favorable: (1) the concentration of neutron absorbers, which prevent the neutrons 
from being available for the chain fission, was low; (2) water played the role of moderator (the 
zones RZ2 and RZIO operated at a depth of several thousand metres, so that the water pressure 
and temperature was close to the pressurized water reactors of 20 Mpa and 300 C) and slowed 
down fast neutrons so that they can interact with other ^35^ g^j^^j ^3-1 ^j^g reactor was large enough 
so that the neutrons did not escape faster than they were produced. It is estimated that the Oklo 
reactor powered 10 to 50 kW. This explanation is backed up by the substantial depletion of 235u 
well as a correlated peculiar distribution of some rare-earth isotopes. These rare-earth isotopes are 
abundantly produced during the fission of uranium and, in particular, the strong neutron absorbers 
like gl^Sm, la^Eu, gt^Gd and gl^Gd are found in very small quantities in the reactor. 

From the isotopic abundances of the yields, one can extract informations about the nuclear 
reactions at the time the reactor was operational and reconstruct the reaction rates at that time. 
One of the key quantity measured is the ratio g2^Sm/g|''Sm of two light isotopes of samarium which 
are not fission products. This ratio of order of 0.9 in normal samarium, is about 0.02 in Oklo ores. 
This low value is interpreted |461| by the depletion of g^^Sm by thermal neutrons produced by the 
fission process and to which it was exposed while the reactor was active. The capture cross-section 
of thermal neutron by g2^Sm 

n + lfSm^ll'Sm + j (37) 

is dominated by a capture resonance of a neutron of energy of about 0.1 eV (£",- = 97.3 meV 
today). The existence of this resonance is a consequence of an almost cancellation between the 
electromagnetic repulsive force and the strong interaction. 

Shlyakhter |461j pointed out that this phenomenon can be used to set a constraint on the time 
variation of fundamental constants. His argument can be summarized as follows. 

• First, the cross-section (J[n,-y) strongly depends on the energy of a resonance at Er = 
97.3 meV. 
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• Geochemical data allow to determine the isotopic composition of various element, such as ura- 
nium, neodynium, gadolinium and samarium. Gadolinium and neodium allow to determine 
the fluence (integrated flux over time) of the neutron while both gadolinium and samarium 
are strong neutron absorbers. 

• From these data, one deduces the value of the averaged value of the cross-section on the 
neutron flux, ang. This value depends on hypothesis on the geometry of the reactor zone. 

• The range of allowed value of (T149 was translated into a constraint on Er- This step involves 
an assumption on the form and temperature of the neutron spectrum. 

• Er was related to some fundamental constant, which involve a model of the nucleus. 
In conclusion, we have different steps, which all involve assumptions: 

• Isotopic compositions and geophysical parameters are measured in a given set of bore-hold 
in each zone. A choice has to be made on the sample to use, in order e.g. to ensure that 
they are not contaminated. 

• With hypothesis on the geometry of the reactor, on the spectrum and temperature of the 
neutron flux, one can deduce the effective value of the cross-sections of neutron absorbers 
(such as samarium and gadolinium). This requires to solve a network of nuclear reaction 
describing the fission. 

• One can then infer the value of the resonance energy Er, which again depends on the as- 
sumptions on the neutron spectrum. 

• Er needs to be related to fundamental constant, which involves a model of the nucleus and 
high energy physics hypothesis. 

We shall now detail the assumptions used in the various analysis that have been performed since 
the pioneering work of Ref . |461| . 

3.2.2 Constraining the shift of the resonance energy 

Cross sections The cross-section of the neutron capture ([57]) strongly depends on the energy of 
a resonance at Er = 97.3 meV and is well described by the Breit-Wigner formula 



where go = (2J + l)(2s + 1)^^(2/ + 1)^^ is a statistical factor which depends on the spin of 
the incident neutron s = 1/2, of the target nucleus /, and of the compound nucleus J. For 
the reaction (157)) . we have 170 = 9/16. The total width F = Fn -I- F^ is the sum of the neutron 
partial width Fn = 0.533 meV (at Er = 97.3 meV and it scales as ^/E in the center of mass) 
and of the radiative partial width F-y — 60.5 meV. gt^Gd has a resonance at Er — 26.8 meV with 
Fn = 0.104 meV, F^ = 108 meV and g ^ 5/8 while gt^Gd has a resonance a.t Er — 31 A mcV with 
Fn = 0.470 meV, F^ = 106 meV and g = 5/8. 

As explained in the previous section, this cross-section cannot be measured from the Oklo 
data, which allow only to measure its value averaged on the neutron flux n(v,T), T being the 
temperature of the moderator. It is conventionally defined as 



f^(n,7)(-^) 



goTT FnF. 



(38) 



2 rUnE {E- Er^ +T^/4 




(39) 
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where the velocity vq = 2200 m • s ^ corresponds to an energy Eq — 25.3 meV and v — y^2E/mn, 
instead of 

_ _ / cr{n.'()niv,T)vdv 
J n{v, T)vdv 

When the cross-section behaves as cr = cfqVq/v, which is the case for nuclei known as 
absorbers", a — (Tq and does not depend on the temperature, whatever the distribution n{v). In a 
similar way, the effective neutron flux defined 



(/) = Wo / n{v,T)dv (40) 

which differs from the true flux 



'P ~ J "(^' T)vdv. 

However, since d-(f> = acj), the reaction rates are note affected by these definitions. 

Extracting the effective cross-section from the data To "measure" the value of a from the 
Oklo data, we need to solve the nuclear reaction network that controls the isotopic composition 
during the fission. 

The samples of the Oklo reactors were exposed 1378) to an integrated effective fluence J (pdt of 
about 10^^ neutron-cm^^ = I kb^^. It implies that any process with a cross-section smaller than 
1 kb can safely be neglected in the computation of the abundances. This includes neutron capture 
by li'^Sm and gf^Sm, as weU as by If'Gd and gl^Gd. On the other hand, the fission of i2^U, the 
capture of neutron by go'^Nd and by 62^Sm with respective cross-sections (T5 ~ 0.6 kb, 0-143 ^ 0.3 kb 
and (7149 > 70 kb are the dominant processes. It follows that the equations of evolution for the 
number densities iVi47, iVi48, ^^149 and iV235 of 62''Sm, |2*Sm, g|^Sm and takes the form 



dTV 



147 



(/)d< 
diVl48 

4>dt 

dNu9 

(fidt 

diV235 

4)dt 



= -(Ti47iV^47 -I- (T/235 2/147^235 (41) 

= <7i47A^147 (42) 

= -O-ugNug + 0'/235yi49A'^235 (43) 

= -'T5iV235, (44) 



where yi denotes the yield of the corresponding element in the fission of 12'^U and is the fission 
cross-section. This system can be integrated under the assumption that the cross-sections and the 
neutron fiux are constant and the result compared with the natural abundances of the samarium 
to extract the value of (T149 at the time of the reaction. Here, the system has been closed by 
introducing a modified absorption cross-section [122] to take into account both the fission, 
capture but also the formation from the a-decay of gl^Pu. One can instead extend the system by 
considering gl^Pu, and 92^!! (see Ref. [235] ). While most studies focus on the samarium, Ref. |221j 
also includes the gadolinium even though it is not clear whether it can reliably be measured [122] . 
They give similar results. 

By comparing the solution of this system with the measured isotopic composition, one can 
deduce the effective cross-section. At this step, the different analysis |461[ I410| [T22l I22H [3021 14111 
I235j differ from the choice of the data. The measured values of 1T149 can be found in these articles. 
They are given for a given zone (RZ2, RZIO mainly) with a number that correspond to the number 
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Ore 


neutron spectrum 


Temperature (°C) 


(T149 (kb) 


Ai;^ (meV) 


Ref. 


? 


Maxwell 


20 


55 ±8 


0± 20 


[461] 


RZ2 (15) 


Maxwell 


180-700 


75± 18 


-1.5 ± 10.5 


[122] 


RZIO 


Maxwell 


200-400 


91 ±6 


4± 16 


[221] 


RZIO 








-97 ±8 


[221] 




Maxwell + epithermal 


327 


91 ±6 




[302] 


RZ2 


Maxwell + epithermal 




73.2 ± 9.4 


-5.5 ±67.5 


|4TT1 


RZ2 


Maxwell + epithermal 


200-300 


71.5 ±10.0 




[255] 


RZIO 


Majcwell + epithermal 


200-300 


85.0 ±6.8 




[235] 


RZ2+RZ10 








7.2 ± 18.8 


[235] 


RZ2+RZ10 








90.75 ± 11.15 


[535] 



Table 8: Summary of the analysis of the Oklo data. The principal assumptions to infer the value 
of the resonance energy Er are the form of the neutron spectrum and its temperature. 



of the bore-hole and the depth (e.g. in Table 2 of Ref. [122J . SC39-1383 means that we are dealing 
with the bore- hole number 39 at a depth of 13.83 m). Recently, another approach |411[ 1235] was 
proposed in order to take into account of the geometry and details of the reactor. It relies on a 
full-scale Monte-Carlo simulation and a computer model of the reactor zone RZ2 |411| and both 
RZ2 and RZIO f235' and allows to take into account the spatial distribution of the neutron flux. 

Determination of Er To convert the constraint on the effective cross-section, one needs to 
specify the neutron spectrum. In the earlier studies [46 ![ 1410) , a Maxwell distribution, 

"th(«,r) = (^)'^%-^, 

was assumed for the neutron with a temperature of 20° C, which is probably too small. Then uq 
is the mean velocity at a temperature Tq = m,iWQ/2fcB = 20.4° C. Refs. |122l 1221] also assume a 
Maxwell distribution but let the moderator temperature vary so that they deduce an effective cross- 
section <7{Rr,T). They respectively restricted the temperature range to 180° C< T < 700° C and 
200° C< T < 400° C, based on geochemical analysis. The advantage of the Maxwell distribution 
assumption is that it avoids to rely on a particular model of the Oklo reactor since the spectrum 
is determined solely by the temperature. 

It was then noted [302[ 1411] that above an energy of several eV, the neutrons spectrum shifted 
to a l/E tail because of the absorption of neutrons in uranium resonances. The distribution was 
thus adjusted to include an epithermal distribution 

n{v) = (1 - f)nth{v,T) + /ncpi(w), 

with Uepi = vl/v"^ for v > Vc and vanishing otherwise. Vc is a cut-off velocity that also needs to be 
specified. The effective cross-section can then be parameterized |235] as 

a = g{T)<j„ + rol, (45) 

where g{T) is a measure of the departure of a from the 1/v behavior, / is related to the resonance 
integral of the cross-section and ro is the Oklo reactor spectral index. It characterizes the contri- 
bution of the epithermal neutrons to the cross-section. Among the unknown parameters, the most 
uncertain is probably the amount of water present at the time of the reaction. Ref. [235] chooses 
to adjust it so that tq matches the experimental values. 
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These hypothesis on the neutron spectrum and on the temperature, as well as the constraint 
on the shift of the resonance energy, are summarised in Table [51 Many analysis |221[|41lTl235j find 
two branches for AEr = Er — Ero, with one (the left branch) indicating a variation of Er- Note 
that these two branches disappear when the temperature is higher since a{Er,T) is more peaked 
when T decreases but remain in any analysis at low temperature. This shows the importance of 
a good determination of the temperature. Note that the analysis of Ref. 14111 indicates that the 
curves cr{T^ Er) lie appreciably lower than for a Maxwell distribution and that Ref. [221] argues 
that the left branch is hardly compatible with the gadolinium data. 



3.2.3 Prom the resonance energy to fundamental constants 

The energy of the resonance depends a priori on many constants since the existence of such reso- 
nance is mainly the consequence of an almost cancellation between the electromagnetic repulsive 
force and the strong interaction. But, since no full analytical understanding of the energy levels of 
heavy nuclei is available, the role of each constant is difficult to disentangle. 

In his first analysis, Shlyakhter |461) stated that for the neutron, the nucleus appears as a 
potential well with a depth Vq ~ 50 MeV. He attributed the change of the resonance energy to a 
modification of the strong interaction coupling constant and concluded that Ag^/g^ ~ AEr/Vo. 
Then, arguing that the Coulomb force increases the average inter-nuclear distance by about 2.5% 
for A ~ 150, he concluded that Aa^.^, /a^.^, ~ 20Ag^/g^, leading to laEM/Q^EMl < 10~^''yr~^, 
which can be translated to 

\Aa^^/a^J< 1.8x10-^. (46) 



The following analysis focused on the fine-structure constant and ignored the strong interaction. 
Damour and Dyson [122] related the variation of Er to the fine-structure constant by taking into 
account that the radiative capture of the neutron by g|^Sm corresponds to the existence of an 
excited quantum state of gl^Sm (so that Er — E^^q — Eug — rrin) and by assuming that the nuclear 
energy is independent of a^.^ . It follows that the variation of 0^,,^^ can be related to the difference of 
the Coulomb binding energy of these two states. The computation of this latter quantity is difficult 
and requires to be related to the mean-square radii of the protons in the isotopes of samarium. In 
particular this analysis |122| showed that the Bethe-Weizacker formula overestimates by about a 
factor the 2 the a^.^ -sensitivity to the resonance energy. It follows from this analysis that 

AEr , . 

- -1.1 MeV, (47) 



which, once combined with the constraint on AEr-, implies 

- 0.9 X 10-^ < Aa^Ja^^, < 1.2 x 10^^ (48) 

at 2a level, corresponding to the range —6.7 x 10~^''yr~^ < a^^/a^^^ < 5.0 x 10~-^''yr~-^ if a^^ is 
assumed constant. This tight constraint arises from the large amplification between the resonance 
energy (~ 0.1 eV) and the sensitivity 1 McV). The re- analysis of these data and also including 
the data of Ref. [221] with gadolinium, found the favored result a^^^ /a^j^j = (— 0.2±0.8)xl0^^''yr^^ 
which corresponds to 

Aa™/«EM = (-0.36 ± 1.44) x IQ-* (49) 

and the other branch (indicating a variation; see Table [S|) leads to Oem/'^em — (4-9 ± 0.4) x 
10~^^yr^^. This non-zero result cannot be eliminated. 

The more recent analysis, based on a modification of the neutron spectrum lead respectively 

to gn] 

Aa^M /«EM = (3-85 ± 5.65) x 10^^ (50) 
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and p35] 



Aa,,K,, = (-0.65 ±1.75) X 10- 



(51) 



at a 95% confidence level, both using the formahsm of Ref. |122j . 

Olive et al. [395) . inspired by grand unification model, reconsider the analysis of Ref. [122' by 
letting all gauge and Yukawa couplings vary. Working within the Fermi gas model, the over-riding 
scale dependence of the terms which determine the binding energy of the heavy nuclei was derived. 
Parameterizing the mass of the hadrons as rrii cx Aqcd(1 + '^i'Tiq/AqcD + ■ ■ Oj they deduce that 
the nuclear Hamiltonian was proportional to TOq/AqcD at lowest order, which allows to estimate 
that the energy of the resonance is related to the quark mass by 



(2.5-10) X lO^^Aln 



A 



QCD 



(52) 



Using the constraint (H7|) . they first deduced that 



A In 



A, 



QCD 



< (1-4) X 10" 



cx TOq° on the basis of grand unification (see § 16.31 for details), they 



Then, assuming that a 
concluded that 

|Aa,M/a,„|<(2-8)xlO-i°. (53) 

Similarly, Refs. |2061l463ll210j related the variation of the resonance energy to the quark mass. 
Their first estimate [206] assumes that it is related to the pion mass, m^r, and that the main 
variation arises from the variation of the radius R ~ 5fm +1/ m^r of the nuclear potential well of 
depth Vq, so that 



3 X 10«- 



assuming that R ~ 1.2^^/'^ro, tq being the inter-nucleon distance. 

Then, in Ref. [463] , the nuclear potential was described by a Walecka model which keeps only 
the a (scalar) and lu (vector) exchanges in the effective nuclear force. Their masses was related to 
the mass TOs of the strange quark to get mo- c>c to|?'^'* and oc m^'^^. It follows that the variation 
of the potential well can be related to 
The constraint (|47|) then implies that 



of the potential well can be related to the variation of mg. and to^j and thus on mq by y oc 



A In 



Aqcd/ 



< 1.2 X 10 
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By extrapolating from light nuclei where the iV-body calculations can be performed more accu- 
rately, it was concluded |207j that the resonance energy scales as /S.Er — 10(A In ATq- O.lA InOj^^J, 
so that the the constraints from Ref. [411] would imply that lS.\n{X^/ a^^^ < 7 x 10^^. 



In conclusion, this last results illustrate that a detailed theoretical analysis and quantitative 
estimates of the nuclear physics (and QCD) aspects of the resonance shift still remain to be carried 
out. In particular, the interface between the perturbative QCD description and the description in 
term of hadron is not fully understand: we do not know the exact dependence of hadronic masses 
and coupling constant on Aqcd and quark masses. The second problem concerns modelling nuclear 
forces in terms of the hadronic parameters. 

At present, the Oklo data, while being stringent and consistent with no variation, have to be 
considered carefully. While a better understanding of nuclear physics is necessary to understand 
the full constant-dependence, the data themselves require more insight, particularly to understand 
the existence of the left-branch. 
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3.3 Meteorite dating 

Long-lived a- or /3-decay isotopes may be sensitive probes of the variation of fundamental constants 
on geological times ranging typically to the age of the Solar system, t ^ (4—5) Gyr, corresponding 
to a mean redshift of z ~ 0.43. Interestingly, it can be compared with the shallow universe 
quasar constraints. This method was initially pointed out by Wilkinson |535) and then revived 
by Dyson |167) . The main idea is to extract the ajj,„ -dependence of the decay rate and to use 
geological samples to bound its time variation. 

The sensitivity of the decay rate of a nucleus to a change of the fine-structure constant is 
defined, in a similar way as for atomic clocks [Eq. ()22|) ]. as 

^^"^ ■• (54) 



9 In aj,M ' 

A is a function of the decay energy Q. When Q is small, mainly due to an accidental cancellation 
between different contributions to the nuclear binding energy, the sensitivity Sa maybe strongly 
enhanced. A small variation of the fundamental constants can either stabilize or destabilize certain 
isotopes so that one can extract bounds on the time variation of their lifetime by comparing 
laboratory data to geophysical ans Solar system probes. 

Assume some meteorites containing an isotope X that decays into Y are formed at a time i*. 
It follows that 

Nx{t) = Nx*e-^^'-*'\ Nvit) = Nx. [l - e-^^*"**)] + Ny, (55) 
if one assumes the decay rate constant. If it is varying then these relations have to be replaced by 

7V^(t) = iVx*e-^'-^(*'''^*' 

so that the value of Nx today can be interpreted with Eq. (|55p but with an effective decay rate 

A=-^ f\{t')dt'. (56) 

From a sample of meteorites, we can measure {A^x(io), A'^y (to)} for each meteorite. These two 
quantities are related by 

Nvih) = - l] Nxih) + Ny., 

so that the data should lie on a line (since Nx* is a priori different for each meteorite), called an 
"isochron", the slope of which determines A(io — i*)- It follows that meteorites data only provides 
an average measure of the decay rate, which complicates the interpretation of the constraints (see 
Refs. |220[ 1219] for explicit examples). To derive a bound on the variation of the constant we also 
need a good estimation of tg — , which can be obtained from the same analysis for an isotope 
with a small sensitivity Sq, as well as an accurate laboratory measurement of the decay rate. 

3.3.1 Long lived a-decays 

The a-decay rate. A, of a nucleus of charge Z and atomic number A, 

^+^X^^X + ^He, (57) 

is governed by the penetration of the Coulomb barrier that can be described by the Gamow theory. 
It is well approximated by 

- A(a^M,w)exp (^-47rZa^„^) , (58) 
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Element 


Z 


A 


Lifetime (yr) 


Q (MeV) 




Sm 


62 


147 


1.06 X 10^^ 


2.310 


774 


Gd 


64 


152 


1.08 X 10" 


2.204 


890 


Dy 


66 


154 


3 X 10*^ 


2.947 


575 


Pt 


78 


190 


6.5 X 10" 


3.249 


659 


Th 


90 


232 


1.41 X 10^" 


4.082 


571 


U 


92 


235 


7.04 X 10® 


4.678 


466 


U 


92 


238 


4.47 X 10^ 


4.270 


548 



Table 9: Summary of the main nuclei and their physical properties that have been used in a-decay 
studies. 



where v/c ^ ^/Qj2m^ is the escape velocity of the a particle. A is a function that depends 
slowly on a^^^ and Q. It follows that the sensitivity to the fine-structure constant is 

s^ -A.Z^^M^ ( 1 - i^^) . (59) 
70/2^ V 2dlna,,J 

The decay energy is related to the nuclear binding energies B{A, Z) of the different nuclei by 

Q = B{A, Z) + Bc,- B{A + 4, Z + 2) 

with Ba = B{4:, 2). Physically, an increase of aj,„ induces an increase in the height of the Coulomb 
barrier at the nuclear surface while the depth of the nuclear potential well below the top remains 
the same. It follows that a-particle escapes with a greater energy but at the same energy below 
the top of the barrier. Since the barrier becomes thiner at a given energy below its top, the 
penetrability increases. This computation indeed neglects the effect of a variation of on the 
nucleus that can be estimated to be dilated by about 1% if increases by 1%. 

As a first insight, when focusing on the fine-structure constant, one can estimate Sa by varying 
only the Coulomb term of the binding energy. Its order of magnitude can be estimated from the 
Bethe-Weizacker formula 

E,,, = 98.25^^^y^a,„ MeV. (60) 

Table [5] summarizes the most sensitive isotopes, with the sensitivities derived from a semi-empirical 
analysis for a spherical nucleus [395 . They are in good agreement with the ones derived from 
Eq. dnni) (e.g., for ^^^U, one would obtain Sa = 540 instead of Sa = 548). 

The sensitivities of all the nuclei of Table [S] are similar, so that the best constraint on the time 
variation of the fine-structure constant will be given by the nuclei with the smaller A A/ A. 

Wilkinson |535J considered the most favorable case, that is the decay of II^U for which Sa ~ 548 
(see Table [9]). By comparing the geological dating of the Earth by different methods, he concluded 
that the decay constant A of ^^^U, ^ssy g,^^^ '^^'^Th have not changed by more than a factor 3 or 4 
during the last 3 — 4 x 10^ years from which it follows 

|Aa,„/a,,J <8x 10-3. (61) 

This constraint was revised by Dyson 167 who claimed that the decay rate has not changed by 
more than 20%, during the past 2 x 10^ years, which implies 

|Aa™/aEMl <4x 10-^ (62) 

Uranium has a short lifetime so that it cannot be used to set constraints on a longer time scales. 
It is also used to calibrate the age of the meteorites. It was thus suggested [395j to consider ^'*''Sm. 
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Assuming that AA147/A147 is smaller than the fractional uncertainty of 7.5 x 10 ^ of its half-life 

|Aa™/a,Ml ^ xlO-5. (63) 

As for the Oklo phenomena, the effect of other constants has not been investigated in depth. 
It is clear that at lowest order both Q and mp scales as Aqcd so that one needs to go beyond such 
a simple description to determine the dependence in the quark masses. Taking into account the 
contribution of the quark masses, in the same way as for Eq. ([52]) . it was argued that A cx Xq^"^^"*^", 
which leads to |AlnXq| < 10~^. In a grand unify framework, that could lead to a constraint of 
the order of |Alna^„| < 2 x lO^'^. 



3.3.2 Long lived /3-decays 

Dicke [150] stressed that the comparison of the rubidium-strontium and potassium-argon dating 
methods to uranium and thorium rates constrains the variation of . 

As long as long-lived /3-decay isotopes are concerned for which the decay energy Q is small, we 
can use a non-relativistic approximation for the decay rate 

A = A±Qf± (64) 

respectively for /3~-decay and electron capture. A± are functions that depend smoothly on a^^^ 
and which can thus be considered constant, p+ = i + 2> and p_ = 2£ -I- 2 are the degrees of 
forbiddenness of the transition. For high-Z nuclei with small decay energy Q, the exponent p 
becomes p = 2 + yj\ — a^^Z"^ and is independent of i. It follows that the sensitivity to a variation 
of the fine-structure constant is 

d In aj,M 

The second factor can be estimated exactly as for a-decay. We note that A± depends on the 
Fermi constant and on the mass of the electron as A± (x G^m^QP . This dependence is the same 
for any /3-decay so that it will disappear in the comparison of two dating methods relying on two 
different /3-decay isotopes, in which case only the dependence on the other constants appear again 
through the nuclear binding energy. Note however that comparing a a- to a /3- decay may lead to 
interesting constraints. 



We refer to § III. A. 4 of FVC |495| for earlier constraints derived from rubidium-strontium, 
potassium- argon and we focus on the rhenium-osmium case, 

Re ]fOs + De + e- (66) 

first considered by Peebles and Dicke [401] . They noted that the very small value of its decay 
energy Q = 2.6 keV makes it a very sensitive probe of the variation of a^-^. In that case p ~ 2.8 
so that Sa — —18000; a change of 10~^% of a^.^ will induce a change in the decay energy of order 
of the keV, that is of the order of the decay energy itself. Peebles and Dicke [401] did not have 
reliable laboratory determination of the decay rate to put any constraint. Dyson |166j compared 
the isotopic analysis of molybdenite ores (Ai87 = (1.6 ±0.2) x 10~^^yr~^), the isotopic analysis of 
14 iron meteorites (A187 = (1.4± 0.3) x 10~^^ yr~^) and laboratory measurements of the decay rate 
(Ai87 = (1.1 ± 0.1) X 10~^^ yr^^). Assuming that the variation of the decay energy comes entirely 
from the variation of a^^^, he concluded that \Aa^^/a^^\ < 9 x lO"* during the past 3 x 10^ 
years. Note that the discrepancy between meteorite and lab data could have been interpreted as a 
time- variation of a^.^ , but the laboratory measurement were complicated by many technical issues 
so that Dyson only considered a conservative upper limit. 
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The modelisation and the computation of Sq were improved in Ref. [395], foUowing the same 
hnes as for a-decay. 



AAi87 „AQ _ _ |^20MeV^^ Aa^ 



187 



-2.2 X 10'' 



Q \ Q 

if one considers only the variation of the Coulomb energy in Q. A similar analysis |146) leads to 
AlnAi87 ^ 10'A\n[a-^-^X-^-^{X^ - 

The dramatic improvement in the meteoric analysis of the Re/Os ratio |464j led to a recent 
reanalysis of the constraints on the fundamental constants. The slope of the isochron was deter- 
mined with a precision of 0.5%. However, the Re/Os ratio is inferred from iron meteorites the 
age of which is not determined directly. Models of formation of the Solar system tend to show 
that iron meteorites and angrite meteorites form within the same 5 million years. The age of 
the latter can be estimated from the 207p]-,_208p|-, niethod which gives 4.558 Gyr [3351 so that 



Ai87 = (1.666 ±0.009) X 10 



-11 



We could thus adopt [395] 



AA 



187 



A 



187 



< 5 X 10"^ 



However, the meteoritic ages are determined mainly by 238y dating so that effectively we have 
a constraint on the variation of Ai87/A238- Fortunately, since the sensitivity of ^38^ jg much 
smaller than the one of the rhenium, it is safe to neglect its effect. Using the recent laboratory 
measurement [331] (Ai87 — (—1.639 ± 0.025) x 10^^^ yr~^), the variation of the decay rate is not 
given by the dispersion of the meteoritic measurement, but by comparing to its value today, so 
that 

AAi87 (g^) 



A 



187 



= -0.016 ±0.016. 



The analysis of Re. |396| . following the assumption of Ref. |395J . deduced that 

Aa™/aBM =(-8±16) x 10"^ 



(68) 



at a 95% confidence level. 

As pointed out in Ref. j 220l 1219] . this constraints really represents a bound on the average 
decay rate A since the formation of the meteorites. This implies in particular that the redshift 
at which one should consider this constraint depends on the specific functional dependence A(t). 
It was shown that well-designed time dependence for A can obviate this limit, due to the time 
average. 



3.3.3 Conclusions 

Meteorites data allow to set constraints on the variation of the fundamental constants which are 
comparable to the ones set by the Oklo phenomenon. Similar constraints can also bet set from 
spontaneous fission (see § HI. A. 3 of FVC |495) ) but this process is less well understood and less 
sensitive than the a- and /3- decay processes and. 

From an experimental point of view, the main difficulty concerns the dating of the meteorites 
and the interpretation of the effective decay rate. 

As long as we only consider a^.^ , the sensitivities can be computed mainly by considering the 
contribution of the Coulomb energy to the decay energy, that reduces to its contribution to the 
nuclear energy. However, as for the Oklo phenomenon, the dependencies in the other constants, 
Xq, Gp, /i. . . , require a nuclear model and remain very mo del- dependent. 
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3.4 Quasar absorbtion spectra 
3.4.1 Generalities 

Quasar (QSO) absorption lines provide a powerful probe of the variation of fundamental constants. 
Absorption lines in intervening clouds along the line of sight of the QSO give access to the spectra 
of the atoms present in the cloud, that it is to paleo-spectra. The method was first used by 
Savedoff [441] who constrained the time variation of the fine-structure constraint from the doublet 
separations seen in galaxy emission spectra. For general introduction to these observations, we 
refer to Refs. [4(171 l470l [276] . 

Indeed, one cannot use a single transition compared to its laboratory value since the expansion 
of the universe induces a global redshifting of all spectra. In order to tackle down a variation of 
the fundamental constants, one should resort on various transitions and look for chromatic effects 
that can indeed not be reproduce by the expansion of the universe which acts chromatically on all 
wavelengths. 

To achieve such a test, one needs to understand the dependencies of different types of transitions, 
in a similar way as for atomic clock experiments. Refs. [174| I168J suggested to use the convenient 
formulation 



in order to take into account the dependence of the spectra on the fine-structure constant, u 
is the energy in the rest-frame of the cloud, that is at a redshift z, is the energy measured 
today in the laboratory, q and (72 are two coefficients that determine the frequency dependence 
on a variation of and that arise from the relativistic corrections for the transition under 
consideration. The coefficient q is typically an order of magnitude larger than 52 so that the 
possibility to constrain a variation of the fine-structure constant is mainly determined by qi. 
This coefficients were computed for a large set of transitions, first using a relativistic Hartree- 
Fock method and then using many-body perturbation theory. We refer to Refs. [1741 l44l [M] for 
an extensive discussion of the computational methods and a list of the g-coefficients for various 
transitions relevant for both quasar spectra and atomic clock experiments. Fig. [3] summarizes some 
of these results. The uncertainty in q are typically smaller than 30cm~^ for Mg, Si, Al and Zn, 
but much larger for Cr, Fe and Ni due to their more complicated electronic configurations. The 
accuracy for from dedicated laboratory measurements now reach 0.004 cm~^. It is important 
to stress that the form ([69|) ensures that errors in the g-coefficients cannot lead to a non zero 
detection of Aa^„ . 

The shift between two lines is easier to measure when the difference between the g-coefficients 
of the two lines is large, which occurs e.g. for two levels with large q of opposite sign. Many 
methods were developed to take this into account. The alkali doublet method (AD) focuses on 
the fine-structure doublet of alkali atoms. It was then generalized to the many-multiplet method 
(MM) which uses correlations between various transitions in different atoms. As can be seen on 
Fig. O some transitions are almost insensitive to a variation of a^.^ . This is the case of Mgii, 
which can be used as an anchor, i.e. a reference point. To obtain strong constraints one can 
either compare transitions of light atoms with those of heavy atoms (because the a^-^^ dependence 
of the ground state scales as Z^) or compare s — p and d — p transitions in heavy elements (in 
that case, the relativistic correction will be of opposite signs). This latter effect increases the 
sensitivity and strengthens the method against systematic errors. The results of this method relies 
however on two assumptions: (i) ionization and chemical homogeneity and (ii) isotopic abundance 
of Mgii close to the terrestrial value. Even though these are reasonable assumptions, one cannot 
completely rule out systematic biases that they could induce. The AD method completely avoids 
the assumption of homogeneity because, by construction, the two lines of the doublet must have 
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Figure 3: Summary of the values of some coefficients entering the parameterization i69\) and 
necessary to interpret the QSO absorption spectra data. From Ref. \87Cf^ 



the same profile. Indeed the AD method avoids the implicit assumption of the MM method that 
chemical and ionization inhomogeneities are negligible. Another way to avoid the influence of 
small spectral shift due to ionization inhomogeneities within the absorber and due to possible 
non-zero offset between different exposures was to rely on different transitions of a single ion in 
individual exposure. This method has been called the Single ion differential alpha measurement 
method (SIDAM). 

Most studies are based on optical techniques due to the profusion of strong UV transitions that 
are redshifted into the optical band (this includes AD, MM, SIDAM and it implies that they can 
be applied only above a given redshift, e.g. SilV at z > 1.3, FeiiA1608 at z > 1) or on radio 
techniques since radio transitions arise from many different physical effects (hyperfine splitting 
and in particular Hi 21 cm hyperfine transition, molecular rotation. Lambda-doubling, etc). In 
the latter case, the line frequencies and their comparisons yield constraints on different sets of 
fundamental constants including a-^„ , t/p and /z. These techniques are thus complementary since 
systematic effects are different in optical and radio regimes. Also the radio techniques offer some 
advantages: (1) to reach high spectral resolution (< 1 km/s), alleviating in particular problems 
with line blending and the use of e.g. masers allow to reach a frequency calibration better than 
roughly 10 m/s; (2) in general, the sensitivity of the line position to a variation of a constant is 
higher; (3) the isotopic lines are observed separately, while in optical there is a blend with possible 
differential saturations (see e.g. Ref. [108] for a discussion). 

Let us first emphasize that the shifts in the absorption lines to be detected are extremely small. 
For instance a change of a-^^^ of order 10~^ corresponds a shift of at most 20 mA for a redshift 
of z ~ 2, which would corresponds to a shift of order ~ 0.5 km/s, or to about a third of a pixel 
at a spectral resolution of i? 40000, as achievd with Keck/HIRES or VLT/UVES. As we shall 
discuss later, there are several sources of uncertainty that hamper the measurement. In particular, 
the absorption lines have complex profiles (because they result from the propagation of photons 
through a highly inhomogeneous medium) that are fitted using a combination of Voigt profiles. 
Each of these components depends on several parameters including the redshift, the column density 
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and the width of the hne (Doppler parameter) to which one now needs to add the constants that 
are assumed to be varying. These parameters are constrained assuming that the profiles are the 
same for all transitions, which is indeed a non-trivial assumption for transitions from different 
species (this was one of the driving motivation to use transition from a single species and of the 
SID AM method). More important, the fit is usually not unique. This is not a problem when the 
lines are not saturated but it can increase the error on aj,„ by a factor 2 in the case of strongly 
saturated lines [90] . 

3.4.2 Alkali doublet method (AD) 

The first method used to set constraint on the time variation of the fine-structure constant relies 
on fine-structure doublets splitting for which 

It follows that the relative separation is proportional a^MJ ^vjv cx c?^^ so that the variation of 
the fine structure constant at a redshift z can be obtained as 




where c,. ^ 1 is a number taking into account the relativistic corrections. This expression is indeed 
a simple approach of the alkali doublet since one should, as for atomic clocks, take into account 
the relativistic corrections more precisely. Using the formulation one can deduce that 

_ (5g -I- bq2 
~ Sq + 2Sq2 ' 

where the Sq are the differences between the g-coefficients for the doublet transitions. 

Several authors have applied the AD method to doublets of several species such as e.g. Civ, 
Nv, Ovi, Mgll, Ahll, Sin, Siiv. We refer to § III.3 of FVC |495| for a summary of their results 
(see also Ref. |316j ) and focus on the three most recent analysis, based on the Siiv doublet. In 
this particular case, q = 766 (resp. 362) cm~^ and q2 = 48 (resp. —8) cm~^ for Siiv A1393 (resp. 
A1402) so that c,- = 0.8914. The method is based on a minimization of multiple component 
Voigt profile fits to the absorption features in the QSO spectra. In general such a profile depends 
on three parameters, the column density N, the Doppler width (6) and the redshift. It is now 
extended to include Aa^.^ /agj, . The fit is carried out by simultaneously varying these parameters 
for each component. 

• Murphy et al. [374] analyzed 21 Keck/HIRES Siiv absorption systems toward 8 quasars to 
obtain the weighted mean of the sample, 

^"EM/ttEM = (-0-5± 1.3) X 10"^ 2.33<z< 3.08, (70) 

with a mean redshift of z = 2.6. The S/N ratio of these data is in the range 15-40 per pixel 
and the spectral resolution is i? ^ 34000. 

• Chand et al. [90] analyzed 15 Siiv absorption systems selected from a ESO-UVES sample 
containing 31 systems (eliminating contaminated, saturated or very broad systems; in par- 
ticular a lower limit on the column density was fixed so that both lines of the doublets are 
detected at more than da) to get the weighted mean, 

Aa^^/a^^ = (-0.15 ±0.43) X 10"^ 1.59 < z < 2.92. (71) 
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The improvement of the constraint arises mainly from a better S/N ratio, of order 60-80 
per pixel, and resolution R ^ 45000. Note that combining this result with the previous 
one (|70l in a weighted mean would lead to Aa^^Ja^^^.^ — (— 0.04± 0.56) x 10^^ in the range 
1.59 < z < 3.02 

• The analysis |346| of seven CiV systems and two Siiv systems in the direction of a single 
quasar, obtained by the VLT-VES (during the science verification) has led to 

^"EM/aEM = (-3.09 ± 8.46) x 10"^ 1.19 < z < 1.84. (72) 

This is less constraining than the two previous analysis, mainly because the g-coefhcients are 
smaller for Civ (see Ref. |406] for the calibration of the laboratory spectra) 

One limitation may arise from the isotopic composition of silicium. Silicium has three naturally 
occurring isotopes with terrestrial abundances ^^Si : ^^Si : ^°Si = 92.23 : 4.68 : 3.09 so that each 
absorption line is a composite of absorption lines from the three isotopes. It was shown that this 
effect of isotopic shifts |374| is however negligible in the case of Siiv. 

3.4.3 Many multiplet method (MM) 

A generalization of the AD method, known as the many-mulptiplet was proposed in Ref. [175] . 
It relies on the combination of transitions from different species. In particular, as can be seen 
on Fig. [21 some transitions are fairly unsensitive to a change of the fine-structure constant (e.g. 
Mgii or Mgi, hence providing good anchors) while others such as Fell are more sensitive. The first 
implementation |518j of the method was based on a measurement of the shift of the Fell (the rest 
wavelengths of which are very sensitive to aj,^^) spectrum with respect to the one of Mgii. This 
comparison increases the sensitivity compared with methods using only alkali doublets. Two series 
of analysis were performed during the past ten years and lead to contradictory conclusions. The 
accuracy of the measurements depends on how well the absorption line profiles are modelled. 

Keck/HIRES data. The MM-method was first applied in Ref. [518^ who analyzed one transition 
of the Mgii doublet and five Fell transitions from three multiplets. Using 30 absorption systems 
toward 17 quasars, they obtained 

Aa^^Ja^^, = (-0.17 ±0.39) x 10-^ 0.6 < z < 1 
Aa^^/a^^ = (-1.88 ± 0.53) x 10"^ l<z< 1.6. 

This was the first claim that a constant may have varied during the evolution of the universe. It 
was later confirmed in a re-analysis [372| 1519] of the initial sample and by including new optical 
QSO data to reach 28 absorption systems with redshift z = 0.5 — 1.8 plus 18 damped Lyman- 
a absorption systems towards 13 QSO plus 21 Siiv absorption systems toward 13 QSO . The 
analysis used mainly the multiplets of Niii, Cm and Znii and Mgi, Mgi, Alii, Aim and Fell was 
also included. The most recent analysis |363J relies on 128 absorbtion spectra, later updated [370] 
to include 143 absorption systems. The more robust estimates is the weighted mean 

^"em/"em = (-0-57 ±0.11) X 10-^ 0.2<z<4.2. (73) 

The resolution for most spectra was R ~ 45000 and the S/N per pixel ranges from 4 to 240, with 
most spectral regions with S/N~ 30. The wavelength scale was calibrated by mean of a Thorium- 
argon emission lamp. This calibration is crucial and its quality is discussed in Ref. |369| 1371) for 
the Keck/HIRES (see also Ref. [137]) as well as Ref. [330] for the VLT/UVES measurements. 

The low-z (z < 1.8) and high-z rely on different ions and transitions with very different ckj^,^- 
dependencies. At low-2;, the Mg transitions are used as anchors against which the large positive 
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shifts in the Fell can be measured. At high-z, different transitions are fitted (Fell, Sii, Cm, Mil, 
Znll, Ahi, Ahii). The two sub-samples respond differently to simple systematic errors due to their 
different arrangement of q-coefficients in wavelength space. The analysis for each sample give the 
weighted mean 

Aa^M/ttEM = (-0-54 ± 0.12) x 10-^ 0.2 < z < 1.8 

Aa,,/a,,, = (-0.74 ± 0.17) x 10-^ 1.8 < z < 4.2, (74) 
with respectively 77 and 66 systems. 

Hunting systematics. While performing this kind of observations a number of problems and 
systematic effects have to be taken into account and controlled. (1) Errors in the determination of 
laboratory wavelengths to which the observations are compared. (2) While comparing wavelengths 
from different atoms one has to take into account that they may be located in different regions of the 
cloud with different velocities and hence with different Doppler shifts. (3) One has to ensure that 
there is no transition not blended by transitions of another system. (4) The differential isotopic 
saturation has to be controlled. Usually quasar absorption systems are expected to have lower 
heavy element abundances. The spatial inhomogeneity of these abundances may also play a role. 
(5) Hyperfine splitting can induce a saturation similar to isotopic abundances. (6) The variation 
of the velocity of the Earth during the integration of a quasar spectrum can also induce differential 
Doppler shift. (7) Atmospheric dispersion across the spectral direction of the spectrograph slit 
can stretch the spectrum. It was shown that, on average, thisv can, for low redshift observations, 
mimic a negative Aa^.^^ /a^.^^ , while this is no more the case for high redshift observations (hence 
empahasizing the complementarity of these observations). (8) The presence of a magnetic field will 
shift the energy levels by Zeeman effect. (9) Temperature variations during the observation will 
change the air refractive index in the spectrograph. In particular, fiexures in the instrument are 
dealt with by recording a calibration lamp spectrum before and after the science expossure and 
the signal-to-noise and stability of the lamp is crucial (10) Instrumental effects such as variations 
of the intrinsic instrument profile have to be controlled. 

All these effects have been discussed in details in Refs. [3711 1372] to argue that none of them 
can explain the current detection. This was recently complemented by a study on the calibration 
since adistortion of the wavelength scale could lead to a non-zero value of Aa^.^. The quality of 
the calibration is discussed in Ref. 369 and shown to have a negligible effect on the measurements 
(a similar result has been obtained for the VLT/UVES data |530) ). 

As we pointed out earlier, one assumption of the method concerns the isotopic abundances of 
Mgii that can affect the low-z sample since any changes in the isotopic composition will alter the 
value of effective rest- wavelengths. This isotopic composition is assumed to be close to terrestrial 
24Mg : 25Mg : "^^Mg = 79 : 10 : 11. No direct measurement of rug = C^Mg + 25Mg)/24Mg in 
QSO absorber is currently feasible due to the small separation of the isotopic absorption lines. It 
was however shown [232], on the basis of molecular absorption lines of MgH that rug generally 
decreases with metallicity. It was also argued that ^"^C is a tracer of ^^Mg and was shown to be 
low in the case of HE 0515-4414 [318. . However, contrary to this trend, it was found [548) that 
rjvig can reach high values for some giant stars in the globular cluster NGC 6752 with metallicity 
[De/H]~ —1.6. This led Ashenfelter et al. [T^ to propose a chemical evolution model with strongly 
enhanced population of intermediate (2 — 8Mq) stars which in their asymptotic giant branch phase 
are the dominant factories for heavy Mg at low metallicities typical of QSO absorption systems, 
as a possible explanation of the low-z Keck/HIRES observations without any variation of a^M- 
It would require that ryig reaches 0.62, compared to 0.27 (but then the UVES/VLT constraints 
would be converted to a detection). However, such modified nucleosynthetic history will lead to 
an overproduction of elements such as P, Si, Al, P above current constraints |193j . 

In conclusion, no compelling evidence for a systematic effect has been raised at the moment. 
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VLT/UVES data. The previous results, and their importance for fundamental physics, led 
another team to check this detection using observations from UVES spectrograph operating on 
the VLT. In order to avoid as much systematics as possible, and based on numerical simulations, 
they apply a series of selection criteria [89j on the systems used to constrain the time variation of 
the fine-structure constant: (1) consider only lines with similar ionization potentials (Mgii, Fell, 
Sin and Alii) as they are most likely to originate from similar regions in the cloud; (2) avoid 
absorption lines contaminated by atmospheric lines; (3) consider only systems with hight enough 
column density to ensure that all the mutiplets are detected at more than 5(7; (4) demand than 
at least one of the anchor lines is not saturated to have a robust measurement of the redshift; (5) 
reject strongly saturated systems with large velocity spread; (6) keep only systems for which the 
majority of the components are separated from the neighboring by more than the Doppler shift 
parameter. 

The advantage of this choice is to reject most complex or degenerate systems, which could 
result in uncontrolled systematics effects. The drawback is indeed that the analysis will be based 
on less systems. 

Ref. [89 (l466j analyzed the observations of 23 absorption systems, fulfilling the above criteria, in 
direction of 18 QSO with a S/N ranging between 50 and 80 per pixel and a resolution R > 44000. 
They concluded that 

Aa^^/a^^, = (-0.06 ± 0.06) x 10'^, OA < z < 2.3, 

hence giving a 3cr constraint on a variation of a^^-^ . 

This analysis was challenged by Murphy, Webb and Flambaum f364| 13661 1365] . Using (quoting 
them) the same reduced data, using the same fits to the absorption profiles, they claim to find 
different individual measurements of Aa^^/a^j, and a weighted mean, 

AaEM/^EM = (-0.44 ± 0.16) X 10~^ 0.4 < z < 2.3, 

which differs from the above cited value. The main points that were raised are (1) the fact that some 
of the uncertainties on Aa^^j^^/aj^-^^ are smaller than a minimum uncertainty that they estimated 
and (2) the quality of the statistical analysis (in particular on the basis of the curves). These 
arguments were responded in Ref. ^46 7. The revision [467] of the VLT/UVES constraint rejects 
two more than Aa deviant systems that were claimed to dominate the reanalysis |366| 1365] and 
concludes that 

^"EM/aEM = (0-01 ±0.15) X 10"^ O.4<0<2.3, (75) 

emphasizing that the errors are probably larger. 

On the basis of the articles |3641 13661 1365] and the answer |467) . it is indeed difficult (without 
having played with the data) to engage one of the parties. This exchange has enlightened some 
differences in the statistical analysis. 

To finish, let us mention that Ref. |358] reanalyzed some systems of Refs. [89l 1466] by means 
of the SIDAM method (see below) and disagree with some of them, claiming for a problem of 
calibration. They also claim that the errors quoted in Ref. ^370) are underestimated by a factor 
1.5. 

Regressional MM (RMM). The MM method was adapted to use a linear regression method [421] 
The idea is to measure the redshift Zi deduced from the transition i and plot Zi as a function of the 
sensitivity coefficient. If Aa^-^^^ ^ then there should exist a linear relation with a slope propor- 
tional to Aq!j,„ /q!j,„ . On a single absorption system (VLT/UVES), on the basis of Fell transition, 
they concluded that 

^^^m/^^m - (-0.4 ± 1.9 ± 2.7,yst) X 10-^ z = 1.15, (76) 
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compared to AQ!j.„/aj,„ — (0.1 ± 1.7) x 10^^ that is obtained with the standard MM technique on 
the same data. This is also consistent with the constraint (|78p obtained on the same system with 
the HARPS spectrograph. 

Open controversy. At the moment, we have to face a situation in which two teams have per- 
formed two independent analysis based on data sets obtained by two instruments on two telescopes. 
Their conclusions do not agree, since only one of them is claiming for a detection of a variation 
of the fine-structure constant. This discrepancy between VLT/UVES and Keck/Hires results is 
yet to be resolved. In particular, they use data from a different telescopes observing a different 
(Southern/Northern) hemisphere. 

Note however a recent analysis |237l of the wavelength accuracy of the Keck/HIRES spectro- 
graph. An absolute uncertainty of Az ^ 10^^, corresponding to AA ^ 0.02 A with daily drift 
of Az ~ 5 X 10~^ and multiday drift of Az 2 x lO"*^. While the cause of this drift remains 
unknown, it is argued [237] that this level of systematic uncertainty makes it difficult to use the 
Keck/HIRES to constrain the time variation of a^.^ (at least for a single system or a small sample 
since the distortion pattern pertains to the echelle orders as they are recorded on the CDD, that 
is it is similar from exposure to exposure, the effect on Aajj,„/aj,M for an ensemble of absorbers 
at different redshifts would be random since the transitions fall in different places with respect to 
the pattern of the disortion). This needs to be confirmed and investigated in more details. We 
refer to Ref. [367] for a discussion on the Keck wavelength calibration error and Ref. |530| for the 
VLT/UVES as well as Ref. [85] for a discussion on the ThAr calibration. 

On the one hand, it is sane that one team has reanalyzed the data of the other and challenge 
its analysis. This would indeed lead to an improvement of the robustness of these results. Indeed 
a similar reverse analysis would also be sane. On the other hand both teams have achieved an 
amazing work in order to understand and quantify all sources of systematics. Both developments, 
as well as the new techniques which are appearing, should hopefully set this observational issue. 
Today, it is unfortunately premature to choose one data set compared to the other. 

A recent data [520] set of 60 quasar spectra (yielding 153 absorption systems) for the VLT was 
used and split at z = 1.8 to get 

(Aa,„/a,„)vLT;.<i.8 = (-0.06 ±0.16) x 10-^ 
in agreement with the former study [467] . while at higher redshift 

(Aa,„/a,„)vLT.>i.8 = (+0-61 ± 0.20) x 10"^ 

This higher component exhibits a positive variation of , that is of opposite sign with respect 
to the previous Keck/HIRES detection |370] 

(Aa™/«™)Kcck;.<i.8 = (-0.54±0.12)xl0-^ (A«™/«BM)Kcck; .>i.8 = (-0.74±0.17)xl0-^ 

It was pointed out that the Keck/HIRES and VLT/UVES observations can be made consistent in 
the case the fine structure constant is spatially varying |520] . Indeed, one can note that they do not 
correspond to the same hemisphere and invoque a spatial variation. Ref. |520] concludes that the 
distribution of a^.^ is well represented by a spatial dipole, significant at 4.1ct, in the direction right 
ascension 17.3 ±0.6 hours and declination —61 ±9 deg (see also Ref. [47ll49] ). This emphasizes the 
difficulty to compare different data sets and shows that the constraints can easily be combined as 
long as they are compatible with no variation but one must care about a possible spatial variation 
otherwise. 
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3.4.4 Single ion differential measurement (SIDAM) 



This method |317j is an adaptation of the MM method in order to avoid the influence of small 
spectral shifts due to ionization inhomogeneities within the absorbers as well as to non-zero offsets 
between different exposures. It was mainly used with Fell which provides transitions with positive 
and negative q-coefhcients (see Fig. [3]). Since it relies on a single ion, it is less sensitive to isotopic 
abundances, and in particular not sensitive to the one of Mg. 

The first analysis relies on the QSO HE 0515-4414 that was used in Ref. f421' to get the 
constraint (j76l) . An independent analysis (358j of the same system gave a weighted mean 

A«,M/aEM = (-0.12 ± 1.79) x 10-^ z = 1.15, (77) 

at Icr. The same system was studied independently, using the HARPS spectrograph mounted on 
the 3.6m telescope at La Silla observatory [91]. The HARPS spectrograph has a higher resolution 
that UVES; R ~ 112000. Observations based on Feii with a S/N of about 30-40 per pixel set the 

constraint 

Aa,„/a,„ = (0.5±2.4)xl0-«, z = 1.15. (78) 
The second constraint [32 H 1358] is obtained from an absorption system toward Q 1101-264, 

Aa™/"™ = (5-66 ± 2.67) x 10-^ z = 1.84, (79) 

These constraints do not seem to be compatible with the results of the Keck/HIRES based on the 
MM method. A potential systematic uncertainty which can affect these constraints is the relative 
shift of the wavelength calibration in the blue and the red arms of UVES where the distant Fe 
lines are recorded simultaneously (see e.g. Ref. [359_^ for discussion of systematics of this analysis). 

3.4.5 Hi-21 cm vs UV: x = a^^g^/pL 

The comparison of UV heavy element transitions with the hyperfine Hi transition allows to ex- 
tract [IHD] 

since the hyperfine transition is proportional to a^^gpfi~^Roo while optical transitions are simply 
proportional to i?oo- It follows that constraints on the time variation of x can be obtained from 
high resolution 21cm spectra compared to UV lines, e.g. of Sin, Fell and/or Mgii, as first performed 
in Ref. [544] in z ~ 0.524 absorber. 

Using 9 absorption systems, there was no evidence for any variation of x [489] , 

Aa;/x (-0.63 ±0.99) X 10"^ 0.23 < z < 2.35, (80) 

This constraints was criticized in Ref. [275] on the basis that the systems have multiple components 
and that it is not necessary that the strongest absorption arises in the same component in both 
type of lines. However the error analysis of Ref. [489] tries to estimate the effect of the assumption 
that the strongest absorption arises in the same component. 

Following Ref. [146) . we note that the systems lie in two widely-separated ranges and that the 
two samples have completely different scatter. It can thus be split in two samples of respectively 
5 and 4 systems to get 

Acc/x = (1.02 ± 1.68) X 10"^ 0.23 < z < 0.53, (81) 
Ax/x = (0.58 ± 1.94) X 10"^ 1.7 < z < 2.35. (82) 

In such an approach two main difficulties arise: (1) the radio and optical source must coincide 
(in the optical QSO can be considered pointlike and it must be checked that this is also the case 



47 



for the radio source), (2) the clouds reponsible for the 21cm and UV absorptions must be localized 
in the same place. The systems must thus be selected with care and today the number of such 
systems is small and are activily looked for [405] . 

The recent detection of 21cm and molecular hydrogen absorption lines in the same damped 
Lyman-a system at Zabs = 3.174 towards SDSS J1337+3152 constrains |469] the variation x to 

Ax/x = -(1.7 ± 1.7) X 10"^ z ^ 3.174. (83) 

This system is unique since it allows for 21cm, H2 and UV observation so that in principle one 
can measure a^^, , x and /i independently. However, as the H2 column density was low, only 
Werner band absorption lines are seen so that the range of sensitivity coefficients is too narrow to 
provide a stringent constraint, A/i//i < 4 x 10~^. It was also shown that the H2 and 21cm are 
shifted because of the inhomogeneity of the gas, hence emphasizing this limitation. Ref. [405] also 
mentioned that 4 systems ai z — 1.3 sets Ax/x = (0.0 ± 1.5) x 10~® and that another system 
at z = 3.1 gives As/x = (0.2 ± 0.5) x 10~^. Note also that the comparison [280] with Ci at 
z - 1.4 - 1.6 towards Q0458-020 and Q2337-011, yields Ax/x = (6.8 ± 1.0) x lO'^ over the 
band o redshift << z >< 1.46. It was argued that, using the existing constraints on An/n, 
this measurement is inconsistent with claims of a smaller value of from the many- mult iplet 
method, unless fractional changes in gp are larger than those in and fi. 



^2 

EA4 



3.4.6 Hi vs molecular transitions: y = gpa^ 

The Hi 21 cm hyperfine transition frequency is proportional to gpii~^a^^Roo (see I3.1.ip . On 
the other hand, the rotational transition frequencies of diatomic are inversely proportional to their 
reduced mass M. As on the example of Eq. (j34|l where we compared an electronic transition 
to a vibro-rotational transition, the comparison of the hyperfine and rotational frequencies is 
proportional to 

fht 2 M 2 

— « ffp^EM— - 5pa™ = 

where the variation of M/nip is usually suppressed by a large factor of the order of the ratio 
between the proton mass and nucleon binding energy in nuclei, so that we can safely neglect it. 

The constraint on the variation of y is directly determined by comparing the redshift as deter- 
mined from Hi and molecular absorption lines. 

Ay _ Zmol - Zh 
y 1 + Zmol 

This method was first applied [508] to the CO molecular absorption lines [5321 towards PKS 1413+135 
to get 

Ay/y = (-4±6) X 10"^ z = 0.247. 

The most recent constraint |373j relies on the comparison of the published redshifts of two absorp- 
tion systems determined both from Hi and molecular absorption. The first is a system at z — 0.6847 
in the direction ofTXS 0218+357 for which th spectra of C0(l-2), i3CO(l-2), 01^0(1-2), CO(2-3), 
HC0+(l-2) and HCN(l-2) are available. They concluded that 

Ay/y= (-0.16 + 0.54) X 10"^ z = 0.6847. (84) 

The second system is an absorption system in direction of PKS 1413+135 for which the molecular 
lines of C0(1 - 2), HC0+(l-2) and HCO+(2-3) have been detected. The analysis led to 

A?//2/= (-0.2 + 0.44) X 10"^ z = 0.247. (85) 
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Ref. [12 obtains the constraints |A?//y| < 3.4 x 10~^ at z ~ 0.25 and z ^ 0.685. 

The radio domain has the advantage of heterodyne techniques, with a spectral resohition of 
10^ or more, and deahng with cold gas and narrow lines. The main systematics is the kinematical 
bias, i.e. that the different lines do not come exactly from the same material along the line of 
sight, with the same velocity. To improve this method one needs to find more sources, which may 
be possible with the radiotelescope ALMA 0. 

3.4.7 OH - 18 cm: F = .gplaL^)''" 

Using transitions originating from a single species, like with SIDAM, allows to reduce the systematic 
effects. The f8 cm lines of the OH radical offers such a possibility [94l 1277] . 

The ground state, ^113/2^' = 3/2, of OH is split into two levels by A-doubling and each of 
these doubled level is further split into two hyperfine-structure states. Thus it has 2 "main" 
lines (AF = 0) and two "satellite" lines (AF = f). Since this four lines arise from two different 
physical processes (A-doubling and hyperfine splitting), they enjoy the same Rydberg dependence 
but different gp and dependencies. By comparing the four transitions to the Hi hyperfine line, 
one can have access to 

^^5pKmM)'-'' (86) 

and it was also proposed to combine them with HCO"*" transitions to lift the degeneracy. 

Using the four 18 cm OH lines from the gravitational lens at z ^ 0.765 toward PMN J0134-0931 
and comparing the Hi 21cm and OH absorption redshifts of the different components allowed to 
set the constraint [279] 

AF/F = (-0.44 ±0.36 ± l.Osyst) X 10-^ z = 0.765, (87) 

where the second error is due to velocity offsets between OH and Hi assuming a velocity dispersion 
of 3 km/s. A similar analysis [13_7^ in a system in the direction of PKS 1413+135 gave 

AF/F:= (0.51 ±1.26) X 10"^ z = 0.2467. (88) 

3.4.8 Far infrared fine-structure lines: F' = a 

Another combination [297| of constants can be obtained from the comparison of far infrared fine- 
structure spectra with rotational transitions, which respectively behaves as RooCt^^ and i?ooM — 
^00 /m so that they give access to 

A good candidate for the rotational lines is CO since it is the second most abundant molecule in 
the Universe after H2. 

Using the Cii fine-structure and CO rotational emission lines from the quasars Jl 148+5251 and 
BR 1202-0725, it was concluded that 

AF'/F' (0.1± 1.0) X 10-^, z = 6.42, (89) 
AF'/F' = (1.4 ± 1.5) X 10"^ z = 4.69, (90) 

which represents the best constraints at high redshift. As usual, when comparing the frequencies of 
two different species, one must account for random Doppler shifts caused by non-identical spatial 
distributions of the two species. Several other candidates for microwave and FIR lines with good 
sensitivities are discussed in Ref. [298] . 

^ http: / / www.eso.org/ sci / facilities/alma/ 
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3.4.9 "Conjugate" satellite OH lines: G = 3p(«em/^) 



1.85 



The satellite OH IScm lines are conjugate so that the two lines have the same shape, but with 
one line in emission and the other in absorption. This arises due to an inversion of the level 
of populations within the ground state of the OH molecule. This behaviour has recently been 
discovered at cosmological distances and it was shown |94j that a comparison between the sum 
and difference of satellite line redshifts probes G — gp{a^^fi)^-^^ . 

From the analysis of the two conjugate satellite OH systems at z ~ 0.247 towards PKS 1413+135 
and at z - 0.765 towards PMN J0134-0931, it was concluded [94 that 

\AG/G\ < 1.1 X 10-5. (91) 

It was also applied to a nearby system, Centaurus A, to give \AG/G\ < 1.6 x 10"^ at z ~ 0.0018. 
A more recent analysis [27 8\ claims for a tentative evidence (with 2.6cr significance, or at 99.1% 
confidence) for a smaller value of G 

AG/G = (-1.18 ±0.46) X 10^5 (92) 

for the system at z 0.247 towards PKS 1413+135. 

One strength of this method is that it guarantees that the satellite lines arise from the same 
gas, preventing from velocity offset between the lines. Also, the shape of the two lines must agree 
if they arise from the same gas. 

3.4.10 Molecular spectra and the electron-to-proton mass ratio 

As was pointed out in § 13. H molecular lines can provide a test of the variatiorH |482j of /i since 
rotational and vibrational transitions are respectively inversely proportional to their reduce mass 
and its square-root [see Eq. ([Ml) ]. 

Constraints with H2 

H2 is the most abundant molecule in the universe and there were many attempts to use its 
absorption spectra to put constraints on the time variation of fj, despite the fact that H2 is very 
difficult to detect [383] . 

As proposed in Ref. [507] . the sensitivity of a vibro-rotational wavelength to a variation of fi 
can be parameterized as 

A, = A0(1 + Zabs) (^1 + K,^y 

where A° is the laboratory wavelength (in the vacuum) and A^ is the wavelength of the transition 
i in the rest-frame of the cloud, that is at a redshift Zabs so that the observed wavelength is 
Ai/(1 + Zabs). Ki is a sensitivity coefficient analogous to the g-coefficient introduced in Eq. (|69|). 
but with different normalisation since in the parameterization we would have qi — uj^Ki/2, 

^ _ din A, 
* dln/i 

corresponding to the Lyman and Werner bands of molecular hydrogen. From this expression, one 
can deduce that the observed redshift measured from the transition i is simply 

= Zabs + fei^i, 6 = -(1 + Zabs)— , 



Again, ji is used either from me/mp or mp/rrie. I have chosen to use /i = mp/rrie and fj, = nie/mp. 
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which implies in particular that Zabs is not the mean of the Zi if 7^ . Indeed Zi is measured 
with some uncertainty of the astronomical measurements Xi and by errors of the laboratory mea- 
surements A^. But if A/i 7^ there must exist a correlation between Zi and Ki so that a linear 
regression of Zi (measurement) as a function of Ki (computed) allows to extract (^abs, b) and their 
statistical significance. 

We refer to § V.C of FVC [495] for earlier studies and we focus on the latest results. The 
recent constraints are mainly based on the molecular hydrogen of two damped Lyman-a absorption 
systems at z = 2.3377 and 3.0249 in the direction of two quasars (Q 1232+082 and Q 0347-382) 
for which a first analysis of VLT/UVES data showed |260J a slight indication of a variation, 

A/x//^ = (5.7±3.8) X 10"^ 

at 1.5ct for the combined analysis. The lines were selected so that they are isolated, unsaturated 
and unblended. It follows that the analysis relies on 12 lines (over 50 detected) for the first 
quasar and 18 (over 80) for second but the two selected spectra had no transition in common. 
The authors performed their analysis with two laboratory catalogs and got different results. They 
point out that the errors on the laboratory wavelengths are comparable to those of the astronomical 
measurements . 

It was further improved with an analysis on two absorption systems at z = 2.5947 and z — 
3.0249 in the directions of Q 0405-443 and Q 0347-383 observed with the VLT/UVES spectrograph. 
The data have a resolution R — 53000 and a S/N ratio ranging between 30 and 70. The same 
selection criteria where applied, letting respectively 39 (out of 40) and 37 (out of 42) lines for each 
spectrum and only 7 transitions in common. The combined analysis of the two systems led }261J 

A^//i= (1.65 ±0.74) X lO"'^ or A/i//i = (3.05 ± 0.75) x 10"^ 

according to the laboratory measurements that were used. The same data were reanalyzed with 
new and highly accurate measurements of the Lyman bands of H2, which implied a reevaluation 
of the sensitivity coefficient Ki . It leads to the two constraints 1425) 

A^/^ = (2.78 ±0.88) X 10"^ z = 2.59, (93) 
A/i//^= (2.06 ±0.79) X 10"^ z = 3.02, (94) 

leading to a 3.5cr detection for the weighted mean A/i//x — (2.4 ± 0.66) x 10~^. The authors of 
Ref. [425] do not claim for a detection and are cautious enough to state that systematics dominate 
the measurements. The data of the z = 3.02 absorption system were reanalyzed in Ref. |525) which 
claim that they lead to the bound |A/^//^| < 4.9 x 10"^ at a 2cr level, instead of Eq. ([M| . Adding 
a new set of 6 spectra, it was concluded that A/i//i = (15 ± 14) x 10^^ for the weighted fit |526] . 
These two systems were reanalyzed |289j . adding a new system in direction of Q 0528-250, 

A^//^ = (1.01 ±0.62) X 10-^ z = 2.59, (95) 
A/i//i== (0.82 ±0.74) X 10"^ z = 2.8, (96) 
A^/^ (0.26 ±0.30) X 10"^ z = 3.02, (97) 

respectively with 52, 68 and 64 lines. This gives a weighted mean of (2.6±3.0) x 10^^ at z ~ 2.81. To 
compare with the previous data, the analysis of the two quasars in common was performed by using 
the same lines (this implies adding 3 and removing 16 for Q 0405-443 and adding 4 and removing 
35 for Q 0347-383) to get respectively (-1.02 ± 0.89) x 10"^ (z = 2.59) and (-1.2 ± 1.4) x 10"^ 
(z = 3.02). Both analysis disagree and this latter analysis indicates a systematic shift of A/x//x 
toward 0. A second reanalysis of the same data was performed in Refs. |4841 1483] using a different 
analysis method to get 

A/i//^= (-7±8) X 10"^ (98) 
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Recently discovered molecular transitions at z = 2.059 toward the quasar J2123-0050 observed by 
the Keck telescope allow to obtain 86 H2 transitions and 7 HD transitions to conclude |340) 

A/^/^i = (5.6 ± 5.5stat ± 2.7syst) X 10-^ z = 2.059. (99) 

This method is subject to important systematic errors among which (1) the sensitivity to the 
laboratory wavelengths (since the use of two different catalogs yield different results [425] ). (2) 
the molecular lines are located in the Lyman-a forest where they can be strongly blended with 
intervening Hi Lyman-a absorption lines which requires a carfuU fitting of the lines |289j since it 
is hard to find lines that are not contaminated. From an observational point of view, very few 
damped Lyman-a systems have a measurable amount of H2 so that only a dozen systems is actually 
known even though more systems will be obtained soon ^405^ . To finish, the sensitivity coefficients 
are usually low, typically of the order of 10^^. Some advantages of using H2 arise from the fact 
there are several hundred available H2 lines so that many lines from the same ground state can be 
used to eliminate different kinematics between regions of different excitation temperatures. The 
overlap between Lyman and Werner bands also allow to reduce the errors of calibration. 

To conclude, the combination of all the existing observations indicate that fi is constant at the 
10~^ level during the past 11 Gigayrs while an improvement of a factor 10 can be expected in the 
five coming years. 

Other constraints 

It was recently proposed |202l 1203) that the inversion spectrum of ammonia allows for a better 
sensitivity to /i. The inversion vibro-rotational mode is described by a double well with the first two 
levels below the barrier. The tunnelling implies that these two levels are split in inversion doublets. 
It was concluded that the inversion transitions scale as i^inv ~ M^'^^i compared with a rotational 
transition which scales as frot ^ This implies that the redshifts determined by the two types of 
transitions are modified according to Sz-mv = 4.46(1 -I- Zabs)A/i//i and Sz-^ot ^ (1 + Zabs)A/i//i so 
that 

A/i/M- 0.289^i^^^— 

1 + Zabs 

Only one quasar absorption system, at z = 0.68466 in the direction of B 0218-1-357, displaying NH3 
is currently known and allows for this test. A first analysis |202j estimated from the published 
redshift uncertainties that a precision of ^ 2 x 10^^ on Afx/fj, can be achieved. A detailed mea- 
surement |368| of the ammonia inversion transitions by comparison to HCN and HCO"'" rotational 
transitions concluded that 

|A^/Ai| < 1.8 X 10"^ z = 0.685, (100) 

at a 2(7 level. Recently the analysis of the comparison of NII3 to IIC3N spectra was performed 
toward the gravitational lens system PKS 1830-211 (z ~ 0.89), which is a much more suitable 
system, with 10 detected NH3 inversion lines and a forest of rotational transitions. It reached the 
conclusion that 

|A/i//i| < 1.4 X 10"^ z = 0.89, (101) 

at a 3cr level [249] . From a comparison of the ammonia inversion lines with the NH3 rotational 
transitions, it was concluded [35 IJ 

|A/i//i| < 3.8 X 10-^ z = 0.89, (102) 

at 95% C.L. One strength of this analysis is to focus on lines arising from only one molecular species 
but it was mentionned that the frequencies of the inversion lines are about 25 times lower than 
the rotational ones, which might cause differences in the absorbed background radio continuum. 
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This method was also apphed 1319) in the Milky Way, in order to constrain the spatial vari- 
ation of ^ in the galaxy (see § I6.1.3P . Using ammonia emission lines from interstellar molecular 
clouds (Perseus molecular core, the Pipe nebula and the infrared dark clouds) it was concluded 
that A/x = (4— 14) X 10"^. This indicates a positive velocity offset between the ammonia inversion 
transition and rotational transitions of other molecules. Two systems being located toward the 
galactic center while one is in the direction of the anti-center, this may indicate a spatial variation 
of ^ on galactic scales. 

New possibilities 

The detection of several deuterated molecular hydrogen HD transitions makes it possible to 
test the variation of /i in the same way as with H2 but in a completely independent way, even 
though today it has been detected only in 2 places in the universe. The sensitivity coefficients have 
been published in Ref. |262) and HD was first detected by Ref. [383]. 

HD was recently detected [468j together with CO and H2 in a DLA cloud at a redshift of 2.418 
toward SDSS1439-I-11 with 5 lines of HD in 3 components together with several H2 lines in 7 
components. It allowed to set the ia limit of |A/i//i| < 9 x 10~^ . 4071 . 

Even though the small number of lines does not allow to reach the level of accuracy of H2 it is 
a very promising system in particular to obtain independent measurements. 

3.4.11 Emission spectra 

Similar analysis to constrain the time variation of the fundamental constants were also performed 
with emission spectra. Very few such estimates have been performed, since it is less sensitive and 
harder to extend to sources with high redshift. In particular, emission lines are usually broad as 
compared to absorption lines and the larger individual errors need to be beaten by large statistics. 
The Oiii doublet analysis [22] from a sample of 165 quasars from SDSS gave the constraint 

Aa™/«BM - (12 ± 7) X lO-l 0.16 < z < 0.8. (103) 

The method was then extended straightforwardly along the lines of the MM method and ap- 
plied [238] to the fine-structure transitions in Neiii, Nev, Oiii, Oi and Sii multiplets from a sample 
of 14 Seyfert 1.5 galaxies to derive the constraint 

Aa^^/a^^, = (150 ± 70) x 10"^ 0.035 < z < 0.281. (104) 

3.4.12 Conclusion and prospects 

This paragraph illustrates the diversity of methods and the progresses that have been achieved 
to set robust constraints on the variation of fundamental constants. Many systems are now used, 
giving access to different combinations of the constants. It exploits a large part of the electromag- 
netic spectrum from far infrared to ultra violet and radio bands and optical and radio techniques 
have played complementary roles. The most recent and accurate constraints are summarized in 
Table do] and Fig.H 

At the moment, only one analysis claims to have detected a variation of the fine structure 
constant (Keck/HIRES) while the VLT/UVES points toward no variation of the fine structure 
constant. It has led to the proposition that ajj,„ may be space dependent and exhibit a dipole, the 
origin of which is not explained. Needless to say that such a controversy and hypotheses are sane 
since it will help improve the analysis of this data but it is premature to conclude on the issue of 
this debate and the jury is still out. Most of the systematics have been investigated in details and 
now seem under control. 
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Constant 


Method 


System 


Constraint (xlO"'') 


Redshift 


Ref. 




AD 


21 


(-0.5 ±1.3) 


2.33 - 3.08 






AD 


15 


(-0.15 ±0.43) 


1.59 - 2.92 


m 




AD 


9 


(-3.09 ±8.46) 


1.19 - 1.84 


[MS] 




MM 


143 


(-0.57 ±0.11) 


0.2 - 4.2 


[370] 




MM 


21 


(0.01 ±0.15) 


0.4 - 2.3 


m 




SIDAM 


1 


(-0.012 ±0.179) 


1.15 


[555] 




SIDAM 


1 


(0.566 ± 0.267) 


1.84 


[358] 


y 


Hi - mol 


1 


(-0.16 ±0.54) 


0.6847 


[3731 




Hi - mol 


1 


(-0.2 ±0.44) 


0.247 


[373] 




CO, CHO+ 




(-4 ±6) 


0.247 


[532] 


F 


OH - Hi 


1 


(-0.44 ± 0.36 ±1.0syst) 


0.765 


[279] 




OH - Hi 


1 


(0.51 ± 1.26) 


0.2467 


mi 


X 


Hi - UV 


9 


(-0.63 ±0.99) 


0.23 -2.35 


[489] 




Hi - UV 


2 


-(0.17 ±0.17) 


3.174 


[355] 


F' 


Cii - CO 


1 


(1 ± 10) 


4.69 


[324] 




Cii - CO 


1 


(14± 15) 


6.42 


EMI 


G 


OH 


1 


< 1.1 


0.247, 0.765 


m 




OH 


1 


< 1.16 


0.0018 


m 




OH 


1 


(-1.18 ±0.46) 


0.247 


[278] 


M 


H2 


1 


(2.78 ±0.88) 


2.59 






H2 


1 


(2.06 ± 0.79) 


3.02 






H2 


1 


(1.01 ±0.62) 


2.59 


[289] 




H2 


1 


(0.82 ± 0.74) 


2.8 


[^ 




H2 


1 


(0.26 ±0.30) 


3.02 


[289] 




H2 


1 


(0.7 ±0.8) 


3.02, 2.59 


[484] 




NH3 


1 


< 0.18 


0.685 


[368] 




NH3 


1 


< 0.38 


0.685 


[351] 




HC3N 


1 


< 0.14 


0.89 


[249] 




HD 


1 


< 9 


2.418 


[iPT] 




HD 


1 


(0.56±0.55stat ±0.27syst) 


2.059 





Table 10: Summary of the latest constraints on the variation of fundamental constants obtained 
from the analysis of quasar absorption spectra. We recall that y = gpO^.^^, F = gp{a'^_^^iJ,)^'^'^ , 
^ = "emSp/M; F' = "L^ "''^^ l^ = "^p/"^e, G = gp{afi)^-^^. 



Let us what we can learn on the physics from these measurement. As an example, consider the 
constraints obtained on fj,, y and F in the redshift band 0.6-0.8 (see Table [TO]) . They can be used 
to extract independent constraints on gp, a^.^ and /i 

A/Vm = (0±0.18) X 10"^ Aa^^/a^^, = (-0.27±2.09) x 10"^ Affp/.gp = (0.38±4.73) x 10"^ 

This shows that one can test the compatibility of the constraints obtained from different kind 
of systems. Independently of these constraints, we have seen in § 16.31 that in grand unifica- 
tion theory the variation of the constants are correlated. The former constraints show that if 
Aln/Lt = i?Alnaj,„ then the constraint (|100p imposes that |i?AlnQ;jj,j^J < 1.8 x 10^^. In general 
R is expected to be of the order of 30 — 50. Even if its value its time-dependent, that would mean 
that Alna^j^j ^ (1 — 5) x 10~^ which is highly incompatible with the constraint ([73]) obtained by 
the same team on a^.^, but also on the constraints (|70|1 and ([7T|l obtained from the AD method 
and on which both teams agree. This illustrates how important the whole set of data is since one 
will probably be able to constrain the order of magnitude of i? in a near future, which would be a 
very important piece of information for the theoretical investigations. 
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Figure 4: Summary of the direct constraints on obtained from the AD (blue), MM (red) and 
AD (green) methods (left) and on fj, (right) that are summarized in Table fTOl 



We mention in the course of this paragraph many possibilities to improve these constraints. 

Since the AD method is free of the two main assumptions of the MM method, it seems important 
to increase the precision of this method as well as any method relying only on one species. This can 
be achieved by increasing the S/N ratio and spectral resolution of the data used or by increasing 
the sample size and including new transitions (e.g. cobalt |171l I188J ). 

The search for a better resolution is being investigated in many direction. With the current 
resolution of i? ~ 40000, the observed line positions can be determined with an accuracy of 
a\ ^ 1mA. This implies that the accuracy on Aq;j,m/Q!em of the order of 10~^ for lines with 
typical (/-coefficients. As we have seen this limit can be improved to 10~^ when more transitions or 
systems are used together. Any improvement is related to the possibility to measure line positions 
more accurately. This can be done by increasing R up to the point at which the narrowest lines in 
the absorption systems are resolved. The Bohlin formula [60| gives the estimates 

where AApix is the pixel size, Wobs is the observed equivalent width, is the mean number 
of photoelectron at the continuum level and M is the number of pixel covering the line profile. 
The metal lines have intrinsic width of a few km/s. One can thus expect improvements from 
higher spectral resolution. Progresses consercning the calibration are also expected, using e.g. 
laser comb [474]. Let us just mention, the EXPRESSO (Echelle Spectrograph for PREcision Super 
Stable Observation) project |114j on 4 VLT units or the CODEX (COsmic Dynamics EXplorer) on 
E-ELT projects [357|I355|I504) . They shall provide a resolving power of i? = 150000 to be compared 
to the HARP^ (High Accuracy Radial velocity planet Searcher) spectrograph (R ^ 112000) has 
been used but it is operating on a 3.6m telescope. 

The limitation may then lie in the statistics and the calibration and it would be useful to use 
more than two QSO with overlapping spectra to cross-calibrate the line positions. This means 
that one needs to discover more absorption systems suited for these analysis. Many progresses are 
expected. For instance, the FIR lines are expected to be observed by a new generation of telescopes 
such as HERSCHElfl. While the size of the radio sample is still small, surveys are being carried 

^http;/ /obswww.unige.ch/Instruments/HARPS/ 

*http:/ /sci.esa.int/science-e/www/area/index.cfm?fareaid=16 
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out so that the number of known redshift OH, HI and HCO+ absorption systems will increase. 
For instance the future Square Kilometre Array (SKA) will be able to detect relative changes of 
the order of 10"'' in a^^. 

In conclusion, it is clear that these constraints and the understanding of the absorption systems 
will increase in the coming years. 

3.5 Stellar constraints 

Stars start to accumulate helium produced by the pp-reaction and the CNO cycle in their core. 
Furthermore, the products of further nuclear reactions of helium with either helium or hydrogen 
lead to isotopes with A — 5 or A ^ 8, which are highly unstable. In order to produce elements 
heavier than A > 7 by fusion of lighter isotopes, the stars need to reach high temperatures and 
densities. In these conditions, newly produced ^^C would almost immediately be fused further to 
form heavier elements so that one expects only a tiny amount of ^^C to be produced, in contradic- 
tion with the observed abundances. This led Hoyle [255 to conclude that a then unknown excited 
state of the ^^C with an energy close to the 3a-threshold should exist since such a resonance would 
increase the probability that *Be captures an a-particlc. It follows that the production of ^^C in 
stars relies on the three conditions: 

• the decay lifetime of ®Be, of order 10^^^ s, is four orders of magnitude longer than the time 
for two a particles to scatter, so that a macroscopic amount of beryllium can be produced, 
which is sufficient to lead to considerable production of carbon, 

• an excited state of ^^C lies just above the energy of *Be + a, which allows for 

'^He + "^Hc o ®Be, ^Bc + ''He o ^^C* ^ ^^C + 7.367 MeV, 

• the energy level of ^^O at 7.1197 MeV is non resonant and below the energy of ^^C + a, 
of order 7.1616 MeV, which ensures that most of the carbon synthesized is not destroyed 
by the capture of an a-particle. The existence of this resonance, the O^-state of ^'^C was 
actually discovered |110] experimentally later, with an energy of 372 ± 4 keV [today, Eq+ — 
379.47 ± 0.15 keV], above the ground state of three a-particles (see Fig. [S|). 

The variation of any constant that would modify the energy of this resonance would also 
endanger the stellar nucleosynthesis of carbon, so that the possibility for carbon production has 
often been used in anthropic arguments. Qualitatively, if Eq+ is increased then the carbon would 
be rapidly processed to oxygen since the star would need to be hotter for the triple-a process to 
start. On the other hand, if Eq+ is decreased, then all a-particles would produce carbon so that no 
oxygen would be synthesized. It was estimated |332] that the carbon production in intermediate 
and massive stars is suppressed if the various of the energy of the resonance is outside the range 
-250 keV< AE(y+ < 60 keV, which was further improved [443, to, -5 keV< AE^^ < 50 keV 
in order for the C/0 ratio to be larger than the error in the standard yields by more than 50%. 
Indeed, in such an analysis, the energy of the resonance was changed by hand. We expect however 
that if E^+ is modified due to the variation of a constant other quantities, such as the resonance 
of the oxygen, the binding energies and the cross-sections will also be modified in a complex way. 

In practice, to draw a constraint on the variation of the fundamental constants from the stellar 
production of carbon, one needs to go through different steps, any of them involving assumptions, 

1. to determine the effective parameters, e.g. cross sections, which affects the stellar evolution. 
The simplest choice is to modify only the energy of the resonance but it may not be realistic 
since all cross-sections and binding energies should also be affected. This requires to use a 
stellar evolutionary model; 
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Figure 5: (Left) Level scheme of nuclei participating to the ^lie{aa,jy^C reaction. (Right) Central 
abundances at the end of the CHe burning as a function of 5nn for a GOM© star with Z = 0. From 
Ref m^ . 

2. relate these parameters to nuclear parameters. This involves the whole nuclear physics ma- 
chinery; 

3. to relate the nuclear parameters to fundamental constants. As for the Oklo phenomenon, it 
requires to link QCD to nuclear physics. 



A first analysis [3881 13891 1443) used a model that treats the carbon nucleus by solving the 
12-nucleon Schrodinger equation using a three-cluster wave-function representing the three-body 
dynamics of the state. The NN interaction was described by the Minnesota model |2961 1485] 
and its strength was modified by multiplying the effective NN-potential by an arbitrary number p. 
This allows to relate the energy of the Hoyle level relative to the triple alpha threshold, e = Qaaa, 
and the gamma width, F^, as a function of the parameter p, the latter being almost not affected. 
The modified 3Q!-reaction rate was then given by 



\MakBT J h 



F 



■ exp 



eipY 
fceT 



(105) 



where and Na are the mass and number density of the a-particle, The resonance width 
F = TaTry/ (Fc -l-F-y) ~ F^. Thls was included in a stellar code and ran for red giant stars with 1.3, 
5 and 2OM0 with solar metallicity up to thermally pulsating asymptotic giant branch [388] and in 
low, intermediate and high mass (1.3, 5, 15, 25M0) with solar metallicity also up to TP-AGB |443j 
to conclude that outside a window of respectively 0.5% and 4% of the values of the strong and 
electromagnetic forces, the stellar production of carbon or oxygen will be reduced by a factor 30 
to 1000. 

In order to compute the resonance energy of the beryllium-8 and carbon-12 a microsopic cluster 
model was developed |296] . The Hamiltonian of the system is then of the form H ^ J^T T{ri + 
J2j<i^i^ij)j where A is the nucleon number, T the kinetic energy and V the NN interaction 
potential. In order to implement the variation of the strength of the nuclear interaction with 
respect to the electromagnetic interaction, it was taken as 



V{rij) = Vb(ry) + (1 + 6NN)VNir^j), 
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where Snn is a dimensionless parameter that describes the change of the nuclear interaction, Vn 
being described in Ref. |485| . When A > 4 no exact solution can be found and approximate 
solutions in which the wave function of the beryllium-8 and carbon- 12 are described by clusters of 
respectively 2 and 3 a-particle is well adapted. 

First, Snn can be related to the deuterium binding energy as 

ABd/Bd^ 5.7701 X Snn (106) 

which, given the discussion in § 13.8.31 allows to relate Snn to fundamental constants, as e.g. in 
Ref. |104j . Then, the resonance energy of the beryllium-8 and carbon- 12 scale as 

Eni^Be) = (0.09208 - 12.208 x Snn) Mev, ErC^C) = (0.2877 - 20.412 x Snn) Mev, (107) 

so that the energy of the Hoyle level relative to the triple alpha threshold is Qaaa — Ei^{^Be) + 

This was implemenetd [10311180] to population III stars with typical masses, 15 and 60 Mq with 
zero metallicity, in order to compute the central abundances at the end of the core He burning. 
From Fig. [SJ one can distinguish 4 regimes (I) the star ends the CHe burning phase with a core 
composed of a mixture of carbon-12 and oxygen-16, as in the standard case; (II) if the 3a rate is 
weaker, ^^C is produced slower, the raction ^^C(a,7)^®0 becomes efficient earlier so that the star 
ends the CHe burning phase with a core composed mostly of oxygen-16; (HI) for weaker rates, the 
oxygen-16 is further processed to neon-20 and then ^''Mg so that the star ends the CHe burning 
phase with a core composed of ^''Mg and (IV) if the 3a rate is stronger, the carbon-12 is produced 
more rapidly and the star ends the CHe burning phase with a core composed mostly of carbon-12. 
Typically this imposes that 

- 5 X 10"* < (Jatat < 1.5 X 10"^, -ixlQ-^ < ABd/Bd <'dxlQ-^ (108) 

to ensure the ratio C/0 to be of order unity. 

To finish, a recent study |5] focus on the existence of stars themselves, by revisiting the stellar 
equilibrium when the values of some constants are modified. In some sense, it can be seen as a 
generalization of the work by Gamow [224_ to constrain the Dirac model of a varying gravitational 
constant by estimating its effect on the lifetime of the Sun. In this semi-analytical stellar structure 
model, the effect of the fundamental constants was reduced phenomenologically to 3 parameters, 
G which enters mainly on the hydrostatic equilibrium, a-^^^ which enters in the Coulomb barrier 
penetration through the Gamow energy, and a composite parameter C which describes globally 
the modification of the nuclear reaction rates. The underlying idea is to assume that the power 
generated per unit volume, e(r), and which determines the luminosity of the star, is proportional 
to the fudge factor C, which would arise from a modification of the nuclear fusion factor, or equiv- 
alently of the cross section. It thus assumes that all cross-sections are affected is a similar way. 
The parameter space for which stars can form and for which stable nuclear configurations exist 
was determined, showing that no fine-tuning seems to be required. 

This new system is very promising and will provide new informations on the fundamental 
constants at redshifts smaller than z ~ 15 where no constraints exist at the moment, even though 
drawing a robust constraint seems to be difficult at the moment. In particular, an underlying 
limitation arises from the fact that the composition of the interstellar media is a mixture of ejecta 
from stars with different masses and it is not clear which type of stars contribute the most the 
carbon and oxygen production. Besides, one would need to include rotation and mass loss }181] . 
As for the Oklo phenomenon, another limitation arises from the complexity of nuclear physics. 
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3.6 Cosmic Microwave Background 

The CMB radiation is composed of photons emitted at the time of the recombination of hydrogen 
and hehum when the universe was about 300,000 years old [see e.g. Ref. |404) for details on the 
physics of the CMB] . This radiation is observed to be a black-body with a temperature Tq = 2.723 K 
with small anisotropics of order of the /xK. The temperature fluctuation in a direction (??, 1^9) is 
usually decomposed on a basis of spherical harmonics as 

— (i?,^) =^ ^ a,„r,™(i?,^). (109) 

e m=-£ 

The angular power spectrum multipole — (|a;„iP) is the coefficient of the decomposition of the 
angular correlation function on Legendre polynomials. Given a model of structure formation and 
a set of cosmological parameters, this angular power spectrum can be computed and compared to 
observational data in order to constrain this set of parameters. 



The CMB temperature anisotropics mainly depend on three constants: G, a^^^ and mo. 

The gravitational constant enters in the Friedmann equation and in the evolution of the cosmo- 
logical perturbations. It has mainly three effects [429J that are detailed in § 14.4.11 a^.^ , rric affect 
the dynamics of the recombination. Their influence is complex and must be computed numerically. 
We can however trace their main effects since they mainly modify the CMB spectrum through the 
change in the differential optical depth of photons due to the Thomson scattering 

f = XcficCaT (110) 

which enters in the collision term of the Boltzmann equation describing the evolution of the photon 
distribution function and where Xc is the ionization fraction (i.e. the number density of free 
electrons with respect to their total number density no). 

The first dependence arises from the Thomson scattering cross-section given by 

CTt = TT^Ot (111) 

and the scattering by free protons can be neglected since rac/nip ~ 5 x lO"**. 

The second, and more subtle dependence, comes from the ionization fraction. Recombination 
proceeds via 2-photon emission from the 2s level or via the Ly-a photons which are redshifted 
out of the resonance line |400| because recombination to the ground state can be neglected since 
it leads to immediate reionization of another hydrogen atom by the emission of a Ly-a photons. 
Following Refs. [400|I336] and taking into account, for the sake of simplicity, only the recombination 
of hydrogen, the equation of evolution of the ionization fraction takes the form 



dxo _ 
At 



^ (1 - Xo) exp ( ^ ) - TZn^xl 



fll2) 



where Tm is the temperature. At high redshift, Tm is identical to the one of the photons T^ — 
To(l + z) but evolves according to 

where the radiation constant = 4ctsb/c with ctsb = fcg7r^/(607rc^fi,'^) the Stefan-Boltzmann 
constant. In Eq. (|112p . i?„ = —Eij-r? is the energy of the nth hydrogen atomic level, /? is the 
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ionization coefficient, TZ the recombination coefficient, C the correction constant due to tlie redstiift 
of Ly-a ptiotons and to 2-ptioton decay and Up — Ue is ttie number density of protons. /3 is related 
to TZ by tire principle of detailed balance so that 



The recombination rate to all other excited levels is 



7 JB,^ik^T ev - 1 

where cr„/ is the ionization cross-section for the {n, I) excited level of hydrogen. The star indicates 
that the sum needs to be regularized and the a^.^-, mo-dependence of the ionization cross-section 
is complicated to extract. It can however be shown to behave as f7„; oc f {hv / Bi) . Finally, 

the factor C is given by 

l + KK2s{l-x,) 

l + K{P + h.2s){l-X,) ^ ^' 

where is the rate of decay of the 2s excited level to the ground state via 2 photons; it scales 
as meCt^j^. The constant K is given in terms of the Ly-a photon Aq = IGnh/ (SmcO^-^c) by 
K = npX^/{8TrH) and scales as m~^a~^. 

In summary, both the temperature of the decoupling and the residual ionization after recom- 
bination are modified by a variation of a^.^, or mc. This was first discussed in Ref. 1281) . 
The last scattering surface can roughly be determined by the maximum of the visibility function 
g = f exp(— r) which measures the differential probability for a photon to be scattered at a given 
redshift. Increasing Q!j,„ shifts g to a higher redshift at which the expansion rate is faster so that 
the temperature and Xe decrease more rapidly, resulting in a narrower g. This induces a shift of 
the Ce spectrum to higher multipoles and an increase of the values of the Ce- The first effect can 
be understood by the fact that pushing the last scattering surface to a higher redshift leads to a 
smaller sound horizon at decoupling. The second effect results from a smaller Silk damping. 

Most studies have introduced those modification in the RECFAST code ^450j including similar 
equations for the recombination of helium. Our previous analysis shows that the dependences in 
the fundamental constants have various origins, since the binding energies Bi scale has mcO^^, 
(It as m~^, K as m~^a~^, the ionisation coefficients B a,s , the transition frequencies as 

^ EM G ' e EM ' EM ' ^ 

uiM^ , the Einstein coefficients as TOpQ;^ , the decay rates A as and TZ has complicated 

EA4 ' ^ EM ' ^ EM ^ 

dependence which roughly reduces to a~^m~'^. Note that a change in the fine-structure constant 
and in the mass of the electron are degenerate according to Aa^.^ « 0.39Ame but this degeneracy 
is broken for multipoles higher than 1500 [3S]. In earlier works |243[ 1281] it was approximated by 
the scaling TZ cx a^^-^-^^^ with f ^ 0.7. 

o EM ^ 

The first studies |243} I281j focused on the sensibility than can be reached by WMAfU and 
PlanclH. They concluded that they should provide a constraint on 0;^^, at recombination, i.e. at 
a redshift of about z ~ 1,000, with a typical precision lAag^/a^jJ ~ 10^^ — 10^'^. 

The first attempt [TO] to actually set a constraint was performed on the first release of the data 
by BOOMERANG and MAXIMA. It concluded that a value of a^.^, smaller by a few percents in the 
past was favoured but no definite bound was obtained, mainly due to the degeneracies with other 
cosmological parameters. It was later improved [21] by a joint analysis of BBN and CMB data 
that assumes that only a^-^j^ varies and that included 4 cosmological parameters {Qriia.t,^h,h,ns) 



^http: //map. gsfc .nasa.gov/ 

^http : //astro . estec . esa. nl/SA-general/Projects/Planck/ 
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assuming a universe with Euclidean spatial section, leading to —0.09 < Aa^M < 0-02 at 68% 
confidence level. A similar analysis |307| . describing the dependence of a variation of the fine- 
structure constant as an effect on recombination the redshift of which was modelled to scale as 
2* = 1080[1 + 2Aq;j,m/Q!bm]: set the constraint —0.14 < Aa,^^ < 0.02, at a 2cr level, assuming 
a spatially flat cosmological models with adiabatic primordial fluctuations that. The effect of 
reionisation was discussed in Ref. [347]. These works assume that only a^.^ is varying but, as can 
been seen from Eqs. (|109mi5|) . assuming the electron mass constant. 

With the WMAP first year data, the bound on the variation of a^.^, was sharpened [432] 
to —0.05 < Aajj,^, /ofjj.^ < 0.02, after marginalizing over the remaining cosmological parameters 
{^ma.th^ ,^hh^ ,^h^ ,ns, as,T) assuming a universe with Euclidean spatial sections. Restricting to 
a model with a vanishing running of the spectral index (as = dns/dlnk = 0), it gives —0.06 < 
Aaj,„/Q!j,„ < 0.01, at a 95% confidence level. In particular it shows that a lower value of 
makes as = more compatible with the data. This bounds were obtained without using other 
cosmological data sets. This constraint was confirmed by the analysis of Ref. [259] . which got 
-0.097 < Aa^MaEM < 0-034, with the WMAP-lyr data alone and -0.042 < Aa^„/a^„ < 0.026, 
at a 95% confidence level, when combined with constraints on the Hubble parameter from the HST 
Hubble Key project. 

The analysis of the WMAP-3yr data allows to improve |472] this bound to —0.039 < Aa^.^ /q^em < 
0.010, at a 95% confidence level, assuming (r^mat, ^h, h, Us, Zro, As) for the cosmological parameters 
{Q\ being derived from the assumption fix ~ 0, as well as r from the reionisation redshift, Zro) 
and using both temperature and polarisation data (TT, TE, EE). 

The WMAP 5-year data were analyzed, in combination with the 2dF galaxy redshift survey, 
assuming that both ofj^.^, and can vary and that the universe was spatially Euclidean. Letting 6 
cosmological parameters [{flmath'^ , ^hh^, ©, i^s, As), Q being the ratio between the sound horizon 
and the angular distance at decoupling] and 2 constants vary they, it was concluded [4451 1446] 
-0.012 < Aa^M/aEM < 0.018 and -0.068 < Am^/mc < 0.044, the bounds fluctuating slightly 
depending on the choice of the recombination scenario. A similar analyis [376' not including rrio 
gave —0.050 < Aa^^/a^,^ < 0.042, which can be reduced by taking into account some further prior 
from the HST data. Including polarisation data data from ACBAR, QUAD and BICEP, it was 
also obtained [350] -0.043 < Aa^„/a^„ < 0.038 at 95% C.L. and -0.013 < Aa^^Ja^^, < 0.015 
including HST data, also at 95% C.L. Let us also emphasize the work by Ref. [3.48. trying to 
include the variation of the Newton constant by assuming that Aa^^/aEM = QAG/G, Q being a 
constant and the investigation of Ref. [377] taking into account q;em j "^c and /i, G being kept fixed. 
Considering {ilyaa,t,^h, h,ns,T) for the cosmological parameters they concluded from WMAP-5 
data (TT, TE, EE) that -8.28 x 10"^ < Aa^„/a^„ < 1.81 x lO"'"^ and -0.52 < A/i//i < 0.17 

The analysis of Refs. |445[|446] was updated [305] to the WMAP-7yr data, including polarisation 
and SDSS data. It leads to -0.025 < Aa^^/a^„ < -0.003 and 0.009 < Arric/mc < 0.079 at a la- 
level. 

The main limitation of these analysis lies in the fact that the CMB angular power spectrum 
depends on the evolution of both the background spacetime and the cosmological perturbations. It 
follows that it depends on the whole set of cosmological parameters as well as on initial conditions, 
that is on the shape of the initial power spectrum, so that the results will always be conditional 
to the model of structure formation. The constraints on q;em or rric can then be seen mostly 
as constraints on a delayed recombination. A strong constraint on the variation of a-^-^^ can be 
obtained from the CMB only if the cosmological parameters are independently known. Ref. [432] 
forecasts that CMB alone can determine a^^ to a maximum accuracy of 0.1%. 
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Constraint 
(Qem X 10^) 


Data 


Comment 




Rcf. 


[-9,2] 


BOOMERanG-DASI-COBE + BBN 


BBN with Qj,„ only 




mi 






(finiat, ^b, h, Us) 






[-1.4, 2] 


COBE-BOOMERanG-MAXIMA 


(r^mat, f2b, h, Us) 




[307| 


[-5,2] 


WMAP-1 


(f^mat/l^, ilbh^,^Ah^,T, Us 


as) 


[432| 


[-6, 1] 


WMAP-1 


same + 0^—0 




,432] 


[-9.7, 3.4] 


WMAP-1 


(S7mat, ^b, h, ns,T, rric) 




[259] 


[-4.2,2.6] 


WMAP-1 + HST 


same 




[259] 


[-3.9, 1.0] 


WMAP-3 (TT,TE,EE) + HST 


(^mat , ^^b , ^, , ^rc, ) 






[-1.2,1.8] 


WMAP-5 + ACBAR + CBI + 2df 


(f^mat/l^, fib/l^, 8, T, Us, As 


nie) 


p5l 


[-1.9,1.7] 


WMAP-5 + ACBAR + CBI + 2df 


(fimat/l^, fib/l^, 8, T, Us, As 


nie) 


[3461 


[-5.0,4.2] 


WMAP-5 + HST 


(f^mat/l^, fibft^, h, T, Us, As) 




[575| 


[-4.3, 3.8] 


WMAP-5 + ACBAR + QUAD + BICEP 


(r^mat/l^, fib/l^, h, T, Us) 




[350] 


[-1.3,1.5] 


WMAP-5 + ACBAR + QUAD + BICEP -f HST 


(f^mat/l^, fibft^, h, T, Us) 




[350] 


[-0.83,0.18] 


WMAP-5 (TT,TE,EE) 


(f^mat/l^, ilbh^,h, T, Us, As, 


rric, fi) 


[577| 


[-2.5, -0.3] 


WMAP-7 + Ho + SDSS 


(^niat/l^, ribft^, 8, r. Us, As 


nic) 


'305] 



Table 11: Summary of the latest constraints on the variation of fundamental constants obtained 
from the analysis of cosmological data and more particularly of CMB data. All assume Q.k = 0. 



3.7 21 cm 

After recombination, the CMB photons are redshifted and their temperature drops as (1 -|- z). 
The baryons however are prevented from cooUng adiabatically since the residual amount of free 
electrons, that can couple the gas to the radiation through Compton scattering, is too small. It 
follows that the matter decouples thermally from the radiation at a redshift of order z ~ 200. 

The intergalactic hydrogen atoms after recombination are in their ground state which hyperfine- 
structure splits into a singlet and a triple states (lsi/2 with F ~ Q and F = 1 respectively, see 
§ III.B.l of FCV [495] ). It was recently proposed [284) that the observation of the 21 cm emission 
can provide a test on the fundamental constants. We refer to Ref. [285] for a detailed review on 
21 cm. 

The fraction of atoms in the excited (triplet) state versus the ground (singlet) state is conven- 
tionally related by the spin temperature Tg defined by the relation 

3expf-^) (116) 



where = hc/{X2ik^) = 68.2 mK is the temperature corresponding to the 21 cm transition and 
the factor 3 accounts for the degeneracy of the triplet state (note that this is a very simplified 
description since the assumption of a unique spin temperature is probably not correct [285j . The 
population of the two states is determined by two processes, the radiative interaction with CMB 
photons with a wavelength of A21 = 21.1 cm (i.e. 1^21 = 1420 MHz) and spin-changing atomic 
collision. The evolution of the spin temperature is thus dictated by 

The first term corresponds to the collision desexcitation rate from triplet to singlet and the coeffi- 
cient Cio is decomposed as 
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with the respective contribution of H-H and e-H colhsions. The second term corresponds to spon- 
taneous transition and Aio is the Einstein coefficient. The equation of evohition for the gas tem- 
perature Tg is given by Eq. (jll3p with Tm = Tg (we recaU that we have neglected the contribution 
of heUum) and the electronic density satisfies Eq. (|112l) . 

It follows [284] that the change in the brightness temperature of the CMB at the corresponding 
wavelength scales as Tb oc Ai2/j^|]^. Observationally, we can deduce the brightness temperature 
from the brightness 7^, that is the energy received in a given direction per unit area, solid angle 
and time, defined as the temperature of the black-body radiation with spectrum I^. Thus fceTb — 
jlv^ . It has a mean value, T\,{zoh&) at various redshift where l-|-Zobs = 1^21'^'*^ /^obs- Besides, as 
for the CMB, there will also be fiuctuation in Tb due to imprints of the cosmological perturbations 
on Up and Tg. It follows that we also have access to an angular power spectrum C£(zobs) at various 
redshift (see Ref. |327] for details on this computation). 

Both quantities depend on the value of the fundamental constants. Beside the same depen- 
dencies of the CMB that arise from the Thomson scattering cross-section, we have to consider 
those arising from the collision terms. In natural units, the Einstein coefficient scaling is given by 
^12 = I'TaEM'^li'^^^ ~ 2.869 X 10~^^s~-^. It follows that it scales as Axq cx ^p^^a^^TOe- The 
brightness temperature depends on the fundamental constant as Tb cx g^yutx^^^jm^. Note that the 
signal can also be affected by a time variation of the gravitational constant through the expansion 
history of the universe. Ref. [284! (see also Ref. |285j for further discussions), focusing only on 
Ofg^ , showed that this was the dominant effect on a variation of the fundamental constant (the 
effect on Cio is much complicated to determine but was argued to be much smaller). It was esti- 
mated that a single station telescope like LWA0 or LOFAlcl can lead to a constraint of the order 
of Aq!j,„/q!e^5 ^ 0.85%, improving to 0.3% for the full LWA. The fundamental challenge for such 
a measurement is the substraction of the foreground. 

The 21 cm absorption signal in a available on a band of redshift typically ranging from z < 1000 
to z '-^ 20, which is between the CMB observation and the formation of the first stars, that is during 
the so called "dark age". It thus offers an interesting possibility to trace the constraints on the 
evolution of the fundamental constants between the CMB epoch and the quasar absorption spectra. 

As for CMB, cosmological parameters since a change of 1% in respectively the baryon density 
or the Hubble parameter implies a 2% (resp. 3%)on the mean bolometric temperature. The effect 
on the angular power spectrum have been estimated but still require an in depth analysis along 
the lines of e.g. [327. It is motivating since Ci[zo\^s) is expected to depend on the correlators of 
the fundamental constants e.g. (a^^, (x, Zobs)aEM (^'i -^obs)) and thus in principle allows to study 
their fluctuation, even though it will also depend on the initial condition, e.g. power spectrum, of 
the cosmological perturbations. 

In conclusion, the 21cm observation opens a observational window on the fundamental at 
redshifts ranging typically from 30 to 100, but full in-depth analysis is still required (see Refs. |205| 
1286) for a critical discussion of this probe). 

3.8 Big bang nucleosynthesis 
3.8.1 Overview 

The amount of *IIc produced during the big bang nucleosynthesis is mainly determined by the 
neutron to proton ratio at the freeze-out of the weak interactions that interconvert neutrons and 
protons. The result of Big Bang nucleosynthesis (BBN) thus depends on G, a^, a^.^ and 
respectively through the expansion rate, the neutron to proton ratio, the neutron-proton mass 

^http:/ /Iwa. unm.edu 
*http:/ /www. lofar.org 
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difference and the nuclear reaction rates, besides the standard parameters such as e.g. the number 
of neutrino famihes. 

The standard BBN scenario |116[ 1404) proceeds in three main steps: 

1. for T > 1 MeV, (t < 1 s) a first stage during which the neutrons, protons, electrons, positrons 
an neutrinos are kept in statistical equilibrium by the (rapid) weak interaction 

n< — yp + e^+Pe, n-\-Ve< — !-p + e~, n + e^i — >p + Ve. (118) 

As long as statistical equilibrium holds, the neutron to proton ratio is 

[n/p) = c-'^npAe^ (119) 

where Qnp = — mp)c'^ = f .29 MeV. The abundance of the other light elements is given 



by [404] 

Ya = QA 



A-l ^ k T l3('4-l)/2 



( C(3) \ 9(3A-5)/2 ^5/2 



77- iy/y/-^e^^/'^Br, (120) 

where gA is the number of degrees of freedom of the nucleus ;^X, is the nucleon mass, 7] 
the baryon-photon ratio and Ba = {Zra^ + (A — Z)ra-a — raA)c^ the binding energy. 

2. Around T ^ 0.8 MeV [t ^ 2 s), the weak interactions freeze out at a temperature Tf 
determined by the competition between the weak interaction rates and the expansion rate of 
the universe and thus roughly determined by r^(Tf) ~ H{T{) that is 



GliKTt)' y/GKiKTt)^ (121) 

where Gp is the Fermi constant and the number of relativistic degrees of freedom at Tf. 
Below Tf, the number of neutrons and protons change only from the neutron /3-decay between 
Tf to Tn ~ 0.1 MeV when p + n reactions proceed faster than their inverse dissociation. 

3. For 0.05 MeV< T < 0.6 MeV (3s < i < 6min), the synthesis of light elements occurs only 
by two-body reactions. This requires the deuteron to be synthesized (p + n — >■ D) and the 
photon density must be low enough for the photo-dissociation to be negligible. This happens 
roughly when 

— ~r?2exp(-Bz3/rN)-l (122) 
with - 3 X 10-1". The abundance of '^He by mass, Fp, is then well estimated by 

- (123) 

1 + [n/p)^ 

with 

{n/p)^ = {n/p)f exp(-tN/T„) (124) 

with tN oc and t'^ = 1.6366^(1 + 3g\)ml / {2tt^) , with gA 1.26 being the 

axial/vector coupling of the nucleon. Assuming that Bjj oc a^, this gives a dependence 

4. The abundances of the light element abundances, Yi, are then obtained by solving a series of 
nuclear reactions 

Y, = j- ry„ 

where J and F are time-dependent source and sink terms. 
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From an observational point of view, the light elements abundances can be computed as a 
function of rj and compared to their observed abundances. Fig. [S] summarizes the observational 
constraints obtained on helium-4, helium-3, deuterium and lithium-7. On the other hand, rj can 
be determined independently from the analysis of the cosmic microwave background anisotropics 
and the WMAP data [295] have led to to the conclusion that 

'y^Wp = (6-14 ±0.25) X 10-1°. 

This number being fixed, all abundances can be computed. At present, there exists a discrepancy 
between the predicted abundance of lithium-7 based on the WMAP results [10711102) for rj, ''Li/H = 
(5.14±0.50) X IQ-i" and its values measured in metal-poor halo stars in our Galaxy [61], ^Li/H — 
(1.26 ±0.26) X IQ-i" which is factor 3 lower, at least |115| (see also Ref. [465] ). than the predicted 
value. No solution to this Lithium-7 problem is known. A back of the envelope estimates shows 
that we can mimic a lower 77 parameter, just by modifying the deuterium binding energy, letting 
Tn unchanged, since from Eq. (|122p . one just need ABd/T^ ~ —In 9 so that the effective rj 
parameter, assuming no variation of constant, is three times smaller than jy^vMAP- This rough rule 
of thumb explains that the solution of the lithium-7 problem may lie in a possible variation of the 
fundamental constants (see below for details). 

3.8.2 Constants everywhere... 

In complete generality, the effect of varying constants on the BBN predictions is difficult to model 
because of the intricate structure of QCD and its role in low energy nuclear reactions. A solution is 
thus to proceed in two steps, first by determining the dependencies of the light element abundances 
on the BBN parameters and then by relating those parameters to the fundamental constants. 

The analysis of the previous section, that was restricted to the helium-4 case, clearly shows that 
the abundances will depend on: (1) which will affect the Hubble expansion rate at the time of 
nucleosynthesis in the same way as extra-relativistic degrees of freedom do, so that it modifies the 
freeze-out time Tf. This is the only gravitational sector parameter. (2) Tn, the neutron lifetime 
dictates the free neutron decay and appears in the normalisation of the proton-neutron reaction 
rates. It is the only weak interaction parameter and it is related to the Fermi constant Gp, or 
equivalently the Higgs vev. (3) a^^, the fine-structure constant. It enters in the Coulomb barriers 
of the reaction rates through the Gamow factor, in all the binding energies. (4) Qnp, the neutron- 
proton mass difference enters in the neutron-proton ratio and we also have a dependence in (5) 
rriN and rUe and (6) the binding energies. 

Clearly all these parameters are not independent but their relation is often model-dependent. 
If we focus on helium-4, its abundance mainly depends on Qnp, Tf and Tn (and hence mainly 
on the neutron hfetime, Tn). Early studies (see § III.C.2 of FVC [495] ) generally focused on one 
of these parameters. For instance, Kolb et al. [294] calculated the dependence of primordial ^He 
on G, Gp and Qnp to deduce that the helium-4 abundance was mostly sensitive in the change in 
Qnp and that other abundances were less sensitive to the value of Qnp, mainly because ^He has 
a larger binding energy; its abundances is less sensitive to the weak reaction rate and more to 
the parameters fixing the value of (n/p). To extract the constraint on the fine-structure constant, 
they decomposed Qnp as Qnp = ap^Qa + PQp where the first term represents the electromagnetic 
contribution and the second part corresponds to all non-electromagnetic contributions. Assuming 
that Qa and Q^ are constant and that the electromagnetic contribution is the dominant part of 
Q, they deduced that |Aap„/ap„| < 10~^. Campbell and Olive ^ kept track of the changes in 
Tf and Qnp separately and deduced that ~ — while more recently the analysis [309] 

focused on a^^ and v. 

Let us now see how the effect of all these parameters are now accounted for in BBN codes. 
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Bergstrom et al. [50] started to focus on the a^.^, -dependence of the thermonuclear rates. In 
the non-relativistic hmit, it is obtained as the thermal average of the product of the cross, the 
relative velocity and the the number densities. Charged particles must tunnel through a Coulomb 
barrier to react. Changing a-,^,^ modifies these barriers and thus the reaction rates. Separating the 
Coulomb part, the low-energy cross section can be written as 

a{E) = ^e-^-^iE) (^25) 

where rj{E) arises from the Coulomb barrier and is given in terms of the charges and the reduced 
mass Mr of the two interacting particles as 



n{E)^a^,,Z,Z2sj^. (126) 

The form factor S{E) has to be extrapolated from experimental nuclear data but its a j,^, -dependence 
as well as the one of the reduced mass were neglected. Keeping all other constants fixed, as- 
suming no exotic effects and taking a lifetime of 886.7 s for the neutron, it was deduced that 
|Aq;j,„/q;j,^J < 2 X 10~^. This analysis was then extended [381] to take into account the a^^- 
dependence of the form factor to conclude that 



Ref. [381] also took into a account (1) the effect that when two charged particles are produced 
they must escape the Coulomb barrier. This effect is generally weak because the Qi-values (energy 
release) of the different reactions are generally larger than the Coulomb barrier at the exception 
of two cases, '^He(ri,p)'^H and ^Be(n,p)^Li. The rate of these reactions must be multiplied by a 
factor (1 -f aiAaj,„/aj.„). (2) The radiative capture (photon emitting processes) are proportional 
to a^j, since it is the strength of the coupling of the photon and nuclear currents. All these rates 
need to be multiplied by (1 -I- lS.a-^^^ / a-^„) . (3) The electromagnetic contribution to all masses 
was taken into account, which modify the Q^-values as Qi ^ Qi + g^Aa^^/aj^^j). For helium-4 
abundance these effects are negligible since the main ctj^.^ -dependence arises from Q^p- Equiped 
with these modifications, it was concluded that Aa^-^^/a^-^ — — O.OOT^q q]^" using only deuterium 
and helium-4 since the lithium-7 problem was still present. 

Then the focus fell on the deuterium binding energy, Bjj. Flambaum and Shuryak |206[ 12071 
11571 1156] illustrated the sensitivity of the light element abundances on Bd ■ Its value mainly sets 
the beginning of the nucleosynthesis, that is of Tn since the temperature must low-enough in order 
for the photo-dissociation of the deuterium to be negligible (this is at the origin of the deuterium 
bottleneck). The importance of Bu is easily understood by the fact that the equilibrium abundance 
of deuterium and the reaction rate p(n,7)D depends exponentially on Bjj and on the fact that 
the deuterium is in a shallow bound state. Focusing on the TN-dependence, it was concluded |206] 
that ABd/Bd < 0.075. 

This shows that the situation is more complex and that one cannot reduce the analysis to a 
single varying parameter. Many studies then tried to determinate the sensitivity to the variation 
of many independent parameters. 

The sensitivity of the helium-4 abundance to the variation of 7 parameters was first investigated 
by Miiller et al. |361] considering the dependence on the parameters {Xi} = {G, a^.^, w, me, Tn, Qnp, 
B]j} independently, 

AlnyHc = ^c[^'AlnX, 
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Figure 6; (Left): variation of the light element abundances in function of rj compared to the 
spectroscopic abundances. The vertical line depicts the constraint obtained on rj from the study of 
the cosmic microwave background data. The lithium- 7 problem lies in the fact that rj^^^^^^^ < rj^^Ar- 
From Ref. fl02]. (right): Dependence of the light element abundance on the independent variation 
of the BBN parameters, assuming rj — rj^f^p^-p . From Ref. J 105^ 
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and assuming Aqcd fixed (so that the 7 parameters are in fact dimensionless quantities). The c^- 
are the sensitivities to the BBN parameters, assuming the six others fixed. It was concluded that 
Fho oc a-^ °43^2.4^o.024^o.24Q-i.8^o.53(^o.405 fo^. independent variations. They further related 
{Tn,Qnp, Bjj) to (aj.^ji'jTOejTONjTOd — Wu), as we shall discuss in the next section. 

This was generalized by Landau et al. [306] up to lithium-7 considering the parameters {a^M i j 
■'^QCD, ^bh^}, assuming G constant where the variation of Tn and the variation of the masses where 
tied to these parameters but the effect on the binding energies were not considered. 

Coc et al. |104) considered the effect of a variation of (Qnp, Bn^Tn, rric) on the abundances of the 
light elements up to lithium-7, neglecting the effect of a-^^^ on the cross-section. Their dependence 
on the independent variation of each of these parameters is depicted on Fig. [51 It confirmed the 
result of Refs. |206[I391] that the deuterium binding energy is the most sensitive parameter. From 
the helium-4 data alone, the bounds 

-8.2 X 10^2 < — < 6 X 10~^ -4 X 10^2 < < 2.7 X 10"^ (127) 

and 

- 7.5 X 10^2 < —-^^ < 6.5 X 10"^ (128) 
Bd 

at a 2(7 level, were set (assuming yy^^MAp)- The deuterium data set the tighter constraint —4 x 
10~^ < A In So < 3 X 10"^. Note also on Fig. |6] that the lithium-7 abundance can be brought in 
concordance with the spectroscopic observations provided that Bd was smaller during BBN 

-7.5 X 10"^ < — — ^ < -4 X 10"^ 
Bd 

so that Bd may be the most important parameter to resolve the lithium-7 problem. The effect of 
the quark mass on the binding energies was described in Ref. [45]. They then concluded that a 
variation of Amq/mq = 0.013 ±0.002 allows to reconcile the abundance of lithium-7 and the value 
of rj deduced from WMAP. 

This analysis was extended |145) to incorporate the effect of 13 independent BBN parameters 
including the parameters considered before plus the binding energies of deuterium, tritium, helium- 
3, helium-4, lithium-6, lithium-7 and beryllium-7. The sensitivity of the light element abundances 
to the independent variation of these parameters is summarized in Table I of Ref. |145] . These 
BBN parameters were then related to the same 6 "fundamental" parameters used in Ref. |361j . 

All these analysis demonstrate that the effects of the BBN parameters on the light element 
abundances are now under control. They have been implemented in BBN codes and most results 
agree, as well as with semi-analytical estimates. As long as these parameters are assume to vary 
independently, no constraints sharper than 10~^ can be set. One should also not forget to take 
into account standard parameters of the BBN computation such as 7] and the effective number of 
relativistic particle. 

3.8.3 Ftom BBN parameters to fundamental constants 

To reduce the number parameters, we need to relate the BBN parameters to more fundamental 
ones, keeping in mind that this can usually be done only in a model-dependent way. We shall 
describe some of the relations that have been used in many studies. They mainly concern Qnp, 7"n 
and Bd- 

At lowest order, all dimensional parameters of QCD, e.g. masses, nuclear energies etc., are to 
a good approximation simply proportional to some powers of Aqcd- One needs to go beyond such 
a description and takes the effects of the masses of the quarks into account. 
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Qnp can be expressed in terms of the mass on the quarks u and d and the fine-structure constant 

as 

Qnp = aa j^^^Aqcb + ("M - "lu), 

where the electromagnetic contribution today is (aaj,„AQCD)o = —0.76 MeV and therefore the 
quark mass contribution today is (toj — m^) = 2.05 |231j so that 



AQ„F 



Inp 



-0.59^^ 



1.59 



A(md 



(md - TOu) 



(129) 



All the analysis cited above agree on this dependence. 
The neutron lifetime can be well approximated by 



1 



1207r3 



Vq^ - l(2g^ - 9g2 - 8) + 15 In (q + - l) 



with q = Qnp/mc and Gp = l/^/2v'^. Using the former expression for Qnp we can express Tn in 
terms of a^.^ , v and the u, d and electron masses. It follows 



^ ^ 3.86^ 



1.52- 



Ar 



10.4 



A(md - TOu) 
(md - rriu) 



(130) 



Again, all the analysis cited above agree on this dependence. 

Let us now turn to the binding energies, and more particularly to Bjj that, as we have seen, 
is a crucial parameter. This is one the better known quantities in the nuclear domain and it is 
experimentally measured to a precision better than 10^^ 18 . Two approaches have been followed. 

• Pion mass. A first route is to use the dependence of the binding energy on the pion mass |1891 
[37], which is related to the u and d quark masses by 



mq(uu + dd) ^ ~ TO A 



QCD, 



where TOq = ^(tou + md) and assuming that the leading order of < 



dd) ^ depends only 



on Aqcd, Itt being the pion decay constant. This dependence was parameterized (549) as 



AB 



D 



At 



Bd 



where r is a fitting parameter found to be between 6 |189J and 10 [37] • Prior to this result, 
the analysis of Ref. |206] provides two computations of this dependence which respectively 
lead to r = —3 and r = 18 while, following the same lines, Ref. [57] got r — 0.082. 

Ref. |361) . following the computations of Ref. [420J . adds an electromagnetic contribution 
— 0.0081 Aq;j,j,/q;j,j, so that 



ABp 

Bd 



0.0081^. 



2 nin 



(131) 



but this latter contribution has not been included in other works. 

Sigma model. In the framework of the Walecka model, where the potential for the nuclear 
forces keeps only the a and uj meson exchanges. 



g2 g2 

V = --j-^ exp(-TOCTr) + exp(-TOtjr), 
47rr 47rr 
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where gs and gv are two coupling constants. Describing cr as a SU(3) singlet state, its mass 
was related to the mass of the strange quark. In this way one can hope to take into account 
the effect of the strange quark, both on the nucleon mass and the binding energy. In a second 
step Bo is related to the meson and nucleon mass by 

— - — = -48 1- 50 1- 6 

Bd rricr ttin 

so that AB]j /Bjj ~ — ITATOgmg |207| . Unfortunately, a complete treatment of all the nuclear 
quantities on mg has not been performed yet. 

The case of the binding energies of the other elements has been less studied. Ref. |145j follows a 
route similar than for B^ and relates them to pion mass and assumes that 

where fi are unknown coefficients assumed to be of order unity and Ai is the number of nucle- 
ons. No other estimates has been performed. Other nuclear potentials (such as Reid 93 potential, 
Nijmegen potential, Argonne wl8 potential and Bonn potential) have been used in Ref. |100] to 
determine the dependence of Bd on v and agree with previous studies. 

These analysis allow to reduce all the BBN parameter to the physical constants {a-^„ , w, mo, md— 
Wu, ?Tiq) and G that is not affected by this discussion. This set can be further reduce since all the 
masses can be expressed in terms of w as = /i^w, where hi are Yukawa couplings. 

To go further, one needs to make more assumption, such as grand unification, or by relating 
the Yukawa coupling of the top to v by assuming that weak scale is determined by dimensional 
transmutation jl04j . or that the varition of the constant is induced by a string dilaton [7S]. At 
each step, one gets more stringent constraints, which can reach the 10~^ |145J to lO^'"^ |104j level 
but indeed more model-dependent! 



3.8.4 Conclusion 

Primordial nucleosynthesis offers a possibility to test almost all fundamental constants of physics 
at a redshift of z ~ 10*. It is thus very rich but indeed the effect of each constant is more difficult 
to disentangle. The effect of the BBN parameters has been quantified with precision and they 
can be constrained typically at a 10~^ level, and in particular it seems that the most sensitive 
parameter is the deuterium binding energy. 

The link with more fundamental parameters is better understood but the dependence of the 
deuterium binding energy still left some uncertainties and a good description of the effect of the 
strange quark mass is missing. 

We have not considered the variation of G in this section. Its effect is disconnected from the 
other parameters. Let us just stress that assuming the BBN sensitivity on G by just modifying its 
value may be misleading. In particular G can vary a lot during the electron-positron annihilation so 
that the BBN constraints can in general not be described by an effective speed-up factor [1051 1129) . 



4 The gravitational constant 

The gravitational constant was the first constant whose constancy was questioned [154] . From 
a theoretical point of view, theories with a varying gravitational constant can be designed to 
satisfy the equivalence principle in its weak form but not in its strong form |536J . Most theories of 
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gravity that violate the strong equivalence principle predict that the locally measured gravitational 
constant may vary with time. 

The value of the gravitational constant is G = 6.67428(67) x 10~^^ m'^ • kg~^ • s~^ so that 
its relative standard uncertainty fixed by the 

codataEI 

in 2006 is 0.01%. Interestingly, the 
disparity between different experiments led, in 1998, to a temporary increase of this uncertainty 
to 0.15% |241j . which demonstrates the difficulty in measuring the value of this constant. This 
explains partly why the constraints on the time variation are less stringent than for the other 
constants. 

A variation of the gravitational constant, being a pure gravitational phenomenon, does not 
affect the local physics, such as e.g. the atomic transitions or the nuclear physics. In particular, it 
is equivalent at stating that the masses of all particles are varying in the same way to that their 
ratios remain constant. Similarly all absorption lines will be shifted in the same way. It follows 
that most constraints are obtained from systems in which gravity is non-negligible, such as the 
motion of the bodies of the Solar system, astrophysical and cosmological systems. They are mostly 
related in the comparison of a gravitational time scale, e.g. period of orbits, to a non-gravitational 
time scale. It follows that in general the constraints assume that the values of the other constants 
are fixed. Taking their variation into account would add degeneracies and make the constraints 
cited below less stringent. 

We refer to § IV of FVC ^495. for earlier constraints based e.g. on the determination of the 
Earth surface temperature, which roughly scales as G^'^^A/q''^ and gives a constraint of the order 
of |AG/G| < 0.1 [224], or on the estimation of the Earth radius at different geological epochs. 



4.1 Solar systems constraints 

Monitoring the orbits of the various bodies of the Solar system offers a possibility to constrain 
deviations from general relativity, and in particular the time variation of G. This accounts for 
comparing a gravitational time scale (related to the orbital motion) and an atomic time scale and 
it is thus assumed that the variation of atomic constants is negligible on the time of the experiment. 

The first constraint arises from the Earth-Moon system. A time variation of G is then related 
to a variation of the mean motion (n = 27r/P) of the orbit of the Moon around the Earth. A 
decrease in G would induce both the Lunar mean distance and period to increase. As long as the 
gravitational binding energy is negligible, one has 

Earlier constraints rely on paleontological data and ancient eclipses obervations (see § IV.B.l of FVC |495J ) 
and none of them are very reliable. A main difficulty arises from tidal dissipation that also causes 
the mean distance and orbital period to increase (for tidal changes 2h/n + ia/ a = 0), but not as 
in the same ratio as for G. 

The Lunar Laser Ranging (LLR) experiment has measured the relative position of the Moon 
with respect to the Earth with an accuracy of the order of 1 cm over 3 decades. An early analysis 
of this data [538j assuming a Brans-Dicke theory of gravitation gave that |G/G| < 3 x 10^^^ yr^^. 
It was improved |362' by using 20 years of observation to get |G/G| < 1.04 x 10~^^yr~^, the 
main uncertainty arising from Lunar tidal acceleration. With, 24 years of data, one reached [539] 
|G/G| < 6 X 10~^^yr~^ and finally, the latest analysis of the Lunar laser ranging experiment [540] 
increased the constraint to 

I =(4±9)x 10-"yr-i. (133) 





The CODATA is the COmmittee on Data for Science and Technology, see http://www.codata.org/. 
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Similarly, Shapiro et al. [454] compared radar-echo time delays between Earth, Venus and 
Mercury with a caesium atomic clock between 1964 and 1969. The data were fitted to the the- 
oretical equation of motion for the bodies in a Schwarzschild spacetime, taking into account the 
perturbations from the Moon and other planets. They concluded that \G/G\ < 4 x 10" 



-10 



yr" 



The data concerning Venus cannot be used due to imprecision in the determination of the por- 
tion of the planet reflecting the radar. This was improved to \G/G\ < 1.5 x 10~^'^yi~^ by in- 
cluding Mariner 9 and Mars orbiter data [422]. The analysis was further extended |453j to give 



G/G = (-2± 10) X 10- 
data gives [12] 



•yr 



The combination of Mariner 10 an Mercury and Venus ranging 



(0.0 ±2.0) X lO-'Vr 



(134) 



Reasenberg et al. [423] considered the 14 months data obtained from the ranging of the Viking 
spacecraft and deduced, assuming a Brans-Dicke theory, \G/G\ < 10^^^ yr^^. HeUings et al. |248| 
using all available astrometric data and in particular the ranging data from Viking landers on Mars 
deduced that 

G 



G 



(2 ±4) X 10 



-12 



yr 



(135) 



The major contribution to the uncertainty is due to the modeling of the dynamics of the aster- 
oids on the Earth-Mars range. Hellings et al. [248] also tried to attribute their result to a time 
variation of the atomic constants. Using the same data but a different modeling of the asteroids, 
Reasenberg [424j got \G/G\ < 3 x lO^^^yr"^, which was then improved by Chandler et al. [92] to 
|G/G| < 10-"yr-i. 



4.2 Pulsar timing 

Contrary to the Solar system case, the dependence of the gravitational binding energy cannot be 
neglected while computing the time variation of the period. Here two approaches can be followed; 
either one sticks to a model (e.g. scalar-tensor gravity) and compute all the effects in this model 
or one has a more phenomenological approach and tries to put some model-independent bounds. 

Eardley |176j followed the first route and discussed the effects of a time variation of the gravi- 
tational constant on binary pulsar in the framework of the Brans-Dicke theory. In that case, both 
a dipole gravitational radiation and the variation of G induce a periodic variation in the pulse 
period. Nordtvedt |382] showed that the orbital period changes as 

P 
P 

where Ci = Slnmi/SlnG. He concluded that for the pulsar PSR 1913-fl6 (mi ~ m2 and ci ~ C2) 
one gets 

p = -[2 + 5c]g, (137) 

the coefficient c being model dependent. As another application, he estimated that c-^^^^-^ ^ — 5 x 
10~^°, c„„„„ ~ — 10~* and Cg^^ ~ —4 x 10~^ justifying the formula used in the Solar system. 

Damour et al. |133j used the timing data of the binary pulsar PSR 1913-1-16. They imple- 
mented the effect of the time variation of G by considering the effect on P/P. They defined, 
in a phenomenological way, that G/G — —0.5SP/P, where SP is the part of the orbital pe- 
riod derivative that is not explained otherwise (by gravitational waves radiation damping). This 



2(r7iici -f- m2C2) + 3(miC2 + m2Ci) G 

TOl -|- 7712 G 



(136) 



72 



theory-independent definition has to be contrasted with the theory-dependent resuh ()137p by 
Nordtvedt [352] ■ They got 

G/G= (1.0 ±2.3) X 10"" yr-\ (138) 

Damour and Taylor [132] then reexamined the data of PSR 1913-1-16 and established the upper 
bound 

G/G < (1.10 ± 1.07) X 10"" yr-\ (139) 
Kaspi et al. used data from PSR B1913+16 and PSR B1855-f09 respectively to get 



and 



G/G= (4±5) X 10"i2yj,-i 



G/G = (-9 ±18) X IQ-^'^yv-^, 



(140) 
(141) 



the latter case being more "secure" since the orbiting companion is not a neutron star. 

All the previous results concern binary pulsars but isolated ones can also be used. Heintzmann 
and Hillebrandt [247 related the spin-down of the pulsar JP1953 to a time variation of G. The 
spin-down is a combined effect of electromagnetic losses, emission of gravitational waves, possible 
spin-up due to matter accretion. Assuming that the angular momentum is conserved so that 
I/P =constant, one deduces that 



din/ 
dlnG 



G 
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The observational spin-down can be decomposed as 
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(143) 



Since P/P^^^ and P/P^^ are positive definite, it follows that P/ P^^^ > P/Pg so that a bound on G 
can be inferred if the main pulse period is the period of rotation. Heintzmann and Hillebrandt |247| 
then modelled the pulsar by a polytropic (P cx p") white dwarf and deduced that din //dlnG = 
2 — 3n/2 so that |G/G| < 10~^"yr~^. Mansfield |342j assumed a relativistic degenerate, zero 
temperature polytropic star and got that, when G < 0, < —G/G < 6.8 x 10~-^^yr~^ at a 2a 
level. He also noted that a positive G induces a spin-up counteracting the electromagnetic spin- 
down which can provide another bound if an independent estimate of the pulsar magnetic field can 
be obtained. Goldman |234] . following Eardley [176] . used the scaling relations N oc G~^/^ and 
M oc G"'^/^ to deduce that 2d In //dlnG = -5 + 3dln//dlniV. He used the data from the pulsar 
PSR 0655+64 to deduce that the rate of decrease of G was smaller than 



< -G/G < 5.5 X 10" 



(144) 



The analysis [512] of 10 years high pecision timing data on the millisecond pulsar PSR J0437-4715 
has allowed to improve the constraint to 



|G/G| < 2.3 X 10" 



yr 



(145) 



Recently, it was argued [2651 1426) that a variation of G would induce a departure of the neutron 
star matter from /3-equilibrium, due to the changing hydrostatic equilibrium. This would force non- 
equilibrium /3-processes to occur, which release energy that is invested partly in neutrino emission 
and partly in heating the stellar interior. Eventually, the star arrives at a stationary state in which 
the temperature remains nearly constant, as the forcing through the change of G is balanced 
by the ongoing reactions. Comparing the surface temperature of the nearest millisecond pulsar. 
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PSR J0437-4715, inferred from ultraviolet observations, two upper limits for this variation were 
obtained, \GIG\ < 2 x 10^^" yr^^, direct Urea reactions operating in the neutron star core are 
allowed, and \G/G\ < 4 x 10~^^ yr"^, considering only modified Urea reactions. This was extended 
in Ref. [300] in order to take into account the correlation between the surface temperatures and 
the radii of some old neutron stars to get \G/G\ < 2.1 x 10^^^ yr^^. 

4.3 Stellar constraints 

Early works, see § IV. C of FVC [495', studied the Solar evolution in presence of a time varying 
gravitational constant, concluding that under the Dirac hypothesis, the original nuclear resources 
of the Sun would have been burned by now. This results from the fact that an increase of the 
gravitational constant is equivalent to an increase of the star density (because of the Poisson 
equation). 

The idea of using stellar evolution to constrain the possible value of G was originally proposed by 
Teller [481] . who stressed that the evolution of a star was strongly dependent on G. The luminosity 
of a main sequence star can be expressed as a function of Newtons gravitational constant and its 
mass by using homology relations |224[ 1481] . In the particular case that the opacity is dominated 
by free- free transitions, Gamow [224' found that the luminosity of the star is given approximately 
by L oc G'^^^M^■^. In the case of the Sun, this would mean that for higher values of G, the burning 
of hydrogen will be more efficient and the star evolves more rapidly, therefore we need to increase 
the initial content of hydrogen to obtain the present observed Sun. In a numerical test of the 
previous expression, Delg'Innocenti et al. [139] found that low-mass stars evolving from the Zero 
Age Main Sequence to the red giant branch satisfy L cx G^-^Af^-^, which agrees to within 10% of 
the numerical results, following the idea that Thomson scattering contributes significantly to the 
opacity inside such stars. Indeed, in the case of the opacity being dominated by pure Thomson 
scattering, the luminosity of the star is given by L cx G*M^. It follows from the previous analysis 
that the evolution of the star on the main sequence is highly sensitive to the value of G. 

The driving idea behind the stellar constraints is that a secular variation of G leads to a variation 
of the gravitational interaction. This would affect the hydrostatic equilibrium of the star and in 
particular its pressure profile. In the case of non-degenerate stars, the temperature, being the only 
control parameter, will adjust to compensate the modification of the intensity of the gravity. It 
will then affect the nuclear reaction rates, which are very sensitive to the temperature, and thus 
the nuclear time scales associated to the various processes. It follows that the main stage of the 
stellar evolution, and in particular the lifetimes of the various stars, will be modified. As we shall 
see, basically two types of methods have been used, the first in which on relate the variation of G 
to some physical characteristic of a star (luminosity, effective temperature, radius), and a second 
in which only a statistical measurement of the change of G can be infered. Indeed, the first class 
of methods are more reliable and robust but is usually restricted to nearby stars. Note also that 
they usually require to have a precise distance determination of the star, which may depend on G. 

4.3.1 Ages of globular clusters 

The first application of these idea has been performed with globular clusters. Their ages, deter- 
mined for instance from the luminosity of the main-sequence turn-off, have to be compatible with 
the estimation of the age of the Galaxy. This gives the constraint ^139] 

G/G= (-1.4 ±2.1) X 10~"yr"^ (146) 

The effect of a possible time dependence of G on luminosity has been studied in the case of 
globular cluster H-R diagrams but has not yielded any stronger constraints than those relying on 
celestial mechanics 
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4.3.2 Solar and stellar sysmology 



A side effect of the change of luminosity is a change in the depth of the convection zone so that the 
inner edge of the convecting zone changes its location. This induces a modification of the vibration 
modes of the star and particularly to the acoustic waves, i.e p- modes [140] . 

Helioseismology. This waves are observed for our star, the Sun, and heliosysmology allows 
to determine the sound speed in the core of the Sun and, together with an equation of state, 
the central densities and abundances of helium and hydrogen. Demarque et al. [140] considered 
an ansatz in which G cx and showed that < 0.1 over the last 4.5 x 10^ years, which 
corresponds to \G/G\ < 2 x 10~^-^yr~^. Guenther et al. |239j also showed that g-modes could 
provide even much tighter constraints but these modes are up to now very difficult to observe. 
Nevertheless, they concluded, using the claim of detection by HiU and Gu [250] , that \G/G\ < 
4.5 X 10~^^yr~^. Guenther et al. |240J then compared the p-mode spectra predicted by different 
theories with varying gravitational constant to the observed spectrum obtained by a network of 
six telescopes and deduced that 

G/G < 1.6 X 10^12 y^-i^ (^47) 

The standard Solar model depends on few parameters and G plays a important role since stellar 
evolution is dictated by the balance between gravitation and other interactions. Astronomical ob- 
servations determines GMq with an accuracy better than lO"'^ and a variation of G with GMq fixed 
induces a change of the pressue (P = GAIq/Rq) and density {p — Mq/Rq). The experimental 
uncertainties in G between different experiments have important implications for helioseismology. 
In particulat the uncertainties for the standard solar model lead to a range in the value of the 
sound speed in the nuclear region that is as much as 0.15% higher than the inverted helioseismic 
sound speed [333] . While a lower value of G is preferred for the standard model, any definite 
prediction is masked by the uncertainties in the solar models available in the literature. Ricci and 
Villante |430j studied the effect of a variation of G on the density and pressure profile of the Sun 
and concluded that present data cannot constrain G better than 10~^%. It was also shown |333j 
that the information provided by the neutrino experiments is quite significant because it consti- 
tutes an independent test of G complementary to the one provided by helioseismology. 



White dwarfs. The observation of the period of non-radial pulsations of white dwarf allows 
to set similar constraints. White dwarfs represent the final stage of the stellar evolution for stars 
with a mass smaller to about IOMq. Their structure is supported against gravitational collapse by 
the pressure of degenerate electrons. It was discovered that some white dwarfs are variable stars 
and in fact non-radial pulsator. This opens the way to use seismological techniques to investigate 
their internal propoerties. In particular, their non-radial oscillations is mostly determined by the 
Brunt- Vaisala frequency 

dlnPi/71/p 

^ = .9 ^ 

ar 

where g is the gravitational acceleration, Fi the first adiabatic exponent and P and p the pressure 
and density (see e.g. Ref. |283] for a white dwarf model taking into account a varying G). A 
variation of G induces a modification of the degree of degeneracy of the white dwarf, hence on the 
frequency N as well as the cooling rate of the star, even though this is thought to be negligible 
at the luminosities where white dwarfs are pulsationally unstable |52j. Using the observation of 
G117-B15A that has been monitored during 20 years, it was concluded [42] that 

- 2.5 X 10"^" yr^i < G/G < 4.0 x 10"" yr-\ (148) 
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at a 2tT-level. The same observations were reanalyzed in Ref. |52 to obtain 



|G/G| < 4.1 X lO^^yr-^ (149) 
4.3.3 Late stages of stellar evolution and supernovae 

A variation of G can influence the white dwarf coohng and the hght curves ot Type la supernovae. 

Garcia-Berro et al. [227\ considered the effect of a variation of the gravitational constant on 
the cooling of white dwarfs and on their luminosity function. As first pointed out by Vila |514] , 
the energy of white dwarfs, when they are cool enough, is entirely of gravitational and thermal 
origin so that a variation of G will induce a modification of their energy balance and thus of their 
luminosity. Restricting to cold white dwarfs with luminosity smaller than ten Solar luminosity, the 
luminosity can be related to the star binding energy B and gravitational energy, E , as 

which simply results from the hydrostatic equilibrium. Again, the variation of the gravitational 
constant intervenes via the Poisson equation and the gravitational potential. The cooling process 
is accelerated if G/G < which then induces a shift in the position of the cut-off in the luminosity 
function. Garcia-Berro et al. [22T concluded that 

0< -G/G < (1±1) X 10"" yr"^ (151) 

The result depends on the details of the cooling theory, on whether the C/0 white dwarf is strat- 
ified or not and on hypothesis on the age of the galactic disk. For instance, with no stratification 
of the C/0 binary mixture, one would require G/G = —(2.5 ± 0.5) x 10~^^ yr~^ if the Solar neigh- 
borhood has a value of 8 Gyr (i.e. one would require a variation of G to explain the data). In the 
case of the standard hypothesis of an age of 11 Gyr, one obtains that < —G/G < 3 x 10^^^ yr~^- 

The late stages of stellar evolution are governed by the Chandrasekhar mass {hc/GY^'^m~'^ 
mainly determined by the balance between the Fermi pressure of a degenerate electron gas and 
gravity. 

Simple analytical models of the light curves of Type la supernovae predict that the peak of 
luminosity is proportional to the mass of nickel synthetized. In a good approximation, it is a 
fixed fraction of the Chandrasekhar mass. In models allowing for a varying G, this would induce a 
modification of the luminosity distance- redshift relation |226l 12331 1429] . It was however shown that 
this effect is small. Note that it will be degenerate with the cosmological parameters. In particular, 
the Hubble diagram is sensitive to the whole history of G{t) between the highest redshift observed 
and today so that one needs to rely on a better defined model, such as e.g. scalar-tensor theory |429j 
(the effect of the Fermi constant was also considered in Ref. }195j ). 

In the case of Type II supernovae, the Chandrasekhar mass also gouvernes the late evolution- 
ary stages of massive stars, including the formation of neutron stars. Assuming that the mean 
neutron star mass is given by the Chandrasekhar mass, one expects that G/G = —2M^^/iM^^. 
Thorsett |486| used the observations of five neutron star binaries for which five Keplerian parame- 
ters can be determined (the binary period P;,, the projection of the orbital semi-major axis ai sini, 
the eccentricity e, the time and longitude of the periastron To and w) as well as the relativistic 
advance of the angle of the periastron lu. Assuming that the neutron star masses vary slowly as 
M„g = M'^^ — Mj^^t^g, that their age was determined by the rate at which Pt, is increasing (so 
that tpfs — SPfe/A) and that the mass follows a normal distribution, Thorsett 486 deduced that, 
at 2a, 

G/G = (-0.6 ±4.2) X 10^12 yj.-i^ (152) 
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4.3.4 New developments 

It has recently been proposed that the variation of G inducing a modification of binary's binding 
energy, it should affect the gravitational wave luminosity, hence leading to corrections in the 
chirping frequency [550) . For instance, it was estimated that a LISA observation of an equal-mass 
inspiral event with total redshifted mass of lO^M© for three years should be able to measure 
G/G at the time of merger to better than 10~^^/yr. This method paves the way to constructing 
constraints in a large band of redshifts as well as in different directions in the sky, which would be 
an unvaluable constraint for many models. 

More speculative is the idea [23] that a variation of G can lead a neutron to enter into the region 
where strange or hydrid stars are the true ground state. This would be associated to Gamma- 
Ray-Burst that are claimed to be able to reach the level of 10^ ^^/yr on the time variation of 
G. 

4.4 Cosmological constraints 

Cosmological observations are more difficult to use in order to set constraints on the time variation 
of G. In particular, they require to have some ideas about the whole history of G as a function 
of time but also, as the variation of G reflects an extension of General relativity, it requires to 
modify all equations describing the evolution (of the universe and of the large scale structure) in 
a consistent way. We refer to Refs. ^499. .4971 1501J for a discussion of the use of cosmological data 
to constrain deviations from general relativity. 

4.4.1 Cosmic microwave background 

A time-dependent gravitational constant will have mainly three effects on the CMB angular power 
spectrum (see Ref. |429j for discussions in the framework of scalar-tensor gravity in which G is 
considered as a field): 

1. The variation of G modifies the Friedmann equation and therefore the age of the Universe 
(and, hence, the sound horizon). For instance, if G is larger at earlier time, the age of the 
Universe is smaller at recombination, so that the peak structure is shifted towards higher 
angular scales. 

2. The amplitude of the Silk damping is modified. At small scales, viscosity and heat conduction 
in the photon-baryon fluid produce a damping of the photon perturbations. The damping 
scale is determined by the photon diffusion length at recombination, and therefore depends 
on the size of the horizon at this epoch, and hence, depends on any variation of the Newton 
constant throughout the history of the Universe. 

3. The thickness of the last scattering surface is modified. In the same vein, the duration of 
recombination is modified by a variation of the Newton constant as the expansion rate is 
different. It is well known that CMB anisotropics are affected on small scales because the 
last scattering "surface" has a finite thickness. The net effect is to introduce an extra, roughly 
exponential, damping term, with the cutoff length being determined by the thickness of the 
last scattering surface. When translating redshift into time (or length), one has to use the 
Friedmann equations, which are affected by a variation of the Newton constant. The relevant 
quantity to consider is the visibility function g. In the limit of an infinitely thin last scattering 
surface, r goes from cx) to at recombination epoch. For standard cosmology, it drops from 
a large value to a much smaller one, and hence, the visibility function still exhibits a peak, 
but it is much broader. 
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In full generality, the variation of G on the CMB temperature anisotropies depends on many 
factors: (1) modification of the background equations and the evolution of the universe, (2) modi- 
fication of the perturbation equations, (3) whether the scalar field inducing the time variaiton of 
G is negligible or not compared to the other matter components, (4) on the time profile of G that 
has to be determine to be consistent with the other equations of evolution. This explains why it 
is very difficult to state a definitive constraint. For instance, in the case of scalar-tensor theories 
(see below), one has two arbitrary functions that dictate the variation of G. As can be seen e.g. 
from Ref. [4291 1375] , the profiles and effects on the CMB can be very different and difficult to com- 
pare. Indeed, the effects described above are also degenerate with a variation of the cosmological 
parameters. 

In the case of Brans-Dicke theory, one just has a single constant parameter uj-^-^ characterizing 
the deviation from general relativity and the time variation of G. It is thus easier to compare the 
different constraints. Chen and Kamionkowski [S3] showed that CMB experiments such as WMAP 
will be able to constrain these theories for ijj^^ < 100 if all parameters are to be determined by 
the same CMB experiment, < 500 if all parameters are fixed but the CMB normalization and 
< 800 if one uses the polarization. For the Planck mission these numbers are respectively, 
800, 2500 and 3200. Ref. 2^ concluded from the analysis of WMAP, ACBAR, VSA and CBI, and 
galaxy power spectrum data from 2dF, that > 120, in agreement with the former analysis 
of Ref. [575 1. An analysis (546j indictates that The WMAP-5yr data and the all CMB data both 
favor a slightly non-zero (positive) G/G but with the addition of the SDSS poser spectrum data, 
the best-fit value is back to zero, concluding that —0.083 < AG/G < 0.095 between recombination 
and today, which corresponds to —1.75 x 10^^^ yr^^ < G/G < 1.05 x 10~^^yr~-^. 

From a more phenomenoloical prospect, some works modelled the variation of G with time in 
a purely ad-hoc way, for instance (88] by assuming a linear evolution with time or a step function. 

4.4.2 BBN 

As explained in details in section [3.8.11 changing the value of the gravitational constant affects the 
freeze-out temperature Tf . A larger value of G corresponds to a higher expansion rate. This rate is 
determined by the combination Gp and in the standard case the Friedmann equations imply that 
Gpt^ is constant. The density p is determined by the number N^, of relativistic particles at the 
time of nucleosynthesis so that nucleosynthesis allows to put a bound on the number of neutrinos 
N^. Equivalently, assuming the number of neutrinos to be three, leads to the conclusion that G 
has not varied from more than 20% since nucleosynthesis. But, allowing for a change both in G 
and Ni, allows for a wider range of variation. Contrary to the fine structure constant the role of G 
is less involved. 

The effect of a varying G can be described, in its most simple but still useful form, by introducing 
a speed-up factor, ^ = H/Hgr, that arises from the modification of the value of the gravitational 
constant during BBN. Other approaches considered the full dynamics of the problem but restricted 
themselves to the particular class of Jordan-Fierz-Brans-Dicke theory [U HSj [24l [83] IIOH I125| 14351 
1547) (Casas et al. [53] concluded from the study of helium and deuterium that oj^^^ > 380 when N^, = 
3 and uj^^ > 50 when = 2.), of a massless dilaton with a quadratic coupling [1051 11061 11291 1440] 
or to a general massless dilaton |447) . It should be noted that a combined analysis of BBN and 
CMB data was investigated in Refs. [1121 1291) . The former considered G constant during BBN 
while the latter focused on a nonminimally quadratic coupling and a runaway potential. It was 
concluded that from the BBN in conjunction with WMAP determination of 77 set that AG/G has 
to be smaller than 20%. We however stress that the dynamics of the field can modify CMB results 
(see previous section) so that one needs to be careful while inferring Ob from WMAP unless the 
scalar-tensor theory has converged close to general realtivity at the time of decoupling. 

In early studies, Barrow [24] assumed that G oc and obtained from the helium abundances 
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that -5.9x10^3 <n< 7 x 10"^ which implies that \G/G\ < (2±9.3) /ix lO^^^ yr ^ , assuming a flat 
universe. This corresponds in terms of the Brans-Dicke parameter to oj-^^^ > 25. Yang et al. |547| 
included the deuterium and lithium to improve the constraint to n < 5 x 10~^ which corresponds 
to > 50. It was further improved by Rothman and Matzner [435] to |n| < 3 x 10~^ implying 
\G/G\ < 1.7 X 10^^'^ yr^^. Accetta et al. p] studied the dependence of the abundances of D, ^He, 
^He and ""Li upon the variation of G and concluded that —0.3 < AG/G < 0.4 which roughly 
corresponds to \G/G\ < 9 x lO"""^^ yr""'^. All these investigations assumed that the other constants 
are kept fixed and that physics is unchanged. Kolb et al. [294] assumed a correlated variation of 
G, a^-^^ and Gp and got a bound on the variation of the radius of the extra-dimensions. 

Although the uncertainty in the helium-4 abundance has been argued to be significantly larger 
that what was assumed in the past [394] , interesting bounds can still be derived [116] . In particular 
translating the bound on extra relativistic degress of freedom (—0.6 < SN^, < 0.82) to a constraint 
on the speed-up factor (0.949 < ^ < 1.062), it was concluded [TTB] , since AG/G = _ ^ = 
7(5iV^/43, that 

AC 

- 0.10 < — < 0.13. (153) 

The relation between the speed-up factor, or an extra number of relativistic degrees of freedom, 
with a variation of G is only approximate since it assumes that the variation of G affects only the 
Friedmann equation by a renormalisation of G. This is indeed accurate only when the scalar field 
is slow-rolling. For instance |105] . the speed-up factor is given (with the notations of Section [5. l.ip 

by 



so that 



G (l + a(^,)^'j2 



= — ^ ' "v^"^^*^ (154) 

^ Go(l + «2)(l_^2,/3)' 

so that AG /Go = — 1 only if a ^ 1 (small deviation from general relativity) and ip'^ <C 1 (slow 
rolling dilaton) . The BBN in scalar-tensor theories was investigated [1051 1129] in the case of a two- 
parameter family involving a non-linear scalar field-matter coupling function. They concluded that 
even in the cases where before BBN the scalar-tensor theory was far from general relativity, BBN 
enables to set quite tight constraints on the observable deviations from general relativity today. In 
particular, neglecting the cosmological constant, BBN imposes ag < 10~^'^ I3~^ {Q^a_th'^ /O-^^)^^^^ 
when /? > 0.5 (with the definitions introduced below Eq. ([163P '). 



5 Theories with varying constants 

As explained in the introduction, Dirac postulated that G varies as the inverse of the cosmic time. 
Such an hypothesis is indeed not a theory since the evolution of G with time is postulated and not 
derived from an equation of evolutiorF^ consistent with the other field equations, that have to take 
into account that G is no more a constant (in particular in a Lagrangian formulation one needs to 
take into account that G is no more constant when varying. 

The first implementation of Dirac's phenomenological idea into a field-theory framework (i.e. 
modifying Einstein gravity and incorporating non-gravitational forces and matter) was proposed 

^"Note that Dirac hypothesis can also be achieved by assuming that e varies as f^/^. Indeed this reflects a choice 
of units, either atomic or Planck units. There is however a difi'erence: assuming that only G varies violates the 
strong equivalence principle while assuming a varying e results in a theory violating the weak equivalence principle. 
It does not mean we are detecting the variation of a dimensionful constant but simply that either /he or Gm^/hc 
is varying. This shows that many implementation of this idea are a priori possible. 
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by Jordan |267j . He realized that the constants have to become dynamical fields and proposed the 
action 



(155) 



77 and ^ being two parameters. It follows that both G and the fine-structure constant have been 
promoted to the status of a dynamical field. 

Fierz |196j realized that with such a Lagrangian, atomic spectra will be space-time-dependent, 
and he proposed to fix ij to the value -1 to prevent such a space-time dependence. This led to the 
definition of a one- parameter (^) class of scalar-tensor theories in which only G is assumed to be 
a dynamical field. This was then further explored by Brans and Dicke |66) (with the change of 
notation £^ uj). In this Jordan- Fierz-Brans-Dicke theory the gravitational constant is replaced 
by a scalar field which can vary both in space and time. It follows that, for cosmological solutions, 
G oc t~"' with = 2 + 3uj-g^/2. Einstein gravity is thus recovered when — > 00. This kind of 
theory was further generalized to obtain various functional dependencies for G in the formalisation 
of scalar-tensor theories of gravitation (see e.g. Damour and Esposito-Farese |123| or Will |536] ). 

5.1 Introducing new fields: generalities 
5.1.1 The example of scalar-tensor theories 

Let us start to remind how the standard general relativistic framework can be extended to make 
G dynamical on the example of scalar-tensor theories, in which gravity is mediated not only by a 
massless spin-2 graviton but also by a spin-0 scalar field that couples universally to matter fields 
(this ensures the universality of free fall). In the Jordan frame, the action of the theory takes the 
form 

/d^ 
[F(^)i? - g^'^Z(^)¥.,^^,, - 2Uip)] + 5„,attcr[V'; 9^.] (156) 

where is the bare gravitational constant. This action involves three arbitrary functions (F, Z 
and U) but only two are physical since there is still the possibility to redefine the scalar field. F 
needs to be positive to ensure that the graviton carries positive energy, i^inatter is the action of 
the matter fields that are coupled minimally to the metric (7^,^. In the Jordan frame, the matter is 
universally coupled to the metric so that the length and time as measured by laboratory apparatus 
are defined in this frame. 

The variation of this action gives the following field equations 



F(^) ( i?^, - Ifj^.R ] = SttG.T^, + Z{^) 



+V^d,F{ip) - g^.nF(^) - g^^Uiifi) , (157) 

2Z(^)n^ = -^i?-^(9„^)2 + 2^ , (158) 
dtp dip dip 

V^Ti; = 0, (159) 

where T = is the trace of the matter energy-momentum tensor T^" = {2/^—g) x SSm/Sg^i:^- 
As expected |184) . we have an equation which reduces to the standard Einstein equation when ip 
is constant and a new equation to describe the dynamics of the new degree of freedom while the 
conservation equation of the matter fields is unchanged, as expected from the weak equivalence 
principle. 
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It is useful to define an Einstein frame action through a conformal transformation of the metric 
= F{Lp)g^^. In the fohowing all quantities labelled by a star (*) will refer to Einstein frame. 
Defining the field 1^9* and the two functions A((y9,) and l^(iy9») (see e.g. Ref. [192] ) by 



dv^^y 3 ^ dlnf(v.) V , 1 



(VJ), 2\/(^,) = f/((p)F-2(^), 



the action (|156p reads as 



1 



S = -r^ / d*a;yq^[i?, - 25r5^<^*a,(p, - 4^] + 5,„attcr [A^ff* ; V-] • (160) 



167rG* 

The kinetic terms have been diagonalized so that the spin-2 and spin-0 degrees of freedom of the 
theory are perturbations of g*^^ and Lp^, respectively. In this frame the field equations are given by 

R%-\r*9% = 87rG,T;, + 2d,^,d,^,~g%{gfd^^,dp^,)-2V{v)g;,, (161) 

□*<y9* = —4^TG^,a{^p^,)T^,+dV{^p)/d(p^. , (162) 
V;r,^, = a(vp,) T,d,^, , (163) 

with a = dhiA/dip^ and /3 = da/dip^. In this version, the Einstein equations are not modified, 
but since the theory can now be seen as the theory in which all the mass are varying in the same 
way, there is a source term to the conservation equation. This shows that the same theory can 
be interpreted as a varying G theory or a universally varying mass theory, but remember that 
whathever its form the important parameter is the dimensionless quantity Gm? /he. 

The action ()156p defines an effective gravitational constant G^a — Gf,/F — Gf,A^ . This constant 
does not correspond to the gravitational constant effectively measured in a Cavendish experiment. 
The Newton constant measured in this experiment is 

Gcav = G.Al{l + al) = ^l^+ 2F+{fi ) ^^^^^ 

where the first term, G^,Aq corresponds to the exchange of a graviton while the second term G^AqUq 
is related to the long range scalar force. The gravitational constant depends on the scalar field and 
is thus dynamical. 

This illustrates the main features that will appear in any such models: (i) new dynamical fields 
appear (here a scalar field), (ii) some constant will depend the value of this scalar field (here G 
is a function of the scalar field). It follows that the Einstein equations will be modified and that 
there will exist a new equation dictating the propagation of the new degree of freedom. 

In this particular example, the coupling of the scalar field is universal so that no violation of 
the universality of free fall is expected. The deviation from general relativity can be quantified in 
terms of the post-Newtonian parameters, which can be expressed in terms of the values of a and 
/3 today as 

' 1 + ag' ^ 2(l + ag)2 ^ ' 

This expression are valid only if the field is light on the Solar system scales. If this is not the case 
then this conclusions may be changed |287j . The Solar system constraints imply ao to be very 
small, typically ofg < 10~^ while /3q can still be large. Binary pulsar observations [12411190] impose 
that /3o > —4.5. The time variation of G is then related to ao, /3o and the time variation of the 
scalar field today 

^ = 2ao (1 + (166) 
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This example shows that the variation of the constant and the deviation from general relativity 
quantified in terms of the PPN parameters are of the same magnitude, because they are all driven 
by the same new scalar field. 

The example of scalar-tensor theories is also very illustrative to show how deviation from general 
relativity can be fairly large in the early universe while still being compatible with Solar system 
constraints. It relies on the attraction mechanism toward general relativity [lTl\ 1128] . 

Consider the simplest model of a massless dilaton — 0) with quadratic coupling (In A = 

a = ^/3ip1). Note that the linear case correspond to a Brans-Dicke theory with a fixed deviation 
from general relativity. It follows that ao = f3ipo^, and /3o = /3. As long a.sV = 0, the Klein-Gordon 
equation can be rewritten in terms of the variable p = In a as 

-^ip': + (1 - w)^'^ = -a(^.)(l - 3^. (167) 

As emphasized in Ref. [127J . this is the equation of motion of a point particle with a velocity 
dependent inertial mass, m(y*) = 2/(3 — (f^) evolving in a potential Q!((^»)(1 — 3w) and subject 
to a damping force, —(1 — w)ip'^. During the cosmological evolution the field is driven toward the 
minimum of the coupling function. If /3 > 0, it drives toward 0, that is a — >■ 0, so that the 
scalar-tensor theory becomes closer and closer to general relativity. When f3 < 0, the theory is 
driven way from general relativity and is likely to be incompatible with local tests unless ip^, was 
initially arbitrarily close from 0. 

It follows that the deviation from general relativity remains constant during the radiation era 
(up to threshold effects in the early universe [10711129] and quantum effects [84]) and the theory is 
then attracted toward general relativity during the matter era. Note that it implies that postulating 
a linear or inverse variation of G with cosmic time is actually not realistic in this class of models. 
Since the theory is fully defined, one can easily compute various cosmological observables (late time 
dynamics [515], CMB anisotropy [HI], weak lensing [Hg, BBN [TUS] UMl [HH]) in a consistent 
way and confront them with data. 



5.1.2 Making other constants dynamical 

Given this example, it seems a priori simple to cook up a theory that will describe a varying 
fine-structure constant by coupling a scalar field to the electromagnetic Faraday tensor as 



S = 



R 



16nG 



- 2(9^0)2 - ii3(0)F2, 



- gd X 



(168) 



so that the fine-structure will evolve according to a = B^^. Such an simple implementation may 
however have dramatic implications. In particular, the contribution of the electromagnetic binding 
energy to the mass of any nucleus can be estimated by the Bethe-Weizacker formula so that 

miA,z) (0) 3 98.25 a(</)) ^^^^JT^ MeV. 

This implies that the sensitivity of the mass to a variation of the scalar field is expected to be of 
the order of 



f{A.z) = d^m^A.z){4>) 



(169) 



It follows that the level of the violation of the universality of free fall is expected to be of the level 
of ?7i2 ^ IQ^'^ X{Ai, Zi] A2, Z2)[d^lnB)Q. Since the factor X{Ai, Zi\ A2, Z2) typically ranges as 



82 



0(0.1 — 10), we deduce that {d^ lni3)o has to be very small for the Solar system constramts to be 
satisfied. It follows that today the scalar field has to be very close to the minimum of the coupling 
function Ini?. This is indeed very simplistic because we only take into account the effect of the 
electromagnetic binding energy (see § 16.31) . 



Let us also note that such a simple coupling cannot be eliminated by a conformal rescaling 
5/^1. = A^{(f)g% since 

j B{ct^)g'^''g^^''F,,F,,^A\ j B{^)A''-\cp)g^/g^:''F,,Fp,^A^x 

so that the action is invariant in Z? = 4 dimensions. 

This example shows that we cannot couple a field blindly to e.g. the Faraday tensor to make 
the fine-structure constant dynamics and that some mechanism for reconciling this variation with 
local constraints, and in particular the university of free fall, will be needed. 

5.2 High-energy theories and varying constants 
5.2.1 Kaluza-Klein 

Such coupling terms naturally appear when compactifying a higher-dimensional theory. As an 
example, let us recall the compactification of a 5-dimensional Einstein- Hilbert action (Ref. [404] . 
chapter 13) 

^^12^/^^'i'- 
Decomposing the 5-dimensional metric ijab as 



M V 



M ^ 



where M is a mass scale, we obtain 



s = tttV I \R' ttt^p' ) <^^d*x, (170) 



167rG* 




where the 4-dimensional gravitational constant is G* = 37rG5/4 / dy. The scalar field couples 
explicitly to the kinetic term of the vector field and cannot be eliminated by a redefinition of the 
metric: again, this is the well-known conformal invariance of electromagnetism in four dimensions. 
Such a term induces a variation of the fine-structure constant as well as a violation of the universal- 
ity of free-fall. Such dependencies of the masses and couplings are generic for higher-dimensional 
theories and in particular string theory. It is actually one of the definitive predictions for string the- 
ory that there exists a dilaton, that couples directly to matter |477j and whose vacuum expectation 
value determines the string coupling constants [542) . 

In the models by Kaluza |268) and Klein |290] the 5-dimensional spacetime was compactified 
assuming that one spatial extra-dimension iS*""^, of radius R^^- It follows that any field xi^^iU) 
can be Fourier transformed along the compact dimension (with coordinate y), so that, from a 
4-dimensional point of view, it gives rise to a tower of of fields x^"^^^'^) of "^n — nRxK- At 
energies small compared to R^j^ only the y-independent part of the field remains and the physics 
looks 4-dimensional. 
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Assuming that the action ()170p corresponds to the Jordan frame action, as the couphng cjiR 
may suggest, it fohows that the gravitational constant and the Yang-Mills coupling associated with 
the vector field A^^ must scale as 



Goc 



-1 



9yIi °^ 'f'/Gcc (f'. 



(171) 



Note that the scaling of G with (f) (or time) is not the one of the gravitational constant that would 



be measured in a Cavendish experiment since Eq. pMl) tells us that Gcav cx G*0"M 1 + 



where the constants Ki depends only on the dimension and topology of the compact space |522| 
so that the only fundamental constant of the theory is the mass scale M/^j^d entering the 4 + 
/^-dimensional theory. A theory on A^4 x A^z? where M-o \s a. /^-dimensional compact space 
generates a low-energy quantum field theory of the Yang-Mills type related to the isometrics of 
M.D [for instance Ref. [541] showed that for _D = 7, it can accommodate the Yang-Mills group 
SU{'i) X SU{2) X U{1)]. The two main problems of these theories are that one cannot construct 
chiral fermions in four dimensions by compactification on a smooth manifold with such a procedure 
and that gauge theories in five dimensions or more are not renormalisable. 

In such a framework the variation of the gauge couplings and of the gravitational constant arises 
from the variation of the size of the extra-dimensions so that one can derives stronger constraints 
that by assuming independent variation, but at the expense of being more model-dependent. Let 
us mention the works by Marciano |343J and Wu and Wang [545] in which the structure constants 
at lower energy are obtained by the renormalisation group. 

Ref. [294] used the variation ()172|) to constrain the time variation of the radius of the extra- 
dimensions during primordial nucleosynthesis to conclude thatj Ai?^^/i?^^ | < 1%. Ref. [55] 
took the effects of the variation of oc R:^^ and deduced from the helium-4 abundance that 
\AR^^/R^^\ < 0.7% and \AR^^/R^J < 1.1% respectively for D = 2 and D = 7 Kaluza-Klein 
theory and that |Ai?j^j^/i?j^j^| < 3.4 x 10"^" from the Oklo data. An analysis of most cosmological 
data (BBN, CMB, quasar etc..) assuming that the extra-dimension scales as -Ro(l + A.t~^^'^) and 
Ro[l + A]{l — cos uj{t — to) concluded that A has to be smaller tha 10~^^ and 10~* respectively [304] . 
while Ref. [328] assumes that gauge fields and matter fields can propagate in the bulk. Ref. [334] 
evaluated the effect of such a couple variation of G and the structures constants on distant super- 
nova data, concluding that a variation similar to the one reported in Ref. [519] would make the 
distant supernovae brighter. 

5.2.2 String theory 

There exist five anomaly free, supersymmetric perturbative string theories respectively known as 
type I, type IIA, type JIB, SO(32) heterotic and Eg x Es heterotic theories (see e.g. Ref. 415 ). 
One of the definitive predictions of these theories is the existence of a scalar field, the dilaton, that 
couples directly to matter [477 and whose vacuum expectation value determines the string coupling 
constant [542] . There are two other excitations that are common to all perturbative string theories, 
a rank two symmetric tensor (the graviton) g^^ and a rank two antisymmetric tensor Bfj^^. The 
field content then differs from one theory to another. It follows that the 4-dimensional couplings 
are determined in terms of a string scale and various dynamical fields (dilaton, volume of compact 
space, . . . ). When the dilaton is massless, we expect three effects: (i) a scalar admixture of a scalar 
component inducing deviations from general relativity in gravitational effects, (ii) a variation of 
the couplings and (iii) a violation of the weak equivalence principle. Our purpose is to show how 




This can be generalized to the case of D extra-dimensions [113] to 



G cx 



D 



(172) 
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the 4-dimensional couplings are related to the string mass scale, to the dilaton and the structure 
of the extra-dimensions mainly on the example of heterotic theories. 

To be more specific, let us consider an example. The two heterotic theories originate from the 
fact that left- and right-moving modes of a closed string are independent. This reduces the number 
of supersymmetry to iV = 1 and the quantization of the left-moving modes imposes that the gauge 
group is either 50(32) or Eg, x Eg depending on the fermionic boundary conditions. The effective 
tree-level action is 



2* 



Ml {i?io + 4n$ - 4(V$)2} - -^FabF 



AB 



(173) 



When compactified on a 6-dimensional Calabi-Yau space, the effective 4-dimensional action takes 
the form 



Sh = 



Ml 



H p2 



(174) 



where (p = Vge couples identically to the Einstein and Yang-Mills terms. It follows that 



Ml 



Mt 



,-2 



(175) 



at tree-level. Note that to reach this conclusion, one has to assume that the matter fields (in the 
'dots' of Eq. (|174l) are minimally coupled to (74; see e.g. Ref. |337| ). 

The strongly coupled SO(32) heterotic string theory is equivalent to the weakly coupled type I 
string theory. Type I superstring admits open strings, the boundary conditions of which divide the 
number of supersymmetries by two. It follows that the tree-level effective bosonic action is iV = 1, 
13 = 10 supergravity which takes the form, in the string frame. 



Si = y"dl"x^/^Mfe- 



e-*M?i?i 



El 
4 



(176) 



where the dots contains terms describing the dynamics of the dilaton, fermions and other form 
fields. At variance with (|173p . the field $ couples differently to the gravitational and Yang- Mills 
terms because the graviton and Yang-Mills fields are respectively excitation of close and open 
strings. It follows that Mj can be lowered even to the weak scale by simply having exp <I> small 
enough. Type I theories require £)9-branes for consistancy. When Vq is small, one can use T- 
duality (to render Vg large, which allows to use a quantum field theory approach) and turn the 
Z39-brane into a Z33-brane so that 



Si = I d^"xy^e-2*Mfi? 



(177) 



where the second term describes the Yang-Mills fields localized on the £)3-brane. It follows that 



Mf = e 



-2* 



9^ 



(178) 



at tree-level. If one compactifies the Z39-brane on a 6-dimensional orbifold instead of a 6-torus, 
and if the brane is localized at an orbifold fixed point, then gauge fields couple to fields Mi living 
only at these orbifold fixed points with a (calculable) tree-level coupling Ci so that 



9^ 



+ c^M,. 



(179) 



The coupling to the field Ci is a priori non universal. At strong coupling, the 10-dimensional Eg x Eg 
heterotic theory becomes M-theory on i?^" x S^/Z2 |254| . The gravitational field propagates in 
the 11-dimensional space while the gauge fields are localized on two 10-dimensional branes. 
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At one-loop, one can derive the couplings by including Kaluza-Klein excitations to get |162] 



5. 



YM 



2 



H 



(180) 



when the volume is large compared to the mass scale and in that case the coupling is no more 
universal. Otherwise, one would get a more complicated function. Obviously, the 4-dimensional 
effective gravitational and Yang-Mills couplings depend on the considered superstring theory, on 
the compactification scheme but in any case they depend on the dilaton. 

As an example, Ref. [337] considered the {N = 1, = 10)-supergravity model derived from the 
heterotic superstring theory in the low energy limit and assumed that the 10-dimensional spacetime 
is compactified on a 6-torus of radius R{x^) so that the effective 4-dimcnsional theory described 
by (|174p is of the Brans-Dicke type with oj = —1. Assuming that (j) has a mass /i, and couples 
to the matter fluid in the universe as S^^tt^ = / d^^x-^— gio cxp{—2^)£^^^^^^^ {gio), the reduced 
4-dimensional matter action is 



The cosmological evolution of (/> and R can then be computed to deduce that Q!em/q^em — 10^° 
eV)~^ yr~^. Ref. considered the same model but assumed that supersymmetry is broken by 
non-perturbative effects such as gaugino condensation. In this model, and contrary to Ref. |337| , 4> 
is stabilized and the variation of the constants arises mainly from the variation of i? in a runaway 
potential. 

To conclude, superstring theories offer a natural theoretical framework to discuss the value of 
the fundamental constants since they become expectation values of some fields. This is a first step 
towards their understanding but yet, no complete and satisfactory mechanism for the stabilization 
of the extra-dimensions and dilaton is known. 

It has paved the way to various models that we detail in § 15.41 

5.3 Relations between constants 

There are different possibilities to relate the variations of different constants. First, in quantum 
field theory, we have to take into account the running of coupling constants with energy and the 
possibilities of grand unification to relate them. It will also give a link between the QCD scale, the 
coupling constants and the mass of the fundamental particles (i.e. the Yukawa couplings and the 
Higgs vev). Second, one can compute the binding energies and the masses of the proton, neutron 
and different nuclei in terms of the gauge couplings and the quark masses. This step involves QCD 
and nuclear physics. Third, one can relate the gyromagnetic factor in terms of the quark masses. 
This is particularly important to interpret the constraints from the atomic clocks and the QSO 
spectra. This allows to set stronger constraints on the varying parameters at the expense of a 
model-dependence. 

5.3.1 Implication of gauge coupling unification 

The first theoretical implication of high-energy physics arises from the unification of the non- 
gravitational interactions. In these unification schemes, the three standard model coupling con- 
stants derive from one unified coupling constant. 

In quantum field, the calculation of scattering processes include higher order corrections of 
the coupling constants related to loop corrections that introduce some intergrals over internal 
4-momenta. Depending on the theory, these integrals may be either finite or diverging as the 




(181) 
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logarithm or power law of a UV cut-off. In a class of theories, called renormalizable, among which 
the standard model of particule physics, the physical quantities calculated at any order do not 
depend on the choice of the cut-off scale. But the result may depend on \nE/m where E is the 
typical energy scale of the process. It follows that the values of the coupling constants of the 
standard model depend on the energy at which they are measured (or of the process in which 
thay are involved). This running arises from the screening due to the existence of virtual particles 
which are polarized by the presence of a charge. The renomalization group allows to compute the 
dependence of a coupling constants as a function of the energy E as 

dii^-^^^^)' 

where the beta functions, depend on the gauge group and on the matter content of the theory 
and may be expended in powers of Qi. For the SU(2) and U(l) gauge couplings of the standard 
model, they are given by 



47r2 V 6 3 y ' 47r2 9 

where rig is the number of generations for the fermions. We remind that the fine-structure constant 
is defined in the limit of zero momentum transfer so that cosmological variation of ctj^.^, arc inde- 
pendent of the issue of the renormalisation group depence. For the SU(3) sector, with fundamental 
Dirac fermion representations, 

/33(53)--^(^^--^ 

nf being the number of quark flavours with mass smaller than E. The negative sign implies 
that (1) at large momentum transfer the coupling decreases and loop corrections become less and 
less significant: QCD is said to be asymptotically free; (2) integrating the renormalisation group 
equation for as gives 

"'^^^ " {33 -nf) HE /A,) 

so that it diverges as the energy scale approaches Ac from above, that we decided to call Aqcd- This 
scale characterises all QCD properties and in particular the masses of the hadrons are expected to 
be proportional to Aqcd up to corrections of order mq/Aqco- 

It was noticed quite early that these relations imply that the weaker gauge coupling becomes 
stronger at high energy, while the strong coupling becomes weaker so that one can thought the 
three non-gravitational interactions may have a single common coupling strength above a given 
energy. This is the driving idea of Grand Unified Theories (GUT) in which one introduces a 
mechanism of symmetry-breaking from a higher symmetry group, such e.g. as SO(IO) or SU(5), 
at high energies. It has two important consequences for our present considerations. First there 
may exist algebraic relations between the Yukawa couplings of the standard model. Second, the 
structure constants of the standard model unify at an energy scale 

a,(Mu) = a2(Mu) - a^iMu) = au{Mu). (182) 

We note that the electroweak mixing angle, i.e. the can also be time dependent parameter, but 
only for E ^ Mu since at E = Mu, it is fixed by the symmetry to have the value sin^ 9 = 3/8, 
from which we deduce that 
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It follows from the renormalisation group relations that 

a-\E) = a~\Mu) - (183) 

where the beta-functfon coefficients are given by hi — (41/10, —19/6, 7) for the standard model (or 
below the SUSY scale Asusy) and by hi = (33/5, 1, —3) for = 1 supersymmetric theory. Given 
a field decoupling at mth, one has 

where bf^'' = 6^+^ — 6'-"^ with 6^+/") the beta-function coefficients respectively above and below 
the mass threshold, with tree-level matching at mth- In the case of multiple thresholds, one must 
sum the different contributions. The existence of these thresholds implies that the running of 
is complicated since it depends on the masses of heavy quarks and coloured superpartner in the 
case of supersymmetry. For non-supersymmetric theories, the low-energy expression of the QCD 
scale is 

(=^)"'-(-5^) 

for E > mt- This implies that the variation of Yukawa couplings, gauge couplings, Higgs vev and 
Aqcd/A^p are correlated. A second set of relations arises in models in which the weak scale is 
determined by dimensional transmutation |186l 1185] . In such cases, the Higss vev is related to the 
Yukawa constant of the top quark by 

z; = AfpCxp(^ ^ j , (185) 

where c is a constant of order unity. This would imply that 5\nv = SSlnh with S ^ 160 |104j . 

The first consequences of this unification were investigated in Refs. [73l [TH I311J where the 
variation of the 3 coupling constants was reduced to the one of au and My /Mp . It was concluded 
that, setting 

i? = 51nAQCD/<51na^„, (186) 

i? ~ 34 with a stated accuracy of about 20% |310|l311j (assuming only au can vary), i? = 38±6 [73] 
and then i? = 46 [74l [75], the difference arising from the quark masses and their associated 
thresholds. However, these results implicitely assume that the electroweak symmetry breaking 
and supersymmetry breaking mechanisms, as well as the fermion mass generation, are not affected 
by the variation of the unified coupling. It was also mentioned in Ref. |74j that R can reach 
—235 in unification based on SU(5) and SO(IO). The large value of R arises from the exponential 
dependence of Aqcd on a^. In the limit in which the quark masses are set to zero, the proton mass, 
as well as all other hadronic masses are proportional to Aqcd, i.e. nip cx Aqci3(1 + 0(111^/ Aqq^))). 
Ref. [311] further relates the Higgs vev to a^.^ by d In w/d In a^.^ = k and estimated that k ^ 70 
so that, assuming that the variation of the Yukawa couplings is negligible, it could be concluded 
that 

TTl 

5 In- --35^1nQ;EM, 

Aqcd 

for the quark and electron masses. This would also implies that the variation of and a^^jj^ are 
correlated, still in a very model-dependent way, typically one can conclude [104) that 

^^_0.8i?^-l-0.6(^-fl)^, 
M "em it- 
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with S ^ 160. The runnmg oiajj can be extrapolated to the Planck mass, Mp. Assumiung au{Mp) 
fixed and letting Mu/Mp vary, it was concluded [TSl] that R = 2n{bu + 3)/ [9a^„ (86(7/3 - 12)] 
where bu is the beta-function coefficient describing the running of aij. This shows that a variation 
of ckg^ and /x can open a windows on GUT theories. A similar analysis |141j assuming that elec- 
troweak symmetry breaking was triggered by nonperturbative effects in such a way that v and ajj 
are related, concludes that ^/x/// = (13 ± 7)(5aj,M/<^EM ^ theory with soft SUSY breaking and 
SiJ,/^ = (— 4 ± 5)(5Q;j,„/ajj,„ otherwise. 

From a phenomenological point of view, Ref. |146) making an assumption of proportionality 
with fixed "unification coefficients" assumes that the variations of the constants at a given redshift 
z depend on a unique evolution factor £{z) and that the variation of all the constants can be derived 
from those of the unification mass scale (in Planck units), Mu, the unified gauge coupling au, the 
Higgs vev, V and in the case of supersymmetric theories the soft supersymmetry breaking mass, 
m. Introducing the coefficients di by 

Aln§^ = dM^, Alnau = dui, Aln^^dni, Aln^ =dsi, 

(ds = for non-supersymmetric theories) and assuming that the masses of the standard model 
fermions all vary with v so that the Yukawa couplings are assumed constant, it was shown that the 
variations of all constants can be related to {dM,dij,dH,ds) and £{z), using the renormalisation 
group equations (neglecting the effects induced by the variation of au on the RG running of fermion 
masses). This decomposition is a good approximation provided that the time variation is slow, 
which is actually backed up by the existing constraints, and homogeneous in space (so that it may 
not be applied as such in the case a chameleon mechanism is at work [67 ). 

This allowed to define 6 classes of scenarios: (1) varying gravitational constant {dn — ds — 
dx = 0) in which only AIjj/Mp or equivalently G'Aq(-,j-, is varying; (2) varying unified coupling 
{du — Ijdn — ds — dpi — 0); (3) varying Fermi scale defined by (dn — i,dij — ds ~ d^i — 
0) in which one has dln/i/dlna^.^ = —325; (4) varying Fermi scale and SUSY-breaking scale 
{ds = du = ^,du = dM = 0) and for which dln/i/dlnaj.^, = —21.5; (5) varying unified coupling 
and Fermi scale {dx = l,rf// = jdx^ds = dM ~ 0) and for which dln/i/dlnaj.^, = (23.2 — 
0.657)7(0.865 + O.O27); (6) varying unified coupling and Fermi scale with SUSY {dx — \,ds — 
dii = ^dx,dM = 0) and for which dln^/dlna^^ = (14 - 0.287)/(0.52 + O.OI37). 

Each scenario can be compared to the existing constraints to get sharper bounds on them |145[ 
I146| 1147} 1361) and emphasize that the correlated variation between different constants (here /i and 
ajj,„) depends strongly on the theoretical hypothesis that are made. 



5.3.2 Masses and binding energies 

The previous section described the unification of the gauge couplings. When we consider "com- 
posite" systems such as proton, neutron, nuclei or even planets and stars, we need to compute 
their mass, which requires to determine their binding energy. As we have already seen, the elec- 
tromagnetic binding energy induces a direct dependence on a^.^, and can be evaluated using e.g. 
the Bethe-Weizacker formula (|60p . The dependence of the masses on the quark masses, via nu- 
clear interactions, and the determination of the nuclear binding energy are especially difficult to 
estimate. 

In the chiral limit of QCD in which all quark masses are negligible compared to Aqcd all 
dimensionful quantities scale as some power of Aqcd- For instance, concerning the nucleon mass, 
mN = cAqcd with c '-^ 3.9 being computed from lattice QCD. This predicts a mass of order 
860 MeV, smaller than the observed value of 940 MeV. The nucleon mass can be computed in 
chiral perturbation theory and expressed in terms of the pion mass as [314] ton = + a2"^^ + 
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a^m^ + aQm% + aN-n + ctatt + ctad (where all coefBcients of this expansion are defined in Ref. I 314j ) , 
which can be used to show [213] that the nucleon mass is scaling as 

mNOcAQCDX°-°'X°-°"- (187) 

It was further extanded [207] by using a sigma model to infer that ttin oc AQCD-^'^q ^^^-'^""'^^- This 
two examples explicitely show the strong dependence in the nuclear modelling. 

To go further and determine the sensitivity of the mass of a nucleus to the various constant, 

m(y4, Z) = Znip + {A — Z)mn + Znic + + E^^ 

one should determine the strong binding energy [see related discussion below Eq. (I16p ] in function 
of the atomic number Z and the mass number A. 

The case of the deuterium binding energy B^i has been discussed in different ways (see § I3.8.3P . 
Many modelisations have been performed. A first route relies on the use of the dependence of 
on the pion mass |189| [37l I420[ 1549 1 . which can then be related to mu, and Aqcd- A second 
avenue is to use a sigma model in the framework of the Walecka model [451] in which the potential 
for the nuclear forces keeps only the cr, p and uj meson exchanges |207j . We also emphasize that 
the deuterium is only produced during BBN, as it is too weakly bound to survive in the regions 
of stars where nuclear processes take place. The fact that we do observe deuterium today sets a 
non-trivial constraint on the constants by imposing that the deuterium remains stable from BBN 
time to today. Since it is weakly bound, it is also more sensitive to a variation of the nuclear force 
compared to the electromagnetic force. This was used in Ref. |144j to constrain the variation of 
the nuclear strength in a sigma-model. 

For larger nuclei, the situation is more complicated since there is no simple modelling. For 
large mass number A, the strong binding energy can be approximated by the liquid drop model 

E as {A-2Zf + 

with(ay,os,aA,op) = (15.7, 17.8,23.7, 11.2) MeV [433]. It has also been suggested [126] that the 
nuclear binding energy can be expressed as 

E^ ~ Aa^ + A^'H^ with ag = af + (189) 

In the chiral limit, 03 has a non- vanishing limit to which we need to add a contribution scaling 
like cx AqcD^q. Ref. [126 also pointed out that the delicate balance between attractive and 
repulsive nuclear interactions [451, implies that the binding energy of nuclei is expected to depend 
strongly on the quark masses 1158] . Recently, a fitting formula derived from effective field theory 
and based of the semi-empirical formula derived in Ref. |223l was proposed [119] as 

where rjs and rjy are the strength of respectively the scalar (attractive) and vector (repulsive) 
nuclear contact interactions normalized to their actual value. These two parameters need to be 
related to the QCD parameters [158] . We also refer to Ref. |211] for the study of the dependence 
of the binding of light [A < 8) nuclei on possible variations of hadronic masses, including meson, 
nucleon, and nucleon-resonance masses. 

These expressions allow to compute the sensitivity coefficients that enter in the decomposition 
of the mass [see Eq. ()200l) ]. They also emphasize one of the most difficult issue concerning the 
investigation about constant related to the intricate structure of QCD and its role in low energy 
nuclear physics, which is central to determine the masses of nuclei and the binding energies, 
quantities that are particularly important for BBN, the universality of free fall and stellar physics. 
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5.3.3 Gyromagnetic factors 

The constraints arising from the comparison of atomic clocks (see § 13. ip involve the fine-structure 
constant a^^ , the proton-to-electron mass ratio /x and various gyromagnetic factors. It is important 
to relate these factors to fundamental constants. 

The proton and neutron gyromagnetic factors are respectively given by gp = 5.586 and — 
—3.826 and are expected to depend on = mq/Aqco [199] . In the chiral limit in which = 
rrid = 0, the nucleon magnetic moments remain finite so that one could have thought that the finite 
quark mass effects should be small. However, it is enhanced by 7r-meson loop corrections which 
are proportional to cx a/to^^Aqcd- Following Ref. |314] . this dependence can be described by 
the approximate formula 

1 + arriTr + brn^ 

The coefficients are given by a = (1.37, 1.85)/GeV and b = (0.452, 0.271)/GeV2 respectively for 
the proton an neutron. This lead 199 to cx m~°-^'''' cx and cx m~^ '^^^ cx 

This was further extended in Ref. 213 to take into the depence with the strange quark mass mg 
to obtain 

5p cx g„ cx (igi) 

This allows to express the results of atomic clocks (see § 13.1.31) in terms of c^emi ^q; -^s and 
Xc- Similarly, for the constants constrained by QSO observation, we have (see Table [T0|) 
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(192) 

once the scaling of the nucleon mass as mN cx AqcbX^'^^'^ X^ '^^^ (see ? I5.3.2p . This shows that 
the 7 observable quantities that are constrained by current QSO observations can be reduced to 
only 4 parameters. 

5.4 Models with varying constants 

The models that can be constructed are numerous and cannot be all reviewed here. We thus focus 
in the string dilaton model in j; l5.4.1l and then discuss the chameleon mechanism in i; l5.4.2l and the 
Bekcnstein framework in § 15.4.31 

5.4.1 String dilaton and Runaway dilaton models 

Damour and Polyakov [130) I131j argued that the effective action for the massless modes taking 
into account the full string loop expansion should be of the form 

S = 



J d*xy/^ A/2 |Bg($)fl + 4B$($) □$-(V$)2 j-BF($)^F2 
^B4<^>)^P4; + ..] (193) 
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in the string frame, Ms being the string mass scale. The functions Bi are not known but can be 
expanded (from the genus expansion of string theory) in the hmit $ — cx) as 



-2<S> 



+ C 



4^)e4* 



(194) 



where the first term is the tree level term. It follows that these functions can exhibit a local 
maximum. After a conformal transformation {gp.^ = CBgg^^^ip = (CBg)~^^'^B^'^'ij}), the action 
in Einstein frame takes the form 

d^x . 
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(195) 



where the field <b is defined as 
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It follows that the Yang-Mills coupling behaves as g^^ = kBp{(t)). This also implies that the QCD 
mass scale is given by 
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(196) 



where b depends on the matter content. It follows that the mass of any hadron, proportional to 
Aq(,j5 in first approximation, depends on the dilaton, mA{Bg, Bp, . . .). 

If, as allowed by the anstaz (|194p . mA{4>) has a minimum 4>m then the scalar field will be driven 
toward this minimum during the cosmological evolution. However if the various coupling functions 
have different minima then the minima of mA{4>) will depend on the particle A. To avoid violation 
of the equivalence principle at an unacceptable level, it is thus necessary to assume that all the 
minima coincide in (/> = (jjrm which can be implemented by setting Bi — B. This can be realized by 
assuming that (j)^ is a special point in field space, for instance it could be associated to the fixed 
point of a Z2 symmetry of the T- or S'-duality [126] . 

Expanding ln_B around its maximum (j>m as InB cx — k(0 — 0m)^/2, Damour and Polyakov |130[ 
1131) constrained the set of parameters (k, (t)o — 4>m) using the different observational bounds. This 
toy model allows one to address the unsolved problem of the dilaton stabilization, to study all the 
experimental bounds together and to relate them in a quantitative manner (e.g. by deriving a link 
between equivalence-principle violations and time- variation of a^.^ ) . This model was compared to 
astrophysical data in Ref. [308] to conclude that |A(/)| < 3.4k10~^. 

An important feature of this model lies in the fact that at lowest order the masses of all nuclei 
are proportional to Aqcd so that at this level of approximation, the coupling is universal and the 
theory reduces to a scalar-tensor theory and there will be no violation of the universality of free 
fall. It follows that the deviation from general relativity are characterized by the PPN parameters 
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with A0O — fpa ~ 4'm and /3s ^ 40 [130j . The variation of the gravitational constant is, from 
Eq. (|166p . simply 
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The value of (po — i?o0o is obtained from the Klein-Gordon equation (|167|) and is typically given 
by 00 — —Z l3sKHQA(j)o were Z is a number that depends on the equation of state of the fluid 
dominating the matter content of the universe in the last e-fold and the cosmological parameters 
so that 

^ = 2fA^o = -2ZHol3^,K.^A<j)l (198) 





G 



is to use the Friedmann equations 



The factor Z is model-dependent and another way to estimate ( 
which imply that (po — Hq^ 1 + qo — |5^mO where q is the deceleration parameter. 

When one considers the effect of the quark masses and binding energies, various composition- 
dependent effects appear. First, the fine-structure constant scales as a,^^ ~ so that 
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The second effect is, as pointed out earlier, a violation of the universality of free fall. In full 
generality, we expect that 



mA(0) = A^Aqcd(0) 
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Using an expansion of the form (|16|) , it was concluded that 
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with B = N + Z,D = N-Z anAE = Z{Z -I) / {N + Zf/'^ and where the value of the parameters 
d are model-dependent. 

It follows from this model that: 

• The PPN parameters, the time variation of a and G today and the violation of the university 
of free-fall all scale as A0q. 

• The field is driven toward 0m during the cosmological evolution, a point at which the scalar 
field decouples from the matter field. The mechanism is usually called the least coupling 
principle. 

• Once the dynamics for the scalar field is solved, A0o can be related to A</)i at the end of 
inflation. Interestingly, this quantity can be expressed in terms of amplitude of the density 
contrast at the end of infiation, that is to the energy scale of infiation. 



The numerical estimations |130] indicate that r}u,H ~ —5.4 x 10^^ 
in such a class of models, the constraint on 10"^'^ implies 1 — 7 
better constraint that the one obtained direcly. 
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This model was extended }135j the a case where the coupling functions have a smooth finite 
limit for infinite value of the bare string coupling, so that Hi = Gi + ©(e^"^), foiling Ref. [230] . 
The dilaton runs away toward its attractor at infinity during a stage of inflation. The late time 
dynamics of the scalar field is similar as in quintessence models, so that the model can also explain 
the late time acceleration of the cosmic expansion. The amplitude of residual dilaton interaction 
is related to the amplitude of the primordial density fluctuations and it induces a variation of the 
fundamental constants, provided it couples to dark matter or dark energy. It is concluded that, in 
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this framework, the largest ahowed variation of is of order 2 x 10 ^, which is reached for a 
violation of the universality of free fall of order 10^^^ and it was established that 

~ ±10-1^1 + go - ^n^oV^O^yT-\ (202) 

where the first square-root arises from the computation of (jjQ. The formalism was also used to 
discuss the time variation of 0;^,^, and ^ [96] . 

The coupling of the dilaton to the standard model fields was further investigated in Refs. |120[ 
1121) . Assuming that the heavy quarks and weak gauge bosons have been integrated out and that 
the dilaton theory has been matched to the light fields below the scale of the heavy quarks, the 
coupling of the dilaton has been parameterised by 5 parameters: de and dg for the couplings to 
the electromagnetic and gluonic field-strength terms, and dm^ , dm„ and dm^ for the couplings to 
the fermionic mass terms so that the interaction Lagrangian is reduces to a linear coupling (e.g. 
oc de4>F'^ for the coupling to electromagnetism etc.) It follows that lS.a-^^/a-^„ — deK.4> for the 
fine structure constant, AAqcd/Aqcd = ddK<t> for the strong sector and ^milmi = dmiK4> for 
the masses of the fermions. These parameters can be constrained by the test of the equivalence 
principle in the Solar system [see § 16.3] . 

In these two string- inspired scenarios, the amplitude of the variation of the constants is related 
to the one of the density fluctuations during inflation and the cosmological evolution. 

5.4.2 The Chameleon mechanism 

A central property of the least coupling principle, that is at the heart of the former models, is 
that all coupling functions have the same minimum so that the effective potential entering the 
Klein-Gordon equation for the dilaton has a well-defined minimum. 

It was realized [287] that if the dilaton has a coupling A'^{(j)) to matter while evolving in a 
potential V{(j>) the source of the Klein-Gordon equation (|167p has a an effective potential 

V,g = Vi^) + A^icl))p. 

In the case where y is a decreasing function of 4>, e.g. a runaway potential, and the coupling is 
an increasing function, e.g. = exp/30/Mp, the effective potential has a minimum the value of 
which depends on the matter local density p (see also Ref. [187] ). The field is thus expected to be 
massive on Earth where the density is high and light in space in the Solar system. It follows that 
the experiment on the universality of free fall in space may detect violations of the universality of 
free fall larger than the bounds derived by laboratory experiments [2881 1493] . It follows (1) that 
the constraints on the time variation of the constants today can be relaxed if such a mechanism is 
at work and (2) that is the constants depend on the local value of the chameleon field, their value 
will be environment dependent and will be different on Earth and in space. 

The cosmological variation of a^.^, in such model was investigated in Ref. [69j [70] . Models based 
on the Lagrangian (|208[) and exhibiting the chameleon mechanism were investigated in Ref. [393] . 

The possible shift in the value of n in the Milky Way (see § KT^ was related [211 IS13] 
to the model of Ref. |393| to conclude that such a shift was compatible with this model. 

5.4.3 Bekenstein and related models 

Bekenstein |3H1 ISS] introduced a theoretical framework in which only the electromagnetic sector 
was modified by the introduction of a dimensionless scalar field e so that all electric charges vary 
in unison Ci = eoie(a;") so that only a^.^ is assumed to possibly vary. 
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To avoid the arbitrariness in the definition of e, which can be rescaled by a constant factor 
while eoi is inversely rescales, it was postulated that the dynamics of e be invariant under global 
rescaling so that its action should be of the form 

he f g^'^df.ed^e — 4 
= / % V-5d X, (203) 



2^2 

I being a constant length scale. Then, e is assumed to enter all electromagnetic interaction via 
CiA^ — >■ e^itA^ where is the usual electromagnetic potential and the gauge invariance is then 
preserved only if eA^ eA^ + A^^ for any scalar function A. It follows that the the action for the 
electromagnetic sector is the standard Maxwell action 

Se = --^ / F^^'^F^^^d^x, (204) 



167r _ 

for the generalized Faraday tensor 

F^,. = i [(eA.),^ - {eA^)..] (205) 

To finish the gravitational sector is assumed to be described by the standard Einstein-Hilbert 
action. Finally, the matter action for point particles of mass m takes the form Sm = X] J[—'mc^ + 
{e/c)u^'A^]-f-'^S^{x' - a;*(T))d''x where 7 is the Lorentz factor and t the proper time. Note that 
the Maxwell equation becomes 

{e-'F^'') = A7:f (206) 

which is the same as for electromagnetism in a material medium with dielectric constant and 
permeability (this was the original description proposed by Fierz [196 and Lichnerowicz [330) : 
see also Dicke [151]). 

It was proposed [439] to rewrite this theory by introducing the two fields 



af, = eAf^, 'i/' = lne 



so that the theory takes the form 



^^1^9 I ^^'^'^ ^-'"^rU^V^d'^ (207) 

with K = I / [At: he) and /^j^ the Faraday tensor associated with a^. The model was further extended 
to include a potential for V' [30] and to include the electroweak theory [4571 . 

As discussed previously, this class of models predict a violation of the universality of free fall and, 
from Eq. p3|) . it is expected that the anomalous acceleration is given by (5a — —M~^{dEj^^/de)\7e. 

From the confrontation of the local and cosmological constraints on the variation of e Beken- 
stein |38] concluded, given his assumptions on the couplings, that a^^, "is a parameter, not a 
dynamical variable" (see however Ref. [39j and then Ref. [299) ) . This problem was recently by- 
passed by Olive and Pospelov [392] who generalized the model to allow additional coupling of a 
scalar field — Bp{(l)) to non-baryonic dark matter (as first proposed in Ref. |134) ) and cosmo- 
logical constant, arguing that in supersymmetric dark matter, it is natural to expect that would 
couple more strongly to dark matter than to baryon. For instance, supersymmetrizing Bekenstein 
model, 4> will get a coupling to the kinetic term of the gaugino of the form M~^<l)xdx- Assuming 
that the gaugino is a large fraction of the stable lightest supersymmetric particle, the coupling to 
dark matter would then be of order 10^ — 10^ times larger. Such a factor could almost reconcile 
the constraint arising from the test of the universality of free fall with the order of magnitude of 
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the cosmological variation. This generalization of the Bekenstein model relies on an action of the 
form 



S = -Mi I R^d''^ 



-5d''x (208) 



- / {E ^«[*^- m,BNM]N^ + \xdx + MlBK{(f)K + ]^M^B^{4,)x^ x^ T^ffd^x 

where the sum is over proton [p = ^^{d^ — ieoA^)] and neutron [p = 7^9^]. The functions B 
can be expanded (since one focuses on small variations of the fine-structure constant and thus 
of (j)) as Bx = 1 + Cx0 + Cjf^^/S. It follows that a^^{(l)) = eQ/AirBpict)) so that /S.a^^/a^„ = 
(f4'+{£,f—'^Cf)<I>'^ f"^- This framework extends the analysis of Ref. [38] to a 4-dimensional parameter 
space (M*, C_F, Cm, Ca)- It contains the Bekenstein model (Cf = —2, Ca = Cm ^ IO^^'Cf), a 
Jordan-Brans-Dicke model (Cf = 0, Ca = —2^2/2uj + 3, Cm = — l/-\/4a; + 6), a string-like model 
(Cf = — \/2, Ca = "^/2, Cm ~ V2/2) so that A/ctjj.^ = 3) and supersymmetrized the Bekenstein 
model (Cf = —2, Cx = ~2, Cm = Cx that Aa^^/a^.^ ^ 5/a;). In all the models, the universality 
of free fall sets a strong constraint on Cf/a/w (with uj = M^/2Mi) and the authors showed that 
a small set of models may be compatible with a vairation of a^,^ from quasar data while being 
compatible the equivalence principle tests. A similar analysis |349j concluded that such model can 
reproduce the variation of 0!^,^, from quasar while being compatible with Oklo and meteorite data. 
Note that under this form, the effective theory is very similar to the one detailed in § 15.4.21 

This theory was also used [40l to study the spacetime structure around charged black-hole, 
which corresponds to an extension of dilatonic charged black hole. It was concluded that a cosmo- 
logical growth of a^^ would decrease the black-hole entropy but with half the rate expected from 
the earlier analysis [1381 1339] . 



5.4.4 Other ideas 

Let us mention without details other theoretical models which can accomodate varying constants: 

• Models involving a late time phase transition in the electromagnetic sector [86l [11] ; 

• Brancworld models [3341 [Hi [ID [329] [398] or extra-dimensions [473] : 

• Model with pseudo-scalar couplings |212| : 

• Growing neutrino models [9] 1529] in which the neutrino masses are a function of a scalar 
field, that is also responsible for the late time acceleration of the universe. In these models 
the neutrinos freeze the evolution of the scalar field when they become non-relativistic while 
its evolution is similar as in quintessence when the neutrinos are ultra-relativistic; 

• Models based on discrete quantum gravity |225j or on loop quantum gravity in which the 
Barbero-Immirzi parameter controls the minimum eigenvalue of the area operator and could 
be promoted to a field, leading to a classical coupling of Einstein gravity with a scalar-field 
stress-energy tensor [352, 487 

• "varying speed of light" models for which we refer to the review [338] and our previous 
analysis |184) for a critical view; 

• Quintessence models with a non-minimal coupling of the quintessence field [30] [TO] [55] HIH 

IIII1I2T1[SI21[3I1I3H5]IM]IM1[S2H] [see discussion § MM ; 

• Holographic dark energy models with non-minimal couplings |236j 
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6 Spatial variations 



The constraints on the variation of the fundamental constant that we have described so far are 
mainly related to their cosmological evolution so that, given the Copernican principle, they reduce 
to constrains on the time variation of the fundamental constants. Indeed, spatial variations can 
also occur. They may be used to set constraints in two regimes: 

• On cosmological scales, the fields dictating the variation of the constants have fluctuations 
that can let their imprint in some cosmological observables. 

• On local scales (e.g. Solar system or Milky Way) the fields at the origin of the variation 
of the constants are sourced by the local matter distribution so that one expect that the 
constants are not homogeneous on these scales. 

6.1 Local scales 

In order to determine the profile of the constant in the Solar system, let us assume that their value 
is dictated by the value of a scalar field. As in § 15.4.11 we can assume that at lowest order the 
profile of the scalar field will be obtained from the scalar-tensor theory, taking into account that 
all masses scale as Aqcd (0* ) where 0* is the value of the field in the Einstein frame. 

6.1.1 Generalities 

We restrict to the weakly self-gravitating (V*/c^ ^ 1) and slow moving {T^^ ^ yoo-j localized 
material systems and follow Ref. |123| . Using harmonic coordinates, defined with respect to the 
metric the Einstein frame metric can be expanded as 

50*0 = - exp (^-2^^ + 0(c-6), g*^ = -±V* + 0{c''% g*. = - exp {2^^ % + 0{c-^), 
so that Eqs. (|16mi62p take the form 

= -47rG,f7* -KO(c-''), D^K' = -47rG*cr: 0(c-2), D*,^* = -47rG,4 + C(c"^) (209) 

where □* is the fiat d'Alembertian and where the scalar field has been assumed to be light so 
that one can neglect its potential. The source terms can be expressed in terms of the matter 
stress-energy tensor in the Einstein frame as 

Restricting to the static case with a single massive point source, the only non-vanishing source 
terms are cr*(r) = M*(5'^(r*) and S'(r) = — Q!(0*)M*(5'^(r*) so that the set of equations reduces to 
two Poisson equations 

C M 

= -47rG*A/*(5^(r*) + 0(c~*), A*0* = 47r^3-^<5^(r*) ©(c"^). (210) 

This set of equations can be solved in the of the retarded Green function. It follows that the 
Einstein frame gravitational potential is 14 (r*) = G^M^/r*. The equation for 0, can be solved 
iteratively, since at lowest order in G*/c^ it has solution 

0* = = 00 TrV^ir^). 



97 



This can be used to determine the Jordan frame metric and the variation of the scalar field in 
function of the Jordan frame coordinates. It follows that at lowest order the Newton potential and 
the scalar field are given by 

where we have neglected the corrections —a{(j)){(t> — (t>o) for the gravitational potential which, given 
the Solar system constraints on ao, is a good approximation. 

Now let us consider any constant ai function of 0. Its profile is thus given by ai{r) ~ ai{(j)Q) — 
aQai{(j)o)^N{r)/c'^ so that 

(r) = -Si((t)o)ao 5 — (212) 

where Si((/'o) is the sensitivity of the constant to a variation of the scalar field, Si = dlnai/di/). 
For laboratory in orbit on an elliptic trajectory. 



a(l — e ) , cosS — e / - . 

cos^=- — — , —(E-esmE) 



l + ecosV)' l-ecos£;' V GM 

where a is the semi-major axis, e the excentricity and -ip the true anomaly. It follows that 

Aa, , . , GM GM 



-(a, -i/;) = — soao — g — '^oao — :r-ecQsil} + 0{e ) 



The first term represents the variation of the mean value of the constant on the orbit compared 
with its cosmological value. This shows that local terrestrial and Solar system experiments do 
measure the effects of the cosmological variation of the constants |123l 14571 1458) . The second term 
is a seasonal modulation and it is usually parameterized [208J as 



Ac 



= h^. (213) 

seasonal ^ 



defining the parameters ki 



6.1.2 Solar system scales 

The parameters hi can be constrained from laboratory measurements on Earth. Since e ~ 0.0167 
for the Earth orbit, the signal should have a peak-to-peak amplitude of 2GMe/ac^ ^ 3.3 x 10^^° 
on a period of 1 year. This shows that the order of magnitude of the constraints will be roughly 
of 10~^^/10~^'' ^ 10~^ since atomic clocks reach an accuracy of the order of 10^^^. The data of 
Refs. |215J and [35] lead respectively to the two constraints |208) 

ko,^^ + 0.17/ce = (-3.5 ± 6) X 10"^ \ka^^^ + O.Uke\ < 2.5 x 10"^ (214) 

for ajjjjj and TOo/Aqcd respectively. The atomic dyprosium experiment |99! allowed to set the 
constraint |194| 

fca^^j - (-8.7 ±6.6) X 10-6, (215) 
which, combined with the previous constraints, allows to conclude that 

fee = (4.9±3.9) X 10"^ fc,, = (6.6 ± 5.2) X 10^^ (216) 

for ttIo/Aqcd and Wq/AqcD respectively. Ref. [59] , using the comparison of caesium and a stron- 
tium clocks derived that 

ka -l-0.36A;e = (1.8±3.2) X 10"^ (217) 
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which, combined with measurement of H-maser [16! , allow to set the 3 constraints 

= (2.5±3.1) X 10"^ fc^ = (-1.3± 1.7) X 10~^ fc, = (-1.9 ± 2.7) x 10"'^ (218) 

Ref. [Sminn] reanalyzed the data by Ref. gOS] to conclude that +0.51fcp = (7.1±3.4) x 10"^. 
Combined with the constraint (|217p . it led to 

fc,, = (3.9±3.1) X 10-^ fc^ = (0.1± 1.4) X 10"^ (219) 

Ref. [22] also used the data of Ref. |434j to conclude 

fca^„ = (-5.4 ± 5.1) X 10"**. (220) 

All these constraints use the sensitivity coefhcients computed in Refs. 1209] . We refer to 
Ref. |264] as an unexplained seasonal variation that demonstrated the difficulty to interpret phe- 
nomena. 

Such bounds can be improved by comparing clocks on Earth and onboard of satellites |208[ 
14381 1341] while the observation of atomic spectra near the Sun can lead to an accuracy of order 
unity [208) . A space mission with atomic clocks onboard and sent to the Sun could reach an 
accuracy of lO"* [MTllMg] . 

6.1.3 Milky Way 

An attempt |319l I356J to constrain fc^ from emission lines due to ammonia in interstellar clouds of 
the Milky Way led to the conclusion that /c^ ^ 1, by considering different transitions in different 
environements. This is in contradiction with the local constraint (j218l) . This may result from rest 
frequency uncertainties or it would require that a mechanism such as chameleon is at work (see 
§ I5.4.2P in order to be compatible with local constraints. The analysis was based on an ammonia 
spectra atlas of 193 dense protostellar and prestellar cores of low masses in the Perseus molecular 
cloud, comparison of N2H+ and N2D+ in the dark cloud L183. 

A second analysis [322] using high resolution spectral observations of molecular core in lines 
of NH3, HC3N and iV2H+ with 3 radio-telescopes showed that |A/i/yLt| < 3 x 10~* between the 
cloud environement and the local laboratory environment. An offset was however measured that 
could be interpreted as a variation of of amplitude Ap./p. = (2.2 ± 0.4stat ± 0.3sys) x 10~^. A 
second analysis [323] map four molecular cores L1498, L1512, L1517, and L1400K selected from 
the previous sample in order to estimate systematic effects due to possible velocity gradients. The 
measured velocity offset, once expressed in terms of A/I, gives Afi = (26 ± Igtat ± 3sys) x 10^^. 

A similar analysis [320] based on the fine-structure transitions in atomic carbon Ci and low- 
laying rotational transitions in ^^CO probed the spatial variation of F = a^^j/i over the Galaxy. 
It concluded that 

lAi^'/i^'l < 3.7 X 10"^ (221) 

between high (terrestrial) and low (interstellar) densities of baryonic matter. Combined with the 
previous constraint on ^ it would imply that lAajj-^^/aj-Ml < 2 x 10~^. This was updated |325J to 
\AF'/F'\ < 2.3 X 10-^ so that |Aa^„/a^„| < 1.1 x lO"'^. 

Since extragalactic gas clouds have densities similar to those in the interstellar medium, these 
bounds give an upper bound on a hypothetic chameleon effect which are much below the constraints 
obtained on time variations from QSO absorption spectra. 

6.2 Cosmological scales 

During infiation, any light scalar field develop super-Hubble fiuctuations of quantum origin, with 
an almost scale invaraint power spectrum (see chapter 8 of Ref. |404j ). It follows that if the 
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fundamental constants depend on such a field, their value must fluctuate on cosmological scales 
and have a non- vanishing correlation function. More important these fluctuations can be correlated 
with the metric perturbations. 

In such a case, the fine-structure constant will behave as a^.^ = a-^^^{t) + 6a-^„{'x.,t), the 
fluctuations being a stochastic variable. As we have seen earlier, enters the dynamics of 
recombination, which would then become patchy. This has several consequences for the CMB 
anisotropics. In particular, similarly to weak gravitational lensing, it will modify the mean power 
spectra (this is a negligible effect) and induce a curl component (B mode) to the polarization [462] . 
Such spatial fluctuations also induce non-Gaussian temperature and polarization correlations in 
the CMB [462] I412j . Such correlations have not allowed to set observational constraints yet but 
they need to be included foe consistency, see e.g. the example of CMB computation in scalar- 
tensor theories ^429 . The effect on large the scale structure was also studied in Refs. [29l I360j 
and the Keck/ HIRES QSO absorption spectra showed |374) that the correlation function of the 
fine-structure constant is consistent on scales ranging between 0.2 and 13 Gpc. 

Recently, it has been claimed [47] 1520] that the fine structure constant may have a dipolar 
variation that would explain consistently the data from the Southern and Northern hemispheres 
(see § I3.4.3|) . Let assume a constant, X say, depend on the local value of a dynamical scalar field 
0. The value of X at the observation point is compared to its value here and today, 

AX/Xo = X(0)/X(0o) - 1. 

Decomposing the scalar field as (?i = 0o + A0, one gets that AX/Xq = sx{4')A(j), with sx defined 
in Eq. ()232|) . Now the scalar field can be decomposed into a background and perturbations as 
= (j){t) + 50(x, t) where the background value depends only on t because of the Copernican 
hypothesis. It follows that 

^^pii ^ sxmKt) ~ M + {sx(0) + s'xmm - m}hm = sx{^)a^ + sx{^)54>m. 

^0 

(222) 

The first term of the r.h.s. depends only on time while the second is space-time dependent. It is 
also expected that the second term in the curly brackets is negligible with respect to the first, i.e. 
Sx((j)) ^ sx{<f>)- It follows that one needs 6cf>{x,t) not to be small compared to the background 
evolution term Acf) for the spatial dependence to dominate over the large scale time dependence. 
This can be achieved for instance if is a seed field whose mean value is frozen. Because of 
statistical isotropy, and in the same way as for CMB anisotropics (see e.g. Ref. |404j ) . one can 
express the equal-time angular power spectrum of AX/Xq for two events on our past lightcone as 

AX(ni,r,0 AX(n2,r,t) \ >^ 2£+l ^^xx), , „ 

Yo Yo /=1.^^^ m(ni-n2). (223) 

If S(j> is a stochastic field characterized by its power spectrum, ((5(/)(ki, t)50(k2, t)) — P^iik, t)6(ki + 
k2) in Fourier space, then 

Cf ^)(z) = ^5i[^(z)] J P^{k,z)n[k(Tjo - v)]k^dk, (224) 

ji being a spherical Bessel function. For instance, if oc where rig is a spectral index, — 1 

corresponding to scale invariance, one gets that £{£ + l)cf^^ cx ^ on large angular scales. 
The comparison of the amplitude of the angular correlation and the isotropic (time) variation is 
model-dependent and has not yet been investigated. 

Another possibility would be that the Copernican principle is not fully statisfied, such as in 
various void models. Then the background value of (j) would depend e.g. on r and t for a spheri- 
cally symmetric spacetime (such as a Lemaitre-Tolman-Bondi spacetime). This could give rise to a 
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dipolar modulation of the constant if the observer (us) is not located at the center of the universe. 
Note however that such a cosmological dipole would also reflect itself e.g. on CMB anisotropies. 
Similar possibilities are also offered within the chameleon mechanism where the value of the scalar 
field depends on the local matter density (see ^ I5.4.2|) . 

More speculative, is the effect that such fluctuations can have during preheating after inflation 
since the decay rate of the inflaton in particles may fluctuate on large scales |292| 1293] . 

6.3 Implication for the universality of free fall 

As we have seen in the previous sections, the tests of the universality of free fall is central in 
contraining the model involving variations of the fundamental constants. 

From Eqs. ([T^ . the amplitude of the violation of the universality of free fall is given by 77^3 
which takes the form 

VAB = — y2\-fAi ~ fsil |Vai|. 
9N 

In the case in which the variation of the constants arises from the same scalar field, the analysis 
of § 16.11 implies that Va^ can be related to the gravitational potential by |VQi| ~ oiiSi{(j))acji_tgN 
so that 

VAB = ^ I/aj - /sil Si(0)azaext = ^ lA^i - Asil Sj(0)aext- (225) 

i i 

This can be expressed in terms of the sensitivity coefhcient ki defined in Eq. (|213p as 

VAB^J2\^A'-^Bz\h, (226) 

i 

since |Vai| = UikigN- This shows that each experiment will yield a constraint on a linear combi- 
nation of the coefficients ki so that one requires at least as many independent pairs of test bodies 
as the number of constants to be constrained. 

While the couplings to mass number, lepton number and the electroamgnetic binding energy 
have been considered [118] [see the example of § 15.4. 1| the coupling to quark masses remains 
a difhcult issue. In particular, the whole difficulty lies in the determination of the coefficients 
Aai [see § 15.3.2] . In the formalism developed in Refs. |120[ 1121] . see § 15.4. 1[ one can relate the 
expected deviation from the universality of free fall to the 5 parameters d and get constraints on 
Drh = d*g{drh - dg) and = d*de where d* = dg + 0.093(dm - dg) + 0.00027de. For instance, 
Be-Ti EotWash experiment and LRR experiment respectively imply 

\Dri-, + 0.22De| < 5.1 X 10"", lAn + 0.28i:»e| < 24.6 X 10-". 

This shows that while the Lunar experiment has a slightly better differential-acceleration sensi- 
tivity, the laboratory-based test is more sensitive to the dilaton coefficients because of a greater 
difference in the dilaton charges of the materials used, and of the fact that only one-third of the 
Earth mass is made of a different material. 

The link between the time variation of fundamental constants and the violation of the univer- 
sality of free fall have been discussed by Bekenstein [3H] in the framework described in § 15.4.21 and 
by Damour-Polyakov [1301 1131] in the general framework described in § 15.4.11 In all these models, 
the two effects are triggered by a scalar field. It evolves according to a Klein-Gordon equation 
(^ -t- 3iJ(/) -I- rn^cj) + . . . = 0), which implies that (f) is damped as (/) oc if its mass is much smaller 
than the Hubble scale. Thus, in order to be varying during the last Hubble time, (j) has to be very 



101 



light with typical mass m ^ Hq ^ 10^'^'^ eV. As a consequence, (j> has to be very weakly coupled 
to the standard model fields to avoid a violation of the universality of free fall. 

This link was revisited in Ref. [HSl 11651 1527) in which the dependence of on the scalar field 
responsible for its variation is expanded as 



The cosmological observation from QSO spectra implies that AA0/Af4 ~ 10"^ at best during 
the last Hubble time. Concentrating only on the electromagnetic binding energy contribution 
to the proton and of the neutron masses, it was concluded that a test body composed of 
neutrons and Up protons will be characterized by a sensitivity A(j^pi?p + i/n-Bn)/mN where (resp. 
z/p) is the ratio of neutrons (resp. protons) and where it has been assumed that rrin ^ nip ^ 
rriN- AssumineF^ that v^^^^ ^ 1/2 and using that the compactness of the Moon-Earth system 
91n(mEarth/™Moon)/c^liiaEM 10~^, ouc gcts ryi2 ~ 10~'^A^. Dvah and Zaldarriaga ^165] obtained 
the same result by considering that Au^ p ^ 6 x 10^^ — 10^^. This implies that A < 10~^ which is 
compatible with the variation of ajj,„ if A(j)/M4 > 10~^ during the last Hubble period. From the 
cosmology one can deduce that (A(^/M4)^ ^ (p^ + Pij,) / ptotai- If dominates the matter content 
of the universe, ptotai, then A0 ~ so that A ^ 10^^ whereas if it is sub-dominant A0 <C M4 
and A 10~^. In conclusion 10^^ < A < 10~^. This explicits the tuning on the parameter 
A. Indeed, an important underlying approximation is that the (^-dependence arises only from the 
electromagnetic self-energy. This analysis was extended in Ref. |143j who included explicitely the 
electron and related the violation of the universality of free fall to the variation of /i. 

In a similar analysis [527] , the scalar field is responsible for both a variation of a^^ and for the 
acceleration of the universe. Assuming its equation of state is Wh ^ 1, one can express its time 
variation (as long as it has a standard kinetic term) as 



It follows that the expected violation of the universality of free fall is related to the time variation 
of 

'^EM today by 



where Q is a parameter taking into account the influence of the mass ratios. Again, this shows 
that in the worse case in which the Oklo bound is saturated (so that d\na^„/dz ~ 10~^), one 
requires 

1 + ^ 10"^ for ri < 10-13, hence providing a string bond between the dark energy 
equation of state and the violation of the universality of free fall. This was extended in Ref. |147j 
in terms of the phenomenological model of unification presented in § 15.3.11 In the case of the string 
dilaton and runaway dilaton models, one reaches a similar conclusion [see Eq. (|202p in I5.4.1| . A 
similar result [345] was obtained in the case of pure scalar-tensor theory, relating the equation of 
state to the post-Newtonian parameters. In all these models, the link between the local constraints 
and the cosmological constraints arise from the fact that local experiments constrain the upper 
value of 00: which quantify both the deviation of its equation of state from —1 and the variation 
of the constants. It was conjectured that most realistic quintessence models suffer from such a 
problem [BSj. 

One question concerns the most sensitive probes of the equivalence principle. This was inves- 
tigated in Ref. |143j in which the coefficients XAi are estimated using the model (|188p . It was 

^^For copper Up = 0.456, for uranium Up = 0.385 and for lead Up = 0.397. 




(227) 
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concluded that they are 2-3 order of magnitude over cosmic clock bounds. However, Ref. |148j 
concluded that the most sensitive probe depends on the unification relation that exist between the 
different couplings of the standard model. Ref. [459] concluded similarly that the univerality of 
free fall is more constraining that the seasonal variations. The comparison with QSO spectra is 
more difficult since it involves the dynamics of the field between z ~ 1 and today. To finish, let us 
stress that these results may be changed significantly if a chameleon mechanism is at work. 

7 Why are the constants just so? 

The numerical values of the fundamental constants are not determined by the laws of nature in 
which they appear. One can wonder why they have the values we observe. In particular, as pointed 
by many authors (see below) , the constants of nature seem to be fine tuned [315] . Many physicists 
take this fine-tuning to be an explanandum that cries for an explanans, hence following Hoyle |256j 
who wrote that "one must at least have a modicum of curiosity about the strange dimensionless 
numbers that appears in physics." 

7.1 Universe and multiverse approaches 

Two possible lines of explanation are usually envisioned: a design or consistency hypothesis and 
an ensemble hypothesis, that are indeed not incompatible together. The first hypothesis includes 
the possibility that all the dimensionless parameters in the "final" physical theory will be fixed 
by a condition of consistency or an external cause. In the ensemble hypothesis, the universe we 
observe is only a small part of the totality of physical existence, usually called the multiverse. 
This structure needs not be fine-tuned and shall be sufficiently large and variegated so that it 
can contain as a proper part a universe like the one we observe the fine-tuning of which is then 
explained by an observation selection effect [62 . 

These two possibilities send us back to the large number hypothesis by Dirac [154] that has 
been used as an early motivation to investigate theories with varying constants. The main concern 
was the existence of some large ratios between some combinations of constants. As we have seen in 
§ 15.3.11 the running of coupling constants with energy, dimensional transmutation or relations such 
as Eq. (|184p have opened a way to a rational explanation of very small (or very large) dimensional 
numbers. This follows the ideas developped by Eddington |177[ I178j aiming at deriving the values 
of the constants from consistency relations, e.g. he proposed to link the fine-structure constant 
to some algebraic structure of spacetime. Dicke [149J pointed out another possible explanation to 
the origin of Dirac large numbers: the density of the universe is determined by its age, this age 
being related to the time needed to form galaxies, stars, heavy nuclei... This led Carter [81] to 
argue that these numerical coincidence should not be a surprise and that conventional physics and 
cosmology could have been used to predict them, at the expense of using the anthropic principle. 

The idea of such a structure called the multiverse has attracted a lot of attention in the past 
years and we refer to Ref. 78 for a more exhaustive account of this debate. While many versions 
of what such a multiverse could be, one of them finds its route in string theory. In 2000, it 
was realized |63j that vast numbers of discrete choices, called flux vacua, can be obtained in 
compactifying superstring theory. The number of possibilities is estimated to range between 10^"° 
and 10^*"^, or maybe more. No principle is yet known to fix which of these vaua is chosen. Eternal 
inflation offers a possibility to populate these vacua and to generate an infinite number of regions 
in which the parameters, initial conditions but also the laws of nature or the number of spacetime 
dimensions can vary from one universe to another, hence being completely contingent. It was 
later suggested by Susskind |476] that the anthropic principle may actually constrain our possible 
locations in this vast string landscape. This is a shift from the early hopes |269) that M-theory 
may conceivably predict all the fundamental constants uniquely. 
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Indeed such a possibility radically changes the way we approach the question of the relation of 
these parameters to the underlying fundamental theory since we now expect them to be distributed 
randomly in some range. Among this range of parameters lies a subset, that we shall call the 
anthropic range, which allow for universe to support the existence of observers. This range can 
be determined by asking ourselves how the world would change if the values of the constants 
were changed, hence doing counterfactual cosmology. This is however very restrictive since the 
mathematical form of the law of physics may be changed as well and we are restricting to a local 
analysis in the neighborhood of our observed universe. The determination of the anthropic region 
is not a prediction but just a characterisation of the sensitivity of "our" universe to a change 
of the fundamental constants ceteris paribus. Once this range is determined, one can ask the 
general question of quantifying the probability that we observe a universe as ours, hence providing 
a probabilistic prediction. This involves the use of the anthropic principle, which expresses the 
fact what we observe are not just observations but observations made by us, and requires to state 
what an observer actually is [379]. 

7.2 Fine-tunings and determination of the anthropic range 

As we have discussed in the previous sections, the outcome of many physical processes are strongly 
dependent on the value of the fundamental constants. One can always ask the scientific question 
of what would change in the world around us if the values of some constants were changed, hence 
doing some counterfactual cosmology in order to determine the range within which the universe 
would have developed complex physics and chemistry, what is usually thought to be a prerequisit 
for the emergence of complexity and life (we emphasize the difficulty of this exercice when it goes 
beyond small and local deviations from our observed universe and physics, see e.g. Ref. [244] for a 
possibly life supporting universe without weak interaction). In doing so, one should consider the 
fundamental parameters entering our physical theory but also the cosmological parameters. 

First there are several constraints that the fundamental parameters listed in Table [T] have to 
satisfy in order for the universe to allow for complex physics and chemistry. Let us stress, in a 
non-limitative way, some examples. 

• It has been noted that the stability of the proton requires rud — a^/^rrip. The anthropic 
bounds on m^, and me (or on the Higgs vev) arising from the existence of nuclei, the 
di-neutron and the di-proton cannot form a bound state, the deuterium is stable have been 
investigated in many works [51 [HI 11191 11441 11591 11601 12511 1252] , even allowing for nuclei made 
of more than 2 baryon species |263j . Typically, the existence of nuclei imposes that + 
and V cannot vary by more that 60% from their observed value in our universe. 

• If the difference of the neutron and proton masses where less that about 1 MeV, the neutron 
would become stable and hydrogen would be unstable |4361 1253] so that helium would have 
been the most abundant at the end of BBN so that the whole history of the formation 
and burning of stars would have been different. It can be deduced that [251j one needs 
TOd ~ rriu — TOc ^1-2 MeV so that the universe does not become all neutrons; md ~ rn^ + mc ^ 
3.4 MeV for the pp reaction to be exothermic and > leading to a finite domain. 

• A coincidence emerges from the existence of stars with convective and radiative envelopes, 
since it requires [79] that ~ It arises from the fact that the typical mass that 
seprates these two behaviour is roughly a~'^a^mp while the masses of star span a few 
decades aroung a~^mp. Both stars seem to be needed since only radiative stars can lead 
to supernovae, required to disseminate heavy elements, while only convective stars may 
generate winds in their early phase, which may be associated with formation of rocky planets. 
This relation while being satisfied numerically in our universe cannot be explained from 
fundamental principles. 
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• Similarly, it seems that for neutrinos to eject the envelope of a star in a supernovae explosion, 
one requires [75] ^ . 

• As we discussed in § 13.51 the production of carbon seems to imply that the relative strength 
of the nuclear to electromagnetic interaction must be tuned typically at the 0.1% level. 

Other coincidences involve also the physical properties, not only of the physical theories, but 
also of our universe, i.e. the cosmological parameters summarized in Table [H Let us remind some 
examples 

• The total density parameter 17 must lie within an order of magnitude of unity. If it were 
much larger the universe will have recollapsed rapidly, on a time scale muxh shorter that 
the main-sequence star lifetime. If it were to small, density fluctuations would have frozen 
before galaxies could form. Typically one expects 0.1 < fig < 10. Indeed, most inflationary 
scenarios lead to f2o ^ 1 so that this may not be anthropically determined but in that case 
inflation should last sufficiently long so that this could lead to a fine tuning on the parameters 
of the inflationary potential. 

• The cosmological constant was probably the first one to be questioned in an anthropical 
way |523j . Weinberg noted that if A is too large, the universe will start accelerating before 
structures had time to form. Assuming that it does not dominate the matter content of 
the universe before the redshift 2;* at which earliest galaxy are formed, one concludes that 
pv = A/SttG < (+2*)/9mato- Weinberg |523| estimated ~ 4.5 and concluded that "if it is 
the anthropic principle that accounts for the smallness of the cosmological constant, then we 
would expect the vacuum energy density ~ (10 — 100)pinato because there is no anthropic 
reason for it to be smaller" . Indeed, the observations indicate pv ^ 2pmato 

• Tegmark and Rees [479] have pointed out that the amplitude of the initial density perturba- 
tion, Q enters into the calculation and determined the anthropic region in the plane (A, Q). 
This demonstrates the importance of determinating the parameters to include in the analysis. 

• Different time scales of different origin seem to be comparable: the radiative cooling, galactic 
halo virialization, time of cosmological constant dominance, the age of the universe today. 
These coincidence were interpreted as an anthropic sign |64) . 

These are just a series of examples. For a multi-parameter study of the anthropic bound, we refer 
e.g. to Ref. |480| and to Ref. |242| for a general anthropic investigation of the standard model 
parameters. 

7.3 Anthropic predictions 

The determination of the anthropic region for a set of parameters is in no way a prediction but 
simply a characterisation of our understanding of a physical phenomenon P that we think is 
important for the emergence of observers. It reflects that, the condition C stating that the constants 
are in some interval, C P, is equivalent to \P =>!C. 

The anthropic principle j8lj states that "what we can expect to observe must be restricted 
by the conditions necessary for our presence as observers" . It has received many interpretations 
among which the weak anthropic principle stating that "we must be prepared to take account of the 
fact that our location in the universe in necessarily privilcdged to the extent of being compatible 
with our existence as observers" , which is a restriction of the Copernican principle oftenly used 
in cosmology, and the strong anthropic principle according to which "the universe (and hence the 
fundamental parameters on which it depends) must be such as to admit the creation of observers 
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within it at some stage." (see Ref. [33] for further discussions and a large bibhography on the 
subject). 

One can then try to determine the probabihty that an observer measure the value x of the 
constant X (that is a random variable fluctuating in the multiverse and the density of observers 
depend on the local value of X). According to Bayes theorem, 

PiX = a;|obs) cx P(obs|JSi: = x)PiX = x). 

P(X = x) is the prior distribution which is related to the volume of those parts of the universe in 
which X — X at dx. P{ohs\X — x) is proportional to the density of observers that are going to 
evolve when X = x. P{X — a;|obs) then gives the probability that a randomly selected observer is 
located in a region where X ~ x ±dx. It is usually rewritten as |515] 

P{x)dx = 7iobs(a;)-Ppriorda;. 

This higlights the difficulty to make a prediction. First, one has no idea of how to compute 
nohs{x)- When restricting to the cosmological constant, one can argue [515J that A does not affect 
microphysics and chemistry and then estimate riobs (x) by the fraction of matter clustered in giant 
galaxies and that can be computed from a model of structure formation. This may not be a good 
approximation when other constants are allowed to vary and it needs to be defined properly. Sec- 
ond, Pprior requires an explicit model of multiverse that would generate sub-universes with different 
values Xi (continuous or discrete) for x. A general argument [524] states that if the range over 
which X varies in the multiverse is large compared to the anthropic region X G [A'min, A'max] one 
can postulate that Pprior is flat on [ATmin , ATmax] ■ Indeed, such a statement requires a measure in 
the space of the constants (or of the theories) that are allowed to vary. This is a strong hypothesis 
which is difficult to control. In particular if Pprior peaks outside of the anthropic domain, it would 
predict that the constants should lie on the boundary of the antropic domain [437] . It also requires 
that there are sufficently enough values of Xi in the antrhopic domain, i.e. Sxi <^ Amax — Amin. 
Garriga and Vilenkin '229' stressed that the hypothesis of a flat Pprior for the cosmological constant 
may not hold in various Higgs models, and that the weight can lower the mean viable value. To 
finish, one want to consider P{x) as the probability that a random observer measures the value x. 
This relies on the fact that we are a typical observer and we are implicitely making a self sampling 
hypothesis. It requires to state in which class of observers we are supposed to be typical (and the 
final result may depend on this choice [379] ) and this hypothesis leads to conlusions such as the 
doomsday argument that have be debated actively [62] 1379] . 

This approach to the understanding of the observed values of the fundamental constants (but 
also of the initial conditions of our universe) by resorting to the actual existence of a multiverse 
populated by different "low-energy" theory of some "mother" microscopic theory allows to explain 
the observed fine-tuning by an observational selection effect. It also sets a limit to the Copernican 
principle stating that we do not leave in a particular position in space since we have to leave in a 
region of the multiverse where the constants are inside the anthropic bound. Such an approach is 
indeed not widely accepted and has been criticized in many ways (T] 1182} 11831 13971 14751 15061 1471] . 

Among the issues to be answered before such an approach becomes more rigorous, let us note: 
(1) what is the shape of the string landscape; (2) what constants should we scan. It is indeed 
important to distinguish the parameters that are actually fine-tuned in order to determine those 
that we should hope to explain in this way [5331 1534] . Here theoretical physics is indeed important 
since it should determine which of the numerical coincidences are coincidences and which are 
expected for some unification or symmetry reasons; (3) How is the landscape populated; (4) what 
is the measure to be used in order and what is the correct way to compute anthropically conditioned 
probabihties. 
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While considered as not following the standard scientific approach, this is the only existing 
window on some understanding the value of the fundamental constants. 

8 Conclusions 

The study of fundamental constants has witnessed tremendous progresses in the past years. In 
a decade, the constraints on their possible space and time varitions have flourished. They have 
reached higher precision and new systems, involving different combinations of constants and located 
at different redshifts, have been considered. This has improved our knowledge on the equivalence 
principle and allowed to test it on astrophysical and cosmological scales. We have reviewed them in 
§[3] and §[4l We have emphasized the experimental observational progresses expected in the coming 
years such as the E-ELT, radio observations, atomic clocks in space, or the use of gravitational 
waves. 

From a theoretical point of view, we have described in § [5] the high-energy models that predict 
such variation, as well as the link with the origin of the acceleration of the universe. In all these 
cases, a spacetime varying fundamental constant reflects the existence of an almost massless field 
that couples to matter. This will be at the origin of a violation of the universality of free fall 
and thus of utmost importance for our understanding of gravity and of the domain of validity of 
general relativity. Huge progresses have been made in the understanding of the coupled variation 
of different constants. While more mo del- dependent, this allows to set stronger constraints and 
eventually to open a observational window on unification mechanisms. 

To finish, we have discussed in §[7]the ideas that try to understand the value of the fundamental 
constant. While considered as borderline with respect to the standard physical approach, it reveals 
the necessity of considering a universe larger than our own, and called the multiverse. It will also 
give us a hint on our location in this structure in the sense that the anthropic principle limits the 
Copernican principle at the basis of most cosmological models. We have stressed the limitations 
of this approach and the ongoing debate on the possibility to make it predictive. 

To conclude, the puzzle about the large numbers pointed on by Dirac has led to a better 
understanding of the fundamental constants and of their roles in the laws of physics. They are now 
part of the general tests of general relativity, as well as a Breadcrumbs to understand the origin of 
the acceleration of the universe and to more speculative structures, such as a multiverse structure, 
and possibly a window on string theory. 
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A Notations 



A.l Constants 

The notations and numerical values of the constants used in this review are summarized in Table [T] 
and Tabled 

A. 2 Sensibility coefficients 

The text introduces several sensitivity coefficients. We recall their definition here. 

• Given an observable O the value of which depends on a set of primary parameters Gfc, the 
sensitivity of the measured value of O to these parameters is 

d In O , , 

= ck. (228) 



dhiGfc 

The value of the quantitoes Cfc requires a physical description of the system, 
the parameters can be related to a set of fundamental constant ai and we define 

dlnGfc 



din ai 



dn- (229) 



The computation of the coefficients d^i requires to specify the theoretical framework and 
depends heavily on our knowledge of nuclear physics and the assumptiuons on unification. 

• A particular sets of parameters dki has been sibgled out for the sensitivity of the mass of a 
body A to a variation of the fundamental constants 

—3 = JAi- (230) 

One also introduces 

^ = A., (231) 
d in Qfi 

so that 

Aai = OiifAi- 

• In models where the variation of the fundamental constants are induced by the variation of 
a scalar field with define 

d In a," , , , , 

= s^{4>)■ (232) 

• In this class of models the variation of the constants can be related to the gravitational 
potential by 

-— ^ = k^■ 233 

d$Ar 
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A. 3 Background cosmological spacetime 

We consider that the spacetime is describe by a manifold Ai with metric g^i, with signature 
(-,+,+,+). In the case of a Minkowsky spacetime g^^i/ — T]^iy • 

In the cosmological context, we will describe the universe by a Friedmann-Lemaitre spacetime 
with metric 

ds^ = -dt^ + a^{t)jij dx'dx^ (234) 

where t is the cosmique time, a the scale factor and 7^ the metric on the constant time hyper- 
surfaces. The Hubble function is defined as H = a/ a. We also define the redshift by the relation 
1 + 2 = ao/a, with gq the scale factor evaluated today. 

The evolution of the scale factor is dictated by the Fricdmann equation 

H^ = ^P-- + T, 235 
3 3 

where p = J,pi is the total energy density of the matter components in the universe. Assuming 
the species i has a constant equation of state wi = Pi/ pi, each component evolves as pi = pio(l + 
^^2(1+11);) _ rpj^g Friedmann equation can then be rewritten as 

^^Y,^^{^ + zf^'^""^^ + nK{l + zf + f^A, (236) 
with the densitiy parameters defined by 

"^--3i?f' ""^^JhI- ^'''^ 

They clearly satisfy ^ $7^ + il^ + = 1 • 

Concerning the properties of the cosmological spacetime, I follow the notations and results of 
Ref. [103. 
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